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ABSTRACT 

Steady, incompressible, turbulent, swirl-free flow through a circular-to- 
rectangular transition duct has been studied experimentally. The transition duct 
has an inlet diameter of 20.43 cm, a length- to-diameter ratio of 1.5, and an exit 
plane aspect ratio of three. The cross-sectional area remains the same at the exit 
as at the inlet, but varies through the transition section to a maximum value ap- 
proximately 15% above the inlet value. The cross-sectional geometry everywhere 
along the duct is defined by the equation of a superellipse. Mean and turbulence 
data were accumulated utilizing pressure and hot-wire instrumentation at five 
stations along the test section. Data are presented for operating bulk Reynolds 
numbers of 88,000 and 390,000. Measured quantities include total and static 
pressure, the three components of the mean velocity vector and the six compo- 
nents of the Reynolds stress tensor. The results show that the curvature of the 
transition duct walls induces a relatively strong pressure-driven crossflow that 
produce a contra-rotating vortex pair along the diverging side-walls of the duct. 
The vortex pair significantly distorts both the mean flow and turbulence field. 
Local equilibrium conditions at the duct exit, if they exist at all, are confined to 
a very small region near the wall indicating that care must be taken when using 
wall functions to predict the flowfield. Analysis of the Reynolds stress tensor 
at the exit plane shows that streamline curvature and lateral divergence effects 
distort the non-dimensional turbulence structure parameters. 

In addition to the transition duct measurements, a hot-wire technique which 
relies on the sequential use of single rotatable normal and slant-wire probes has 
been proposed. The technique is applicable for measurement of the total mean 
velocity vector and the complete Reynolds stress tensor when the primary flow is 
arbitrarily skewed relative to a plane which lies normal to the probe axis of rota- 
tion. Measurement of the mean flow has been verified in fully-developed pipe flow 
under simulated pitch angles up to ±20°. Measurements of the Reynolds stress 
tensor has been verified for the unskewed condition. Under skewed conditions, 
systematic deviations of the stress tensor were observed which are attributed to 
the geometry of the slant-wire probe. A new slant-wire probe is proposed which 
is designed to reduce, if not eliminate, the observed deficiencies. 
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CHAPTER 1 


INTRODUCTION 

The term “transition duct” refers to a class of internal flow configurations in 
which the cross-sectional shape changes in the streamwise direction. Examples 
include ducts which change cross section from square-to-round, round-to-square, 
square-to-rectangular, etc. Transition ducts are commonly used in ventilation 
systems, in aircraft propulsion systems, as wind tunnel components and as dif- 
fusers in hydro-electric turbines. High-performance military aircraft often utilize 
transition ducts as components in engine inlet and exhaust systems. A schematic 
of a typical propulsion system is shown in Fig. 1.1. The flow through these ducts 
is three-dimensional and viscous effects are not negligible. Transverse pressure 
gradients induced by streamline curvature, coupled with viscous effects, can in- 
duce streamwise vortices and regions of local flow separation which adversely af- 
fect overall performance. The propulsion system designer .must be able to predict 
the flow through a given configuration accurately in order to avoid undesirable 
flow characteristics. 

One particular transition duct that is currently of interest is the circular- 
to-rect angular (CR) transition duct. Nonaxisymmetric engine exhaust nozzles, 
usually of rectangular cross section, are being developed in order to improve the 
performance of military aircraft. An example is the Two-Dimensional/Conver- 
gent-Divergent (2-D/C-D) multi-function engine exhaust nozzle. This nozzle, 
first proposed by Pratt and Whitney Aircraft Group, provides conventional jet 
area variations as well as thrust vectoring and reversing capabilities [1]. The three 
modes of operation of the 2-D/C-D nozzle are illustrated in Fig. 1.1. To connect 
a 2-D/C-D nozzle to the engine, a circular-to- rectangular transition section is 
needed. A typical transition duct for this purpose is shown in Fig. 1.2. The ideal 
transition duct should provide a uniform, subsonic flow over as short a distance 
as possible, without incurring unacceptable pressure losses. 

Another application of the CR transition duct is as a component in hydro- 
electric turbines. Hydraulic draft tubes, which are CR transition ducts with a 
90° bend, are used to diffuse the water exiting the vertical turbine runner at 
hydro-electric power plants. The inlet flow to the draft tubes typically contains 
a significant degree of swirl. Pressure recovery in the draft tube is the primary 
parameter of interest for this application. 

Developments in computational fluid mechanics have made it possible to cal- 
culate the three-dimensional turbulent flows generated in these transition ducts. 
A duct which changes from a circular to a rectangular configuration, while main- 
taining a nearly constant cross-sectional area, offers a particularly stringent test 
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for code validation purposes because one set of opposite walls converges while 
the other set diverges at a given streamwise location. This behavior leads to 
regions of convex and concave streamline curvature in the flow which tend to 
respectively suppress and amplify turbulence intensities. In addition, transverse 
pressure gradients of opposite sign are generated by the curvature which induce 
secondary flow that can lead to streamwise vortex formation. If swirl is imparted 
to the flow, then additional helical streamline curvature will be present which 
can have a significant effect on the local turbulence structure. 

An ongoing program at NASA-Lewis Research Center is the sponsorship 
of benchmark quality experimental studies for the purpose of validating three- 
dimensional viscous flow solvers. A variety of inlet and exhaust flow components, 
including the circular-to-rectangular transitio n d uct, are being investigated. The 
Computational Fluid Dynamics Branch at NASA-Lewis has designated a series 
of circular-to-rectangular transition ducts of specific geometry to be of inter- 
est. These ducts cover a wide range of aspect ratios (AR), length ratios ( L/D ) 
and cross-sectional area variations. To facilitate mesh generation for numeri- 
cal computations, the cross sections of all these ducts are defined by superel- 
lipses (see Appendix A). A recent study by Patrick and McCormick [2,3], which 
was under sponsorship by NASA-Lewis, has provided a limited experimental 
database for two of the these ducts, the AR310 ( AR = 3, LID = 1.5) and AR630 
( AR = 6, L/D — 3.0) transition ducts. The particular duct configuration cho- 
sen for the present study has an aspect ratio of three and a length- to- diameter 
ratio of 1.5. An isometric view of the duct is shown in Fig. 1.2. At any given 
streamwise location, the cross-sectional configuration is defined by the equation 
of a superellipse. The local cross-sectional area ratio increases from unity at the 
inlet to a maximum of 1.15 at the midpoint before decreasing back to unity at 
the duct exit. This distribution is intended to model the area variation of a duct 
constructed of flat surfaces and conical sections more typical of the manufactured 
product. The exit plane cross-sectioned shape does not have sharp corners but a 
variable radius fillet as seen in Fig. 1.2. A flow visualization study by Reichert 
et al. [4] shows that the flow remains wholly attached within this configuration 
for subsonic flow in the presence of the adverse pressure gradient induced by the 
15% area expansion. The complete geometric description of this duct is given in 
Appendix A. 



Fig. 1.1. Schematic of typical propulsion system. 
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a) Isometric view. 

Fig. 1.2. Circular-to-rectangular transition duct (Continued . 





CHAPTER 2 


PREVIOUS WORK 

Turbulent flow through a circular- to-rectangular transition duct is character- 
ized by streamline cufvature and streamwise vorticity embedded in the boundary 
layer. In addition, the boundary layer is subjected to both lateral convergence 
and divergence. This chapter will begin by reviewing previous transition duct 
studies, followed by a discussion of the effect that streamline curvature, embed- 
ded vorticity and lateral divergence have on turbulent boundary layer flow. 

2.1 Flow Through Transition Ducts 

The earliest study on transition duct flow was an experimental investiga- 
tion done by Mayer [5] in 1939. In his study, he investigated flow through 
two rectangular-to-circular (and vice versa) transition ducts of constant cross- 
sectional area. The ducts had transition lengths of 0.69 and 2.76 hydraulic diam- 
eters. The data included streamwise static pressure distributions, total pressure 
contours and the three-dimensional velocity field. In a similar study, Taylor et al. 
[6] investigated turbulent flow through a square-to-round transition duct which 
had a 21.5% reduction in cross-sectional area over the length of the transition. 
This decrease in area was the result of the hydraulic diameter (40 mm) being 
held constant along the duct. The transition occurred over two hydraulic diam- 
eters. LDV techniques were used to measure streamwise and transverse velocity 
components along the duct for an operating Reynolds number of 35,350. The 
results of these studies have shown that the length of the transition section is 
influential on flow development and that pressure-driven crossflows (~ 10% of 
maximum streamwise velocity), can lead to significant distortion of the primary 
flow. 

During the early development stages of the 2-D/C-D nozzle, Pratt and Whit- 
ney Aircraft developed a design procedure for circular-to-rectangular transition 
ducts intended to minimize both pressure losses and axial length. This procedure 
is governed by the following criteria [1]: 

1) Constant cross-sectional area 

2) Corner radius decreasing linearly with length 

3) Straight sidewalls 

4) Sidewall divergence angle limit = 45 degrees 

In a recent combined experimental and numerical study by Burley et al. [7,8], 
the original PWA guidelines were examined to determine if expanded design 
criteria could be established so that shorter transitions are possible, thereby 
reducing exhaust system weight. Five different circular-to-rectangular transition 
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duct configurations were investigated to explore the effects of duct length, wall 
shape and cross-sectional area distribution on performance. All of the ducts 
were defined by super-elliptic cross sections (see Appendix A). The transition 
ducts were installed in a transonic wind tunnel with a high aspect ratio, non- 
axisymmetric nozzle and the overall internal performance was measured. In 
addition, one duct was tested with swirl vanes installed. Discharge coefficient 
and thrust ratio versus nozzle pressure ratio were used as performance criteria. 
The results of their investigation show that for length ratios less than or equal 
0.75, large regions of separated flow are present. However, because the flow 
reattached before the entrance to the nozzle, only a small decrease in performance 
was observed. They also found that swirling the flow had a positive effect on 
performance for low nozzle pressure ratios, but that performance was decreased 
when the nozzle was near a choked condition. Finally, these researchers reported 
that decreasing the cross-sectional area along the duct reduced flow separation 
and provided a modest increase in performance. 

Patrick and McCormick [2,3] were the first to make turbulence measure- 
ments within a circular-to-rectangular transition duct. LDV and total pressure 
measurements were made at the inlet and outlet stations of two different ducts 
at an operating Reynolds number of 420,000. The first duct, designated the 
AR310, had an aspect ratio of three, a length-to-diameter ratio of one, and con- 
stant cross-sectional area through the duct. The second duct, designated the 
AR630, had an aspect ratio of six and a length-to-diameter ratio of three. The 
local cross-sectional area ratio increased from unity at the inlet to a maximum 
of 1.10 at the midpoint before decreasing back to unity at the duct exit. In order 
to facilitate grid generation for numerical comparisons, the cross-sectional shape 
everywhere along the ducts was prescribed by the equation of a superellipse. 
Measured quantities included all three mean velocity components and the three 
normal Reynolds stress components at the inlet and outlet planes. The results 
for the AR310 duct showed that the axial mean flow did not develop uniformly 
but had a convex profile along the major axis at the duct exit plane. Outward 
transverse velocities, nominally parallel to the major axis, were observed that 
peaked at about 10% of the bulk velocity. No streamwise vorticity was observed 
except deep in the corner region, but the measurement grid was too coarse to 
discern discrete vortical motion. The AR630 duct behaved quite differently. Here 
the flow developed much more uniformly, and a pair of discrete vortices along 
the duct sidewalls, centered about the duct semi-major axis, were observed. The 
origin of these vortices is in the first half of the transition where the wall cur- 
vature creates a pressure gradient which causes a crossflow from the upper and 
lower walls to the sidewalls. The crossflow meets at the duct centerline and turns 
inward along it. In the second half of the duct, where the curvature changes sign, 
the pressure gradient is reversed, counteracting the secondary motion. If a vorti- 
cal pattern was established in the first half of the AR310 duct, then the reversed 
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pressure gradient was effective in stopping it. 

Miau ti al. [9] experimentally investigated three CR ducts with length-to- 
diameter ratios of 1.08, 0.92 and 0.54, under low subsonic flow conditions. The 
aspect ratio was equal to two and the cross-sectional area was constant for all 
three ducts. Mean flow and turbulence data were taken at the inlet and exit 
planes. Secondary flow patterns indicative of streamwise vortex formation were 
observed at the exit plane of the ducts. From these results, all the terms in 
the axial mean vorticity equation were computed. Their analysis showed that 
the generation of streamwise vorticity is due primarily to transverse pressure 
gradients induced by geometrical deformation. 

With the exception of the performance data reported by Burley et al. [7,8], 
none of the above studies considered the case where a swirl velocity component 
is imparted to the inlet flow. The addition of swirl may have several benefits. 
First, swirl will impart a radial velocity component to the flow, thus improving 
it’s ability to follow steeply sloped sidewalls in the transition duct. Secondly, 
Schwartz [10] has observed that noise associated with axisymmetric jet exhaust 
can be reduced by swirling the flow. Finally, in an axisymmetric jet, the rate 
of decay of the axial velocity component can be substantially increased (reduced 
thermal plume) by swirling the flow, with minimal loss of thrust [11]. Der et al. 
[12] performed a water tunnel flow visualization study of swirling flow through a 
CR duct. Later, Chu et al. [13,14] analyzed these data and found that swirling 
the flow dramatically reduced the thermal plume. Recently, Reichert et al. [4], 
using a duct identical to the one in the present study, compared the mean flow 
field for the cases of swirling and non-swirling inlet flow at an operating Mach 
number of 0.35 and a Reynolds number of 1.5 x 10®. 

Related studies have been undertaken which include the effects of a turbine 
centerbody and axial centerline curvature. Sobota and Marble [15] performed a 
detailed experimental and numerical investigation of a CR transition duct with a 
large centerbody and various degrees of inlet swirl. This study provided insight 
into vorticity generation mechanisms and the ability to tailor vorticity distribu- 
tions. CR transition ducts with a 90° bend in the transition section (hydraulic 
turbine draft tubes) have recently been analyzed with and without swirl by Vu 
and Shyy [16]. The flow through the draft tube was predicted using a finite- 
volume approximation to the full Navier-Stokes equations in conjunction with a 
k-e turbulence model and the results were compared to wind and water tunnel 
experimental data. In general, pressure recovery, as well as the three-dimensional 
velocity field, agreed well between experiment and numerical predictions. 

The review of the literature has revealed that there is only a small experi- 
mental database for flow through a circular-to- rectangular transition duct that is 
of sufficient detail to be useful for CFD code calibration/validation purposes. At 
the present, only Miau’s data set include measurements of the complete Reynolds 
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stress tensor in a CR transition duct. The present study is intended to help fill 
this void by providing complete mean flow and Reynolds stress measurements at 
the inlet and outlet stations, supplemented by mean flow data at intermediate 
stations, for a duct with an aspect, ratio larger than that considered by Miau. 
On the basis of previous work in this area, it was anticipated that skew induced 
secondary flow would have a dominating influence on the primary flow and on 
the local turbulence structure. 

2.2 Streamline Curvature Effects 

A great deal of work has been published on the effects of streamline curvature 
on turbulent boundary layer development. The bulk of the studies have been for 
the quasi-two-dimensional case with the curvature induced by a constant radius 
bend in a square or rectangular wind tunnel. Although the transition duct flow 
is considerably more complex, it is useful to examine previous related results in 
order to gain some insight into the mechanisms operating within the transition 
duct. 

In 1973, Bradshaw [17] published a comprehensive review of the effects of 
streamline curvature on turbulent flow. His work was motivated by what he 
referred to as “the surprisingly large effect exerted on shear-flow turbulence by 
curvature of the streamlines in the plane of the mean shear”. Flows with stream- 
line curvature are characterized by the presence of extra rates of strain, that is, 
rates additional to the simple shear dUfdy. When the equations of motion are 
written in semi-curvilinear coordinates ( e.g the s,n system of reference [17]), 
extra explicit terms appear which account for the presence of curvature. Ex- 
perimental measurements have shown, however, that the effects of extra rates 
of strain are an order of magnitude larger than would appear when calculation 
methods for simple shear flows are extended to curved flows. Bradshaw explains 
this discrepancy by concluding that streamline curvature directly causes large 
changes in the higher-order parameters of the turbulence structure. 

Convex and concave curvature are often referred to as stabilizing and desta- 
bilizing curvature, respectively. Laminar flow over a destabilizing (concave) sur- 
face is subject to centrifugal instability which is characterized by the presence 
of streamwise vortices within the boundary layer; the so-called Taylor-Gortler 
vortices. For turbulent flows over a concave surface, the presence of vortices 
analogous to the laminar Taylor-Gortler type have been observed experimentally 
by So and Mellor [18,19], Meroney and Bradshaw [20] and Hoffmann et al. [21], 
as well as others; and numerically, by direct simulation of the Navier-Stokes 
equations, by Moser and Moin [22]. In addition, the stabilizing and destabilizing 
effect on turbulence acts, respectively, to attenuate and amplify the turbulence 
intensities. 
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In regions sufficiently close to curved walls, mean velocity profiles have been 
observed to follow the flat plate law-of-the-wall for both convex and concave 
curvature [18,19,23]. Hoffman and Bradshaw [24] have suggested that the law- 
of-the-wall applies when y / R c is small. This is important in that it allows for the 
use of law-of-the-wall based wall functions in numerical computations. The agree- 
ment with law-of-the-wall behavior is apparently where the similarity between 
convex and concave curved flows end. In contrast to flat plate flows, turbulent 
flow over a convex surface is characterized by slower boundary layer growth, lower 
wall shear stress, reduced turbulence intensities and reduced heat transfer rates. 
Conversely, turbulent flow over a concave surface is characterized by increased 
boundary layer growth, higher wall shear stress, increased turbulence intensities 
and increased heat transfer rates, as well as the aforementioned streamwise vor- 
tices. The most significant difference in the turbulence statistics appears in the 
Reynolds shear stress. Measurements by So and Mellor [18], Gillis and Johnston 
[25] and Smits et al. [26] show a sharp decrease in the turbulent shear stress 
to near-zero levels in the outer region of the boundary layer when a flat plate 
boundary layer is suddenly subjected to a strong ( 6/R c « 0.10) convex curvature. 
For the concave case, the turbulent shear stress dramatically increased to a near 
two-fold level as compared to flat plate results. Hunt and Joubert [27] found 
similar behavior in flows subjected to mild streamline curvature (S/R c « 0.01), 
but to a lesser degree. The two flow cases also respond differently when they are 
subjected to a flat plate recovery region. Whereas the Reynolds stresses recover- 
ing from convex curvature do so in a monotonic fashion, the stresses in concave 
flows drop well below their entry region values before recovering [26]. 

The dramatic differences between the convex and concave curvature cases 
have hindered development of adequate turbulence models because the mecha- 
nisms that produce them are not well understood. Indeed, Muck tt al. [28] have 
concluded that, although governed by the same dimensional analysis, there is 
no other useful connection between the two cases. FVom this conclusion, they 
imply that allowances for the effect of streamline curvature in calculation meth- 
ods for turbulent flows should be formulated separately for the stabilizing and 
destabilizing cases. 

The above discussion serves to illustrate factors which may add to the com- 
plexity of transition duct flow. The degree to which streamline curvature affects 
transition duct flow depends on the thickness of the incoming boundary layer. 
For the present study, a boundary layer thickness of 6/ R ^ 0.25 is anticipated. 
This corresponds to a maximum curvature parameter of 6/R c « 0.085 for both 
convex and concave walls. On the basis of previous experimental results, it was 
expected that streamline curvature would influence the development of the flow 
in the transition duct. 
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2.3 Embedded Streamwise Vortices 

The transition duct studies reported in Refs. 2,3 and 9 have shown that 
streamwise vorticity is generated within CR transition ducts. The results for 
the AR630 duct (Refs. 2 and 3) s how that a discr ete vor tex pair (common flow 
away from the surface) develops in the duct sidewall boundary layer. Stream- 
line vortices embedded in a boundary layer are as often beneficial as they are 
detrimental. On aircraft flight surfaces, s treamw ise vortices are purposefully 
generated to promote mixing between the freestream and the boundary layer in 
order to forestall flow separation. In combustor applica tions, streamwise vor- 
tices are used to enhance mixing between fuel and Oxidant. In turbomachinery, 
however, streamwise vortices generated by blade-hub junctures may sweep away 
the protective film cooling on adjacent blades causing damaging hot spots. In 
transition duct applications, streamwise vortices result in undesirable pressure 
losses, although they may inhibit flow separation. 

Streamwise vortices in boundary layers can be generated by transverse pres- 
sure gradients or by gradients of the Reynolds stresses. Pressure-gradient induced 
streamwise vorticity occurs whenever a shear layer (laminar or turbulent) with 
spanwise vorticity is deflected laterally by transverse pressure forces. If the de- 
flection occurs over a short spanwise distance, then a discrete vortex is formed. 
Vortices generated in this manner occur in strut-endwall (junction) configura- 
tions and in flow through curved ducts. Reynolds stress induced vorticity occurs 
in turbulent flow through non-circulax ducts, even when the ducts are straight, 
i.e., uncurved in the streamwise direction. Because both types of vorticity gener- 
ation occur in many practical engineering flows, a large body of literature exits 
on the subject. Some of the more comprehensive experimental studies on the 
effect that embedded streamwise vortices have on the mean flow and turbulence 
structure include those due to Shabaka et al. [29], Mehta and Bradshaw [30] 
and Pauley and Eaton [31,32]. These researchers studied the effects of single 
and paired vortices in an otherwise two-dimensional boundary layer flow. Mean 
flow and turbulence (one-point double and triple correlations) were measured. 
The vortices in all these studies were generated by half-delta wings, although 
the placement of the generators in the wind tunnels differed. Whereas Pauley 
and Eaton placed the generators on the floor of the 2-D channel, Shabaka and 
Mehta placed the generators upstream of the contraction in the settling chamber. 
By placing the generators in the plenum, the velocity deficit in the wake of the 
delta wing is reduced to a small percentage of the freestream velocity as the flow 
accelerates through the contraction. Pauley and Eaton obtained data in measure- 
ment planes 97 and 188 cm downstream from the generators. Shabaka and Mehta 
present results in measurement planes between 60 and 255 cm downstream. The 
results of these investigations showed that thickening of the boundary layer oc- 
curred in upwash regions and thinning occurred in downwash regions. Paired 
vortices with the common flow away from the surface were attracted and moved 
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away from the surface. In contrast, vortex pairs with the coininon flow towards 
the surface moved away from each other and stayed in close proximity to the 
wall. In the vicinity of the vortex core(s), a concentrated maxima of turbulence 
intensity occurs, but no large-scale unsteadiness in the flow was detected. ^Vhen 
compared to the surrounding 2-D boundary layer, large changes in the dimen- 
sionless turbulence structure parameters were observed and eddy viscosities were 
reported to be very ill-behaved in the vortex region. These observations led the 
researchers to conclude that full Reynolds stress transport modelling would be 
required for prediction purposes. 

Liandrat et al. [33] used the data of Refs. 29 and 30 for comparison with 
numerical simulations based on mixing length and k-c turbulence models. In 
addition, calculations based on two forms of the Reynolds stress transport equa- 
tions were performed. The results of this investigation showed simple turbulence 
models provide good estimations of overall mean flow properties for the case of a 
single embedded vortex. For the case of paired vortices with common flow away 
from the surface, the mean flow results were found to be largely unsatisfactory. 
For both cases, details of the predicted turbulence structure required Reynolds 
stress transport models. However, even these higher-order models did not give 
adequate predictions of the transverse normal stresses and the secondary shear 
stresses that control the diffusion of streamwise vorticity. In a review of turbu- 
lent secondary flows, Bradshaw [34] concludes that the primary inadequacy in 
the Reynolds stress transport models is in the modelling of the pressure- strain 
term. 

As a final note, Patrick and McCormick [2,3] made mention of the similarity 
between the generation of the vortex pair in the AR630 CR transition duct and 
in a circular pipe with an S-shaped bend. Experimental mean flow results and a 
discussion of vorticity generation for the latter case is presented by Bansod and 
Bradshaw [35]. Limited turbulence measurements in a circular S-shaped duct 
with embedded vortices are presented by Taylor et al. [36]. 

2.4 Laterally Diverging Boundary Layer 

Along the converging walls of a CR transition duct, the boundary layer is 
subjected to lateral divergence. Conversely, along the diverging walls, the bound- 
ary layer converges laterally. Boundary layer divergence is another example of a 
shear flow with extra strain rates (dV/dy, dW/dz). As in the case of streamline 
curvature, the effect that lateral divergence has on the turbulence structure of a 
boundary layer is much larger than, and sometimes opposes, what is predicted 
by explicit terms that appear in the Reynolds stress transport equations. The 
strength of divergence in a boundary layer is typically characterized by the rate- 
of-strain parameter (dW/dz) /(dU/dy). Smits et al. [37] reviewed the effect that 
lateral divergence and convergence have on boundary layer flow. In addition, 
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they obtained mean flow and turbulence measurements in a diverging boundary 
layer that develops on a cylinder-flare where (( dW/dz)/(dU /dy ) 0.1) midway 

through the layer. In the transition region between the cylinder and flare, the 
flow was subjected to relatively strong concave longitudinal curvature. Although 
they found that the effects of divergence and curvature in the transition region 
could not be quantitatively separated, they argue that the memory of the cur- 
vature is short-lived and that the downstream flowfield is primarily a result of 
divergence effects alone. They note that boundary layers subjected to lateral 
divergence and convergence tend to, respectively, thin and thicken. Turbulence 
kinetic energy is amplified for the diverging case and attenuated for the converg- 
ing case. The shear stress in a diverging boundary layer is elevated and a peak 
occurs which moves outward from the surface as the flow develops. Pauley and 
Eaton [31,32] obtained mean flow and turbulence measurements in the diverging 
boundary layer that develops between an embedded vortex pair with the com- 
mon :flow towards the surface. Although the strength of divergence was fairly 
weak (( dW/dz)/(dU/dy ) ~ 10 — 3 ), the flow is unique, inasmuch as complicating 
factors such as streamwise pressure gradients and curvature are absent. Unlike 
the results of Smits (and others), these researchers found no significant differ- 
ence in either the mean flow or the Reynolds stresses when compared to their 
counterparts in the two-dimensional boundary layer outside the vortex pair. 

2.5 Concluding Remarks 

Turbulent flow through a circular-to-rectangular transition duct represents 
a practical engineering flow of interest where multiple complicating effects are 
present. In particular, the transition duct geometry imposes extra rates-of-strain 
on an initially two-dimensional boundary layer possessing only the 1 simple strain 
dU/dy. Although Smits et al. [37] have shown that individual effects cannot be 
quantitatively separated due to non-linear interactions, qualitative assessment 
of the flowfield should certainly be possible. By studying this flow two goals 
are hoped to be achieved. First, a sufficiently detailed mean flow and turbulence 
data set will be provided that may be used for direct comparison with numerically 
generated results. Secondly, the results will be analyzed in a way that will aid 
predictors and modellers in determining the level of sophistication required in 
their computational efforts. 



CHAPTER 3 


EXPERIMENTAL PROGRAM 


3.1 Introduction 

The present study is primarily of an experimental nature. The goal of the 
study is to provide a comprehensive set of mean and turbulence measurements 
at inlet, intermediate and outlet stations of a CR transition duct. These data 
are intended for use in CFD code calibration/ validation and turbulence model 
development. The wind tunnel flow facility, test section instrumentation and 
some data reduction methods are described in this chapter. 

3.2 Flow Facility 

The Square Duct Flow Facility in the Heat Power Laboratory in the Mechan- 
ical Engineering Building has been modified to a configuration appropriate for 
the present study. The inlet to the Square Duct Flow Facility has been replaced 
by the new ductwork illustrated in Fig. 3.1. Atmospheric air enters the wind 
tunnel through a 2.8:1 elliptic bellmouth contraction and then passes through a 
settling chamber which consists of an Alfco combination filter/honeycomb flow 
straightener (Model 36” 0 CYL,9FG-B), two fine (20 x 20) mesh screens, and a 
20:1 concentric contraction. This inlet section was designed to provide a uniform, 
low turbulence level, axisymmetric exit flow using design criteria developed by 
Morel [38]. Prior to construction, the flow through the 20:1 contraction was 
computed using NASA-Lewis ’ VISTA program [39], which is an axisymmetric 
subsonic Navier-Stokes flow solver. The inlet boundary layer thickness was var- 
ied from 6i/R a = 0.001 to 0.05, where R a is the settling chamber radius, and 
the inlet velocity was varied from {7, = 0.762 to 1.524 m/s, which corresponds to 
velocities at the exit plane from U e = 15.2 to 30.5 m/s. For all cases, the results 
showed that the inviscid core velocity at the exit plane of the contraction was 
uniform to within 0.1% of the centerline value. 

From the 20:1 contraction, the flow enters a 20.42 cm diameter pipe of 
variable length. Pipe sections are made in lengths of LjD = 3 and can be added 
or removed so as to vary the boundary layer thickness (6) at the inlet to the 
transition duct. To promote transition to turbulent flow, a 2.54 cm wide strip 
of #36 sandpaper was placed at the beginning of the first pipe. At the end of 
the pipe section is a probe access ring which has provisions for making detailed 
measurements of the test section inlet flow. From the probe access ring the flow 
enters the test section which consists of the transition duct and a removable 
transition duct extension. The rectangular duct downstream of the test section 
supports a probe traversing mechanism. Finally, a diffuser section (2 degrees 
divergence) provides the link to the existing 0.254 x 0.254 meter square duct. 
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Air is drawn through the facility by means of a two-speed centrifugal fan 
located at the exit of the Square Duct Flow Facility. The fan discharges the 
air back into the laboratory through a set of remotely actuated shutters which 
provide a means for varying the mass flow rate through the wind tunnel. Details 
of the remaining Square Duct Flow Facility are described by Eppich [40]. 

The primary materials used to build the flow facility are as follows. The 
two inlet contractions are constructed of polyester resin fiberglass with embed- 
ded aluminum mounting flanges. The settling chamber, which houses the fil- 
ter/honeycomb and screens, is constructed of wood and Formica. The 0.203 
meter diameter inlet pipes, the probe access ring, the probe access duct and the 
diffuser are all fabricated of aluminum. The transition duct and the transition 
duct extension are constructed of epoxy resin fiberglass in halves which part in 
the x-y plane. The patterns for molding the duct halves were provided by the 
NASA-Lewis Research Center. 

3.3 Test Section 

A side view of the test section is shown in Fig. 3.2.a. Cross-sectional views 
at each of the six data stations indicated in Fig. 3.2.a are shown in Figs. 3.2.b 
through 3.2.f. The origin of the coordinate system shown in Fig. 3. 2. a. was 
chosen so as to be in agreement with NASA-Lewis’ definition of the transition 
duct. Station 1 is located one inlet duct diameter upstream from where the 
beginning of transition occurs (Station 2). Although the cross section at Station 
2 is still circular, it was anticipated that some distortion of the flow field will 
occur due to the influence of the changing downstream geometry. At Stations 3 
and 4 the geometry is changing rapidly and relatively large transverse velocities 
were expected. Station 5 is located at the end of transition and data Station 6 
is located two diameters downstream from the end of transition. An isometric 
view of the six data station cross sections is shown in Fig. 3.3. 

The probe traversing mechanism located just downstream of the test section 
in Fig. 3. 2. a was used for acquiring data at Stations 5 and 6. This mechanism 
holds the probe axis parallel to the duct centerline and has provisions for rotating 
the probe about it’s longitudinal axis by means of a spring-loaded (anti-backlash) 
bevel gear arrangement. Dial indicators were used to position the probes in each 
direction to within an estimated accuracy of ±0.025 mm, and the probes were 
rotated to fixed angular positions to within an estimated accuracy of ±0.5 degree. 
The vertical and horizontal traversing capabilities are such that the probe can 
be positioned anywhere within the cross section at Stations 5 and 6. However, 
due to the convergence of the upper and lower walls of the transition duct, there 
is an increasingly larger area upstream from Station 5 which cannot be accessed. 
Therefore, at Stations 1 through 4 probes were inserted through holes in the 
duct wall which are normal to the duct axial centerline. These access holes are 
illustrated in Figs. 3.2.b through 3.2.e. 
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For this study, two identical transition ducts were constructed. The first of 
these was used when taking data at Stations 1 , 5 and 6. This duct has only static 
pressure taps along the periphery of the lower half of the duct at Station 5 (see 
Fig. 3.2.f). The second duct has probe access holes along the periphery of the 
upper half of the duct and static pressure taps along the periphery of the lower 
half of the duct at Stations 2, 3 and 4 (see Figs. 3.2.c through 3.2.e). Two ducts 
were constructed to insure that there would be no influence on the flow field at 
Stations 5 and 6 from probe access holes at upstream stations. 

3.4 Instrumentation 

The test section is instrumented to facilitate measurement of both mean 
and fluctuating quantities. Mean quantities were measured by means of static 
and total pressure instrumentation and hot-wire anemometry. All fluctuating 
quantities were measured using hot-wire anemometry. Local skin friction was 
measured with Preston tubes, which are simply circular Pitot tubes resting on 
the duct wall. 

A variety of pressure probes were used to measure the various mean quanti- 
ties of interest. Total pressure contours were measured with two types of probes. 
At Stations 1,2,5 and 6, where streamlines in the cross section are everywhere 
nominally parallel to the axial centerline of the duct, circular Pitot tubes having 
an outside tip diameter of 0.635 mm were used. At intermediate data Stations 
3 and 4, where the streamlines were skewed by as much as 20 degrees relative 
to the probe centerline, a United Sensor Model KAC12 Kiel probe was used. 
According to Chue [41], Kiel probes are able to measure total pressure accu- 
rately for skew angles as high as 40°. Boundary layer profiles were measured by 
means of flattened Pitot tubes having outer dimensions of 0.812 x 0.406 mm at 
the tip. Static pressure distributions on the duct axial centerline were measured 
by means of a static pressure probe which was traversed along the centerline. 
Transverse flow angles were measured by means of a two-tube Conrad probe and 
a normal hot-wire. All pressure data were measured with a 10 torr Barocel elec- 
tronic manometer, Model 571D-10T-1C2-V1, coupled with a Datametrics digital 
display unit, Model 1174. 

Hot-wire probes consisted of single, rotatable normal and slant-wires. The 
sensing element of the hot-wire probes is 0.00381 mm (4/x) diameter platinum 
coated tungsten wire. The ends of the wire are copper plated to facilitate sol- 
dering and to define the sensing element length, which typically has a length-to- 
diameter ratio near 300. The prongs of the probes are similar to configurations 
recommended by Comte- Bellot et al. [42], for minimizing aerodynamic dis- 
turbances. The fluctuating turbulence signal was processed by means of a TSI 
Intelligent Flow Analyzer (IFA) system, which consists of an IFA 100 constant 
temperature anemometer and an IFA 200 analog-to-digital (A/D) converter. The 
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digitized signal is recorded on a VAX Workstation, via a DRQ3B DMA card, 
where the various turbulence correlations are computed. Appendix B contains a 
description of the setup and operation of this equipment. 

The probe access ring (see Fig. 3. 2. a) was used to measure the flow conditions 
at the transition duct inlet (Station 1). The ring contains four probe access 
holes, four static pressure taps and four Preston tubes equally spaced around the 
periphery of the ring as shown in Fig. 3.2.b. To ensure that Preston tube data 
were taken within the law-of-the-wall region of the boundary layer, each Preston 
tube was of a different diameter; 2.769, 3.962, 5.537 and 6.350 mm. The tubes 
were semi-permanently installed and were removed before data were taken at the 
downstream stations. 

3.5 Data Reduction 

This section contains the data reduction methods for the pressure probes, 
hot-wire probes and equations for computing boundary layer parameters. 

3.5.1 Pressure Probe Data Reduction 

3.5. 1.1 Mean Velocity 

The total mean velocity along a flow streamline can be deduced from a 
Pitot probe aligned with the streamline or, if alignment is not practical, from 
measurements with a Kiel probe. The total velocity is related to the probe 
pressure by Bernoulli’s relation and the ideal gas law, namely: 

U 0 = (2 hR air T arnb /P amb y' 2 (3.1) 


where: 


U Q 

= total velocity (m/s) 


h 

= measured pressure head (P t — P ) 

(mmHg) 

Pair 

= gas constant for air (287.0 J/(kg ■ 

°K)) 

Tamb 

= ambient temperature (° K) 


P a mb 

= ambient static pressure ( mmHg ) 



The Cartesian velocity components are related to the total velocity vector 

by: 

U = Uq cos ft cos 7 

V = C/osin/?cos7 (3-2) 

W = U 0 cos ft sin 7 

where ft and 7 are the flow angles in the x-y and x-z planes, respectively. These 
flow angles can be measured directly by performing a differential pressure nulling 
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technique with a two-tube Conrad probe, first in the x-y plane and then in the 
x-z plane. 

3.5. 1.2 Skin Friction 

Head and Vasanta Rain’s [43] tabulated presentation of Patel’s [44] Preston 
tube calibration was used to deduce local skin friction values. Their table presents 
the calibration in the functional form: 

— = /(^) (3-3) 

Tw PV 

where A p is the difference between the Preston tube and local wall static pressure 
and d is the Preston tube outside diameter. Head and Vasanta Ram estimate 
that, even without interpolation, their tables should give values that are accurate 
to within ±1 per cent. For the present study, a linear interpolation was used. 
This method of evaluating shear stress presumes that the two-dimensional form 
of the law-of-the-wall is valid and that streamwise pressure gradients are small. 
It was anticipated that this method would be applicable to data taken at Stations 
1,5 and 6, inasmuch as there is no longitudinal wall curvature or cross-sectional 
area change at these locations. 

3.5.2 Hot-Wire Probe Data Reduction 

For reasons given in section 3.3, the positioning of the hot-wire probe is 
dependent upon the particular station being investigated. At Stations 1 through 
4, the probe body centerline was positioned in a direction normal to the axial 
centerline of the duct, while at Stations 5 and 6 the probe was positioned parallel 
to the axial direction. To simplify the following discussion, the former positioning 
method shall hereafter be designated Method A , while the latter will be referred 
to as Method B 

For the hot-wire measurements, single-wire techniques will be employed 
rather than the more complicated two or three-wire methods. This approach 
eliminates some of the difficulties associated with multi-wire techniques, such 
as extensive calibration requirements, possible wire interference, multi-wire drift 
and poor spatial resolution. Single-wire techniques, however, limit turbulence 
measurements to second-order correlations (Reynolds stresses). 

For Method B , a hot-wire technique developed by Al-Beirutty [45,46] and 
Arterberry [47] was used to relate the three mean velocity components and six 
Reynolds stresses to the mean and mean-square anemometer output voltages. 
This technique uses a fixed normal-wire and a single, rotatable slant-wire and 
utilizes an empirical cooling velocity law. The method is applicable to flows of 
low-to-moderate turbulence intensity and zero-to-moderate flow skewness (up to 
30 degrees total skewness in both pitch and yaw). Validation of the technique 
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for low-intensity turbulent flows was accomplished by analyzing data obtained 
in fully-developed pipe flow under simulated skewed flow conditions [47]. For 
moderate- intensity turbulent flows, validation was accomplished by obtaining 
data under simulated skewed flow conditions in a free jet which issued from 
fully-developed pipe flow [45]. The working forms of the mean and turbulence 
response equations developed by Al-Beirutty are presented in Appendix C with- 
out derivation. For details of the development of these equations, the reader is 
referred to Refs. 45 and 46. 

For Method A , appropriate hot-wire response equations are developed fol- 
lowing a methodology similar to that used by Al-Beirutty. The working forms 
of the response equations will be presented in Appendix C, with the details of 
their derivation included in Appendices D and E. Also, following Al-Beirutty, 
the Method A technique was partially validated in the present study by means 
of data obtained in fully- developed pipe flow under simulated skewed flow con- 
ditions. The verification procedure is described in Appendix F. 

3.5.3 Boundary Layer Parameters 

The following are definitions of some axisymmetric incompressible boundary 
layer parameters which are useful in characterizing the inlet flow condition: 

Boundary Layer Thickness 

6 = y @ UjUci = 0.995 (3.4) 


Displacement Thickness 


6x 



h^ dr 


(3.5) 


Momentum Thickness 







(3.6) 


Energy Thickness 


*3 



V* V r W 

W, )( * )dr 


(3.7) 


B = 


R U c i 


Blockage Factor 


(3.8) 
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First Shape Factor 


(3.9) 

Second Shape Factor 

s 3 

17 0 

(3.10) 


It should be pointed out here that equations (3.5), (3.6) and (3.7) are approx- 
imate, but that equation (3.8) is exact when the displacement thickness Si is 
calculated by equation (3.5). 



Fig. 3.1. Wind tunnel schematic. 



a) Side view. 

Fig. 3.2. Transition duct test section schematic. (Continued . 
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b) Cross section of transition duct, Station 1. 

Fig. 3.2. Transition duct test section schematic. (Continued . 
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c) Cross section of transition duct. Station 2. 

Fig. 3.2. Transition duct test section schematic. (Continued . 
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d) Cross section of transition duct, Station 3. 

Fig. 3.2. Transition duct test section schematic. (Continued . 
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e) Cross section of transition duct, Station 4. 

Fig. 3.2. Transition duct test section schematic. (Continued . 
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0 Cross section of transition duct, Stations 5 and 6. 
Fig. 3.2. Transition duct test section schematic. 
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Fig. 33. Transition duct data stations. 



CHAPTER 4 


INLET CONDITIONS 
(Developing Pipe Flow) 


4.1 Introduction 

The flow condition at Station 1 corresponds to partially developed turbu- 
lent pipe flow. This seemingly simple flow case has been the subject of numerous 
experimental and theoretical studies. Probably the most outstanding feature of 
the work which has been done to date is the general lack of agreement among 
the results of the various experimental investigations. Klein [48], after review- 
ing more than a dozen turbulent developing pipe flow experiments, attributed 
the disparities to the extreme sensitivity of upstream flow conditions on flow 
development. Contraction ratio, boundary layer tripping devices and starting 
conditions (smooth contraction vs. annular bleed) all influence local flow devel- 
opment downstream of the pipe inlet. In addition, the use of a boundary layer 
trip makes specification of the virtual origin of the boundary layer difficult. 

Of the experimental studies, those due to Barbin &: Jones [49], Richman 
& Azad [50] and Reichert & Azad [51] are the most complete. The facility 
used by Barbin Sz Jones incorporated a 4:1 circular contraction with an annular 
bleed. A 2.54 cm wide strip of sand particles placed 5.1 cm downstream from the 
leading edge of the pipe served as a boundary layer trip. The coordinate origin 
was coincident with the leading edge of the pipe. Mean flow and turbulence 
data were accumulated over a development length of 40 diameters for a bulk 
Reynolds number of 388,000. In contrast, the facility used by Richman &; Azad 
utilized a smooth 89:1 circular contraction connected directly to the pipe. The 
boundary layer was tripped with a 9 cm wide strip of # 16 sandpaper located 
at the entrance to the pipe. The coordinate origin was chosen to coincide with 
the downstream edge of the boundary layer trip, although the authors estimate 
the virtual origin of the boundary layer to be 3 cm upstream of this location. 
Data for this study were collected over a development length of 70 diameters 
for bulk Reynolds numbers of 100,000, 200,000 and 300,000. Reichert Sz Azad 
used the same facility as Richman & Azad but report that a 5.1 cm wide strip of 
unspecified sandpaper was used as a trip. Data for this study were collected over 
a development length of 70 diameters for seven bulk Reynolds numbers between 
112,000 and 306,000. 

The experimental data of Barbin & Jones and Richman Sz Azad were recently 
used by Martinuzzi Sz Pollard [52,53] for comparative purposes in a comprehen- 
sive evaluation of 11 turbulence models: 4 algebraic models, 2 k-e models and 
5 Reynolds stress models. All models were implemented in the same computer 
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code using equivalent boundary conditions. Although each model had its own 
merits and drawbacks, the low Reynolds number form of the k-e model overall 
performed best. 

4.2 Preliminary Results 

Preliminary measurements were made at Station 1. The purpose of these 
preliminary measurements was: 1) to establish the range of operating Reynolds 
number attainable and 2) to determine the lowest operating Reynolds number 
where fully turbulent flow exists at the first data station. The results of these 
measurements follow. 

The wind tunnel was initially configured with three inlet pipes ( L/D = 
9) and the boundary layer was allowed to develop naturally within the pipe, 
t.e. no boundary layer trip was present upstream. Pitot tube surveys of the 
boundary layer at Station 1 were obtained for two arbitrary Reynolds numbers, 
Reb = 234,000 and 403,000, the latter Reynolds number being very near the 
upper operating limit of the wind tunnel. Assuming constant static pressure 
across the data plane, mean velocity profiles were computed by equation (3.1). 
These results were plotted in law-of-the-wall coordinates using friction velocities 
deduced from measurements with the four different diameter Preston tubes shown 
in Fig. 3.2.b. Agreement with the theoretical law-of-the-wall profile was found to 
be poor, indicating that the boundary layer was not yet in a fully turbulent state. 
To promote transition to turbulence, a 2.54 cm wide strip of # 36 grit sandpaper 
was applied around the periphery of the first inlet pipe, 1.27 cm downstream 
from the joint, as shown in Fig. 3.1. With the boundary layer trip in place the 
pitot tube surveys were repeated at nominally the same Reynolds numbers. A 
comparison of mean velocity profiles with and without the boundary layer trip 
at the two Reynolds numbers is shown in laboratory coordinates in Fig. 4.1. a. 
and in law-of-the-wall coordinates in Fig. 4.1.b. These results indicate that the 
sandpaper trip was effective in producing a fully-turbulent boundary layer at the 
inlet station. 

4.3 Range of Operating Conditions 

With the boundary layer trip installed, the range of operating Reynolds 
number was found to be 0 < Re c t < 479,000 (0 < Reb < 442,000). The 
presence of the boundary layer trip, however, does not guarantee a turbulent 
boundary layer over the entire operating Reynolds number range. Indeed, the 
Preston tube calibration is valid only when law-of-the-wall behavior is present. 
Since the Preston tubes had to be removed before taking data at downstream 
stations, and the operating Reynolds numbers for data acquisition had as yet not 
been established, a correlation between friction velocity and centerline Reynolds 
number was obtained. This allows the inlet skin friction condition to be known 
for any operating Reynolds number, provided that the flow is fully turbulent. 
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This correlation for the four Preston tubes is shown in Fig. 4.2. Agreement is 
generally good for the different diameter tubes, with the exception of results 
referred to the largest tube at centerline Reynolds numbers above 425,000. The 
reason for this deviation is that the boundary layer thins as the Reynolds number 
is increased, allowing the largest tube to extend beyond the law-of-the-wall region 
where the Preston tube calibration is valid. Another small deviation, which was 
repeatable, occurs in^the neighborhood of JZe c /=220,000. The cause for this is 
unknown, but may be associated with an unsteadiness in the wind tunnel at that 
operating condition. Also indicated in Fig. 4.2 is the approximate region where 
fully turbulent flow begins. The determination of this region will be discussed 
shortly. Using the points outside the transition region, an empirical correlation 
was determined that describes friction velocity behavior at Station 1: 

= 0.255 - 0.0712 log 10 (i?e c /) + 0.005771og? 0 (Re c/ ) (4.1) 

U c i 

which is applicable in the range 97,000 < Re c i < 460,000. 

A simple indication of fully turbulent flow is whether or not law- of-the-wall 
behavior is observed within the boundary layer. Fig. 4.1 shows that law-of-the- 
wall behavior exists at Re c i = 255,000 (JJej, = 234,000). To estimate the lower 
limit of fully turbulent flow, Pitot profiles were measured at two relatively low 
Reynolds numbers, Re cl = 55,000 and 97,000 ( Re b = 52,000 and 88,000). These 
results are plotted in laboratory coordinates in Fig. 4.3.a and in law-of-the-wall 
coordinates in Fig. 4.3.b. The mean profile at the lower Reynolds number is 
considerably thinner than the higher Reynolds number profile and appears to 
have a laminar-like shape. In addition, this profile clearly does not follow law- 
of-the-wall behavior. An examination of the boundary layer shape factor, H 12 , 
indicates, however, that the flow is not fully laminar either. For the present 
profile a shape factor of Hu = 2.0 was measured, whereas typical values for 
the fully laminar and fully turbulent profiles are 2.6 and 1.4, respectively [48]. 
Conversely, the higher Reynolds number profile corresponds to a thick boundary 
layer that agrees well with the law-of-the-wall and has a measured shape factor 
of 1.4. 

Based on the above results, it was decided that data for the turbulent flow 
case would be taken at bulk Reynolds numbers of 88,000 and 390,000, for which 
the inlet flow at data Station 1 (refer to Fig. 3. 2. a) should be fully turbulent. 
The upper value represents the maximum operating speed of the wind tunnel 
throttled back slightly to allow adjustment for variations in ambient conditions. 

4.4 Results and Discussion 

For the present study, data were collected at a development length x p / D 
— 9, where x p = 0 corresponds to the pipe inlet. It is likely that the effective 
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length is slightly larger, inasmuch as a boundary layer is already developing 
in the contraction and the sandpaper trip tends to thicken the boundary layer 
artificially. The transverse flow angle in the azimuthal direction was measured 
using a pressure-nulling technique with a two-tube Conrad probe. The probe was 
first nulled at the pipe centerline and then the flow angle along the yi traverse (see 
Fig. 3.2.b) was measured relative to the centerline null value. For both operating 
Reynolds numbers, the computed transverse flow velocity (U$) was found to be 
less tlian 0.25% of the local axial velocity component across the entire traverse. 

Mean velocity profiles plotted in laboratory coordinates and in law-of-the- 
wall coordinates along four equally-spaced radial traverses at Station 1 are shown 
for Ref, = 88,000 and 390,000 in Figs. 4.4 and 4.5, respectively. The excellent 
peripheral symmetry of the flow and agreement with law-of-the-wall behavior 
are readily apparent. The non-dimensional Preston tube diameters (d + ) used 
to deduce skin friction are indicated on the law-of-the-wall plot. The largest 
tube for the higher Reynolds number case extends slightly into the wake region 
and, as a result, skin friction values deduced from this tube were slightly higher 
than those for the three smallest tubes. Friction velocities deduced from the three 
smallest tubes deviated by less than 0.3% from their mean value. Boundary layer 
thickness, integral parameters, skin friction and centerline turbulence intensity 
at Station 1 are summarized in Table 4.1. The integral parameters shown are 
averages of values computed from the four individual traverses. 


Table 4-1- Flow condition at Station 1 



Re^ = 88, 000 

Re b =390,000 

6/R 

0.3141 

0.2855 

Si JR 

0.0448 

0.0383 

62 /R 

0.0312 

0.0281 

Sz/R 

0.0545 

0.0497 

B 

0.0896 

0.0765 

H12 

1.438 

1.364 

H 32 

1.748 

1.771 

T w / D 

1.00 

0.96 

( u'/U) c i 

0.008 

0.003 


R = 10.214 cm 


As expected, the boundary layer is somewhat thicker for the lower Reynolds 
number case. The Hu shape factors agree well with those reported by Klein 
for fully turbulent, developing pipe flow. The wall shear stress in Table 4.1. has 
been normalized by the fully- developed pipe flow shear stress t Wi pd which was 
determined from a relation given by Eppich [40]. Data presented by Reichert 
& Azad show that when a smooth contraction is used as the starting condition, 
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the skin friction at the entrance to the pipe is below the fully-developed value 
and that the skin friction approaches the fully- developed value in an oscillatory 
manner. In contrast, Barbin k Jones data show that when a sharp edge is used 
as a starting condition, the skin friction at the beginning of the pipe is above the 
fully-developed value but asymptotically reaches the fully-developed value within 
the first 15 diameters. The present results show that the wall shear stress at x p f D 
= 9 is very near the fully-developed value. At the lower Reynolds number, the 
flow was much more susceptible to ambient disturbances and this is reflected in 
the larger data scatter and higher freestream turbulence levels. 

Reynolds normal and shear stress distributions at Station 1 for Re\> = 88,000 
and 390,000 are shown in Fig. 4.6. Experimental results obtained by Barbin k 
Jones and Richman k Azad at a comparable location {x p fD = 10) and operating 
Reynolds numbers are also shown for purposes of comparison. The turbulence 
levels shown in these plots near the pipe centerline should not be considered 
representative of the core flow turbulence intensity, inasmuch as the output from 
the probe used to measure the distributions exhibited a small sinusoidal trace on 
the oscilloscope when positioned in the core flow. This noise was attributed to 
vortex shedding from the probe prongs when the probe was positioned normal 
to the axial mean flow. The turbulence intensity values shown in Table 4.1. are 
correct and were measured independently using a forward-facing probe fixed at 
the duct centerline. 

The Reynolds stress profiles in Fig. 4.6 are seen to be nearly Reynolds num- 
ber independent over the range considered. Based on the turbulence results, it 
appears that the boundary layer in the present study is thinner than that of the 
other investigators. In the near-wall region, the normal stresses are seen to be 
significantly higher than those due to Barbin k Jones. The shear stress in the 
near-wall region agrees well with the results of Richman k Azad and with the 
fully- developed distribution. 

The more rapid boundary layer growth associated with Richman and Azad’s 
experiments can probably be attributed to the rougher and wider sandpaper 
tripping device used in their study, which extended over a streamwise width of 
0.9 pipe diameters in comparison to the width used in the present study (0.12 
pipe diameters). In Barbin k Jones’ experiments, a sand grain trip of the same 
relative width as that used in the present study was employed. By-pass bleed 
through an annular gap between the pipe and a 4:1 contraction ratio nozzle was 
used to promote spanwise flow uniformity across the pipe inlet. In reference to 
Fig. 1 of Ref. 49, the pipe in Barbin k Jones’ experiments extended upstream 
into the nozzle, so that the pipe inlet could have been in a region where the flow 
■was still converging, and the flow may have separated at the wedge-shaped lip of 
the pipe. This event, if it occurred, could have promoted more rapid boundary 
layer growth in comparison to that observed in the present study. If local flow 
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separation and reattachment did occur near the pipe inlet in Barbin & Jones 
experiments, then one would also anticipate that flow in the near-entrance would 
not be in local equilibrium. This conjecture is supported to some extent by the 
behavior of the velocity profiles measured by Barbin & Jones at i p /D=1.5, 4.5, 
7.5, all of which show significant departures from the law-of-the-wall across the 
width of the boundary layer (refer to Fig. 6 of Ref. 49). 

The differences between corresponding distributions in Figs. 4.6.a and 4.6.b 
could be reconciled further if a code were available for predicting developing 
turbulent pipe flow with high accuracy starting from a low- turbulence level, uni- 
form inlet flow condition. Unfortunately, predictions of a given variable in the 
entrance region based on various k-e type transport equation models show dif- 
ferences which are the same order of magnitude as those which exist among the 
various data sets for nominally the same operating conditions [52,53]. Further 
work will be required in order to: (1) provide comprehensive data which charac- 
terize disturbance-free, turbulent boundary layer growth in a circular pipe, and 
(2) develop a code which can predict this behavior, recognizing that the data 
may have to be corrected for virtual origin (streamwise displacement) effects. It 
should be noted here that the present data set at Station 1 appears to be rela- 
tively free of upstream disturbances, inasmuch as the core flow at this station is 
uniform and at a relatively low turbulence level (as indicated in Table 4.1) and 
velocity profiles measured along four radial traverses 90° apart are symmetric 
and in excellent agreement with the law-of-the-wall (refer to Figs. 4.4 and 4.5). 

The turbulence kinetic energy k and its rate of production are shown for 
the two operating Reynolds numbers in Figs. 4.7 and 4.8, respectively. In Fig. 
4.7, the conventional wall function value based on C f <1 =0.09 is shown. The lower 
Reynolds number data agree well with the limiting value, but the higher Reynolds 
number data is slightly high. The higher kinetic energy in the near-wall region 
at the higher Reynolds number is supported by the kinetic energy production 
distributions shown in Fig. 4.8. Also indicated are results based on the data of 
Barbin and Jones, which lie considerably below the present results. The results 
shown in Figs. 4.7 and 4.8 will be used later as reference distributions to as- 
sess the distorting effect of the transition duct geometry on the turbulence field. 
Non-dimensional turbulence structure parameters, which will also be useful for 
latter comparisons, are the shear stress correlation R uv = uv/u'v' and the shear- 
energy ratio parameter al = uv /2k. These parameters are sometimes prescribed 
as constants in Reynolds-averaged turbulence models. Distributions measured 
at Station 1 are shown in Figs. 4.9 and 4.10, respectively. The shear stress corre- 
lation may be interpreted as the ability of the normal stresses to generate shear 
stress. In a two-dimensional boundary layer, this parameter is nearly constant 
within the layer with a typical value of about 0.45 [54]. The shear-energy ratio 
parameter has also been observed to be nearly constant in a two-dimensional 
boundary layer with a value of about 0.15. With the exception of the scatter 
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near the boundary layer edge (at r/R « 0.72) for the lower Reynolds number, 
the present results generally agree well with these values. 
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b) Law-of-the-wall profiles. 

Fig. 4.1. Tripped and non-tripped boundary layer 
profiles at Station 1, Reb = 234,000 and 403,000. 
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Fig. 4.2. Friction velocity vs. centerline Reynolds number at Station 1. 
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a) Mean velocity profiles. 



b) Law-of-the-wall profiles. 

Fig. 4.3. Tripped boundary layer profiles at Station 1, Reb = 52,000 and 88,000. 
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a) Mean velocity profiles. 



b) Law-of-the-wall velocity profiles. 

Fig. 4.4. Boundary layer profiles at Station 1, Reb = 88, 000. 
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a) Mean velocity profiles. 
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b) Law-of- the- wall velocity profiles. 

Fig. 4.5. Boundary layer profiles at Station 1, Re*, = 390, 000. 
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b) Shear stress. 

Fig. 4.6. Reynolds stress distributions at Station 1. 
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Fig. 4.7. Turbulent kinetic energy distributions at Station 1. 
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Fig. 4.8. Production of t.k.e. distributions at Station 1. 
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Fig. 4.9. Shear stress correlation distributions at Station 1. 
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CHAPTER 5 


RESULTS AND DISCUSSION 


5.1 Introduction 

The sequence in which data were collected in the transition duct is as follows. 
The first transition duct was installed and the peripheral wall static pressure 
distributions at each data station were obtained. Data were then accumulated at 
Station 5 to check flow symmetry and determine the extent of Reynolds number 
dependence of the flow. The transition duct extension was then installed and 
data were accumulated at Station 6. Finally, the second transition duct was 
installed and measurements were made at the intermediate Stations 3 and 4. No 
data were taken at Station 2. Total pressure and mean velocity were deduced 
from hot-wire measurements made at Stations 3,4,5 and 6. Turbulence quantities, 
however, were only measured at Stations 5 and 6 since the anticipated vortex 
structure was not well defined until Station 5. 

At all stations, symmetry was assumed about the x-z plane and, as such, 
data were accumulated only in Quadrants 1 and/or 2 shown in Figs. 3.2.d-3.2.f. 
At Stations 3 and 4, due to the arrangement of the probe access holes (see Figs. 
3.2.d and 3.2.e), data were taken only in Quadrant 2. For presentation purposes, 
the data were imaged about the y-axis to show an entire duct half. All subsequent 
plots where data have been imaged are so indicated. The spacing between data 
points along individual traverses varied from 1.27 cm in the core region to 0.127 
cm near the wall surface. At Stations 5 and 6, rather than assuming symmetry 
about the midplane z=0, measurements were made in both Quadrants 1 and 
2. Data were taken at selected points on a 63 x 19 rectangular grid (0.254 cm 
spacing) in each quadrant. The spacing between data points varied from 2.54 
cm in the core region to 0.254 cm in regions of large gradients. Approximately 
400 data points were taken in the duct half. It should be emphasized here that 
all subsequent contour results for Stations 5 and 6 are based on data obtained in 
both quadrants, so that the level of symmetry about the x-y plane of the duct is 
a direct indication of the quality of the flow and of the measurement techniques 
employed in this study. For plotting purposes, the measured results at all stations 
were interpolated onto an evenly spaced mesh by means of a monotonic derivative 
spline interpolant [55]. This method of interpolation performs no smoothing and 
guarantees no overshoot of the data. 

5.2 Static Pressure Distribution 

Wall static pressure distributions were measured along the periphery (s- 
coordinate, refer to Figs. 3.2.d through 3.2.f) of the lower half of the duct at 
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Stations 3,4,5, and 6. Peripheral wall static pressure coefficient distributions are 
shown for Re b = 88,000 and 390,000 in Fig. 5.1. The normalizing dimension s ref 
is 1 /4 of the duct circumference at a given data station. Solid symbols represent 
static pressure measured along the duct centerline. The excellent spanwise regu- 
larity at all stations supports the assumption of symmetry about the x-z plane. 
The primary difference between results referred to the two Reynolds numbers 
is that the net pressure coefficient drop along the duct is larger for the lower 
Reynolds number case and is the result of increased viscous losses associated 
with the more rapid boundary layer growth. 

Local static pressure is a function of cross-sectional area, wall curvature 
and viscous forces. Station 3 is located at an axial position where streamwise 
diffusion of the flow is occurring. In reference to Fig. 5.1, this results in a net 
rise in static pressure above the inlet value (C p =0) which is reflected in the 
centerline pressure. Concave curvature along the upper wall (s/s re /=0) induces 
a positive pressure gradient ( dP /dr > 0) which results in the observed pressure 
peak. Conversely, convex wall curvature along the sidewalls (s/s re /=± 1) induces 
a negative pressure gradient resulting in the pressure minima. A net decrease 
in static pressure occurs between Stations 3 and 4 which is the result of a slight 
decrease in cross-sectional area and viscous losses. Also, the radius of curvature 
of the walls changes sign between Stations 3 and 4 which causes maximum and 
minimum pressures to occur along the side and upper walls, respectively. At 
Station 5 the area has returned to the inlet value so that a net decrease in pressure 
relative to the inlet occurs which is due only to viscous effects. Although there 
is no curvature of the walls at this station, upstream curvature effects (the wall 
curvature peaks between Station 4 and 5) are still strongly present. At Station 6, 
there is no area change or curvature effects, and the static pressure is nominally 
constant across the entire cross section. 

5.3 Effect of Varying Reynolds Number 

Results presented thus far have shown that, aside from a thicker boundary 
layer at the lower operating Reynolds number, no appreciable Reynolds num- 
ber dependence exists over the limited range considered. In order to determine 
the influence of Reynolds number on downstream flow development, mean flow 
measurements and some limited turbulence data were accumulated at Station 
5 for operating bulk Reynolds numbers of 88,000 and 390,000. Total pressure 
contours, axial mean velocity contours and transverse velocity vectors for the 
two operating conditions axe shown in Figs. 5.2, 5.3 and 5.4, respectively. These 
results are shown here only for comparison; the physical significance of the flow- 
field behavior will be discussed in the next section. In Fig. 5.4, the reference 
velocity vector represents the largest measured vector in the plane. The mean 
flowfields at the two operating conditions are observed to be very similar. One 
notable difference is that the vortex pair for the lower Reynolds number is more 
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circular and centered further away from the wall. Overall, though, the differences 
observed were not significant enough to justify repeating all measurements for 
both operating conditions. It was decided, therefore, to restrict the bulk of the 
measurements at the remaining stations to one operating condition. Since, for 
most practical applications, the operating condition for a circular-to-rectangular 
transition duct would be closer to the higher operating condition, most of the 
remaining measurements are for the Rt b = 390,000 case.- 


5.4 Mean Flow Results 
5.4.1 Mean flow contours 

Mean flow variables were measured in the transverse plane at Stations 3,4,5 
and 6. Reference coordinates applicable to theses data stations are shown in 
Figs. 3.2.d through 3.2.f, respectively. Total pressure contours (Pitot probe) 
and axial velocity contours (hot-wire probes) at Stations 3,4,5 and 6 are shown 
in Figs. 5.5 and 5.6, respectively. At Station 3, the total pressure and axial 
velocity contours generally follow the cross- sectional shape of the duct. At 
Station 4, a distortion of the contours is seen to develop in the vicinity of the 
sidewall. The data at Stations 5 and 6 show that the distortion grows in the axial 
direction. The development of the distortion is qualitatively very similar to the 
results presented by Taylor ef o2. [39] for turbulent flow through a non-diffusing 
S-shaped duct. Transverse velocity vectors measured by means of hot-wires are 
shown in Fig. 5.7. The magnitude of the reference velocity vector indicated in the 
plots represents the maximum vector magnitude observed at that station. The 
distortion of the primary flow is due to a secondary flow pattern which develops 
into a discrete vortex pair along the duct sidewalls. This secondary flow arises 
as a result of lateral skewing of the near-wall flow in the vicinity of the sidewall 
induced by transverse pressure gradients (refer to Fig. 5.1) associated with wall 
curvature. At Station 5, the vortices are oblong in shape and their centers are 
positioned relatively near the duct sidewall. At Station 6, the vortices have 
grown in lateral extent, are more circular and are centered further away from the 
sidewall. Between Stations 2 and 5, lateral divergence of the upper (lower) wall 
boundary layer and lateral convergence of the sidewall boundary layer causes 
thinning and thickening of the boundary layer, respectively. Between Stations 
5 and 6, the boundary layer everywhere thickens by natural growth and, in the 
vicinity of the sidewall, by lateral convergence. 

The static pressure distribution at Station 5 was measured directly with a 
static pressure probe and calculated from the total pressure and velocity distri- 
butions: 

P = Pt- ^PambUZ/(R a ,rT amb ) (5.1) 

where the variables Eire defined as in equation (3.1). Measurement of the static 
pressure at Station 5 was possible since the total flow angle is less than 10°. 



50 


The measured and calculated distributions are shown in Fig. 5.8. The measured 
static pressure in the near wall region agrees well with the wall tap measurements 
shown in Fig. 5.1.b. The calculated pressure field qualitatively agrees with the 
measurements, but doesn’t show the saddle- shaped distributions as clearly. 

5.4.2 Boundary layer divergence 

The mean velocity distributions were used to estimate the strength of the 
divergence along the midplanes y=0 and z=0 (negative divergence is the same 
as convergence). On the midplane y=0 (upper wall), the divergence parame- 
ter is defined as (dV / dy) / (dU / dy i) evaluated midway through the boundary 
layer. On the midplane z = 0 (side wall), the divergence parameter is defined as 
( dW/dz)/(dU/dy2 ) evaluated at the first data point from the wall (« 0.5 cm). 
An estimate of the axial distribution of these divergence parameters is shown in 
Fig. 5.9. Along the duct upper (lower) surface, the divergence is, for the most 
part,. confined to the actual transition section, peaking at around 20%. To fur- 
ther illustrate the degree of divergence on the upper and lower surfaces of the 
duct, Fig. 5.10 reproduces surface oil flow results by Reichert et al. [4] obtained 
in an identical transition duct at an operating condition of Re c i t i = 1.57 x 10 6 . 
Along the duct sidewalls, the flow begins to converge in a manner similar to the 
divergence on the upper wall, but then a large jump is observed between x/R = 
2.8 and 4.0 (Stations 4 and 5). This jump is due to a reduction in the primary 
strain rate dJJ/dy2 which is a direct result of the vortex pair. The results at 
x/R = 8.0 shows that the boundary layer continues to converge well into the 
transition duct extension. Based on the strength of the divergence observed, 
it is anticipated that divergence effects will be reflected in the local turbulence 
structure. 

5.4.3 Mean flow profiles 

Total pressure profiles, at Stations 1,3, 4, 5 and 6, measured along the y 2 
(sidewall) and y% (upper wall) traverses (see Fig. 3.2) are shown in Figs. 5. 11. a 
and 5.11.b, respectively. Similarly, axial velocity profiles are shown in Figs. 5. 12. a 
and 5.12.b. These results show that, between Stations 1 and 5, the boundary 
layer substantially thickens along the y 2 traverse and thins along the y$ traverse. 
This is due to the secondary flow transferring boundary layer fluid along the duct 
periphery from the vicinity of the y 3 traverse to the vicinity of the y 2 traverse. 
Beyond Station 5, where the cross-sectional shape is constant, the boundary 
layer along both traverses thickens. Along the y 3 traverse, the thickening is due 
to natural boundary layer growth, but along the y 2 traverse, the thickening is 
due also to the common outward flow associated with the vortex pair. 

At Stations 5 and 6, the total pressure and velocity profiles along the y 2 
traverse exhibit a double inflection behavior. The total pressure contours (Figs. 
5.5.c and 5.5.d) and the axial velocity contours (Figs. 5.6.c and 5.6.d) indicate 
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that along traverses adjacent and parallel to the y 2 traverse, a double peaking 
behavior is observed. This is a result of the vortex pair convecting a ridge of 
high momentum fluid from the potential core flow to the region along the duct 
sidewalls, re-energizing the boundary layer and very probably preventing flow 
separation. A break in the ridge occurs at the midplane (y 2 traverse) due to a 
transfer of low momentum fluid from the boundary layer toward the centerline, 
creating a flat spot in the velocity field. This flat region is seen to be much larger 
at Station 6 (Fig. 5.6. d) than at Station 5 (Fig. 5.6.c). 

5.4.4 Streamwise vorticity 

The presence of the vortex pairs in the exit plane of the transition duct 
is undesirable, inasmuch as they cause significant regions of total pressure loss. 
Streamwise vorticity in non- circular ducts can be generated by two different 
mechanisms. The first is vorticity generation by the lateral deflection (by pressure 
gradients) of a shear layer with spanwise vorticity. This mechanism is often 
referred to as skew-induced streamwise vorticity. An example of this type is the 
horseshoe vortex generated by a blunt obstruction in a 2-D boundary layer. The 
second mechanism is streamwise vorticity created by the Reynolds stresses and 
is referred to as stress-induced vorticity. Streamwise vorticity generation by the 
Reynolds stresses occurs in turbulent flows through straight non-circular ducts. 
Generally, stress-induced vorticity is much weaker than skew-induced vorticity. 
These mechanisms are represented in the steady axial mean vorticity equation. 
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The LHS of equation (5.2) represents the increase in streamwise vorticity by 
convection. The first term on the RHS represents production of streamwise 
vorticity by vortex line stretching (streamwise acceleration causes amplification 
of vorticity). The second and third terms on the RHS represent the increase in 
vorticity due to lateral skewing (by transverse pressure gradients) of vorticity 
in the transverse directions. These are the terms associated with skew-induced 
vorticity. The fourth term represents the production of streamwise vorticity by 
the primary shear stresses and is often neglected since it contains a streamwise 
gradient. The fifth and sixth terms on the RHS represent the production of 
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streamwise vorticity by inhomogeneity of the transverse normal stress anisotropy 
and by the secondary shear stress, respectively. These terms are responsible for 
stressed-induced vorticity. The last term on the RHS represents the diffusion of 
streamwise vorticity by viscous forces. 

The significance of equation (5.2) can now be discussed relative to the present 
tr ansi tion duct configuration. Since the sign of axial vorticity changes across 
planes of symmetry, to avoid confusion, the following discussion will be restricted 
to Quadrant 1 (refer to Fig. 3.2) of the transition duct. The axial vorticity com- 
ponent at Stations 3,4,5 and 6 was calculated by interpolating the transverse 
velocity components onto a uniform grid (0.508 x 0.508 cm) and then evaluating 
the derivatives by central difference approximations. Contours of axial vorticity 
at Stations 3,4,5 and 6 are shown in Fig. 5.13. Negative vorticity is represented 
by dashed contour lines. The filled circles in these plots mark the approximate 
location where the peak vorticity occurs. At Station 1, the vorticity field is com- 
prised only of transverse vorticity components fl v and the axial component 
Q x being zero. Beginning at Station 2, the wall curvature induces the transverse 
pressure gradients which were discussed in Section 5.3 and illustrated in Fig. 5.1. 
In the first half of transition, a global deceleration of the flow occurs due to the 
area expansion. Fig. 5. 13. a (Station 3) shows that in Quadrant 1 the transverse 
pressure gradient creates primarily negative vorticity. The no-slip condition re- 
quires there to be a thin layer of positive vorticity in the very near-wall region 
which was not resolved in the present measurements. The generation of nega- 
tive vorticity will occur in any straight CR transition duct without swirl since 
the pressure gradient is primarily a function of wall curvature. In the second 
half of transition, the wall curvature changes sign which causes a reversal of the 
transverse pressure gradient (see Fig. 5.1). Also, the contracting area causes a 
global acceleration of the flow. It might be expected that the reversal of the 
pressure gradient would effectively cancel the vorticity generated in the first half 
of transition. This is not the case, however, as Figs. 5.13.b and 5.13.C (Stations 4 
and 5), show that the negative vorticity migrates towards the midplane 2=0 and 
intensifies. This strengthening of the vorticity is caused by streamwise accelera- 
tion (vortex stretching) and concave surface curvature along the duct sidewalls. 
Although the vortices are not generated by centrifugal instabilities, concave cur- 
vature will accentuate the vorticity in the same manner that Taylor- Gortler type 
vortices intensify. In addition, the reversed pressure gradient assists the inward 
(negative y direction) flow between the developing vortex pair. In the transition 
duct extension (Fig. 5.13.d, Station 6), the vorticity is diffused by turbulent ac- 
tion. Between Stations 5 and 6, the magnitude of the peak vorticity drops from 
Q X R/U b = 2.4 to 0.09. 

Miau et al. [9] evaluated all the terms in equation (5.2) at the exit plane 
of two CR transition ducts of constant cross- sectional area. Both ducts had 
aspect ratios (Ai?) of two, but the transition lengths differed: L/D= 0.54 and 
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L/D=z 1.08. The results of their analysis showed that the generation of axial 
vorticity at the exit plane was primarily skew- induced (terms 3 and 4), rather 
than turbulence induced (terms 5-7). A notable difference between Miau s results 
and the present study is that the axial vorticity in the corner region at the exit 
plane is of opposite sign. Whereas negative axial vorticity is observed in Quadrant 
1 (upper quadrant, Fig. 5.13.c) of the present study, positive axial vorticity is 
observed in the equivalent quadrant of their ducts. Patrick and McCormick [2,3] 
calculated the axial vorticity from measurements in the exit plane of two CR 
ducts. One duct had an aspect ratio of three and a length-to-diameter ratio 
of one, and the other duct had an aspect ratio of six and a length-to-diameter 
ratio of three. The results of this investigation showed small regions of positive 
vorticity for the AR = 3 duct and large regions of negative vorticity for the AR = 6 
duct. A common feature between the ducts which exhibited positive vorticity 
is that the cross-sectional area through the duct was constant. In contrast, the 
ducts which exhibit negative vorticity in the exit plane had an area expansion 
followed by a contraction through the transition section. For all ducts, the inlet 
area equalled the exit area. Another distinguishing feature is the transition length 
was considerably shorter for the positive vorticity ducts than for the negative 
vorticity ducts: L / D—0.5A and 1.08 vs. 1.5 and 3.0, respectively. These results 
suggest that the condition of the flowfield in the exit plane is very sensitive to the 
geometry of the transition. More specifically, the path that the vortex core takes 
through the duct will determine if the initial vorticity is amplified or attenuated. 

5.5 Turbulence Results 

The complete Reynolds stress tensor was measured at Stations 5 and 6. A 
distinguishing feature of the flowfield at these stations is the turbulence structure 
in the vicinity of the vortex pair. In Section 2.3, several detailed experimental in- 
vestigations of embedded vortex flows were mentioned. In particular, the studies 
by Mehta and Bradshaw [30] and Pauley and Eaton [31,32] are relevant, inasmuch 
as they present detailed mean flow and turbulence results for embedded vortex 
pairs with the common flow away from the surface, as occurs in the transition 
duct. The present measurements afford the opportunity to make comparisons 
with these data sets to determine if modelling conclusions based on their results 
are applicable to the transition duct flow. For comparative purposes, the origin 
of the vortices in the transition duct is taken as the location where the geometry 
deformation begins (Station 2). Based on this, the development length for the 
vortices is approximately 30 and 70 cm at Stations 5 and 6, respectively. The 
first measurement station reported by Pauley and Eaton was 97 cm downstream 
from the generators and Mehta and Bradshaw present most of their results at 
a development length of 135 cm. In addition to development length, other dif- 
ferences between the studies should be noted. Whereas the other investigators 
studied vortices embedded in a two-dimensional boundary layer, the transition 
duct boundary layer in the region of the vortices is three-dimensional. Also, the 
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velocity deficit in the vicinity of the vortex cores is much larger in the transition 
duct than in the other studies. Mehta and Bradshaw placed their generators in 
the settling chamber so that the velocity deficit was a small percentage of the po- 
tential core flow velocity by the time the vortices entered the test section. Pauley 
and Eaton, on the other hand, generated their vortices at the beginning of the 
test section, so that a significant deficit was initially present, but had diminished 
somewhat at the downstream data station where the turbulence measurements 
were made (97 cm). Velocity deficit in the vortex core region is summarized in 
Table 5.1. 


Table 5.1. Velocity deficit in the vicinity 
of the vortex cores. 



x (cm) 

u/u c 

Present results 

30 

0.40 

Present results 

70 

0.50 

Pauley & Eaton [31] 

66 

0.50 

Pauley & Eaton [31] 

97 

0.75 

Pauley & Eaton [31] 

142 

0.90 

Mehta & Bradshaw [30] 

60 

0.95 

Mehta & Bradshaw [30] 

90 

0.95 

Mehta & Bradshaw [30] 

135 

0.95 


Based on these differences alone, it is expected that the present results will 
qualitatively be more similar to the results of Pauley and Eaton than to those of 
Mehta and Bradshaw. 

5.5.1 Turbulence contours 

Axial turbulence intensity u' /U contours measured at Stations 5 and 6 are 
shown in Fig. 5.14. Note that the turbulence intensity here is defined relative 
to the local axial velocity component and not the bulk velocity. Since the local 
intensity exceeds 10% in places, second-order response equations (see Appendix 
C) were used for reducing all of the hot-wire data. The filled circles in Fig. 5.14 
and subsequent figures mark the location of the vortex cores (peak vorticity) as 
an aid to interpretation. In the vicinity of the sidewall, the vortex pair extrudes a 
tongue of moderate turbulence intensity; (u 1 /U) ma x ~ 14% and 11% at Stations 
5 and 6, respectively. The increase in boundary layer thickness and distortion by 
the vortex pair between Stations 5 and 6 is clearly evident. Along the midplane z 
= 0, the turbulence intensity exhibits a double peak behavior which is a result of 
the double inflection observed in the mean velocity profiles shown in Fig. 5. 12. a. 

Contours of the six Reynolds stress components at Stations 5 and 6 are 
shown in Figs. 5.15—5.20. Negative contour levels are represented by dashed 
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lines. The level of symmetry for each stress component about the midplane 
z fR _ o is very good, even for the difficult-to-measure vw stress component 
shown in Fig. 5.20. As expected, the uw and vw stress components change sign 
between Quadrants 1 and 2. In general, Figs. 5.15—5.20 show that the Reynolds 
stress contours are more distorted at Station 6 than at Station 5, and that peak 
contour values are higher and displaced farther from the duct sidewall. The 
increase in peak magnitude of the stresses is in contrast to the decrease in peak 
axial turbulence intensity (Fig. 5.14) which resulted from an increase in the local 
axial velocity in the vicinity of the high turbulence levels. 

At both_Stations_5 and 6, the peak magnitude of the transverse normal stress 
components v 2 and to 2 are nearly equal and are roughly half the peak magnitude 
of the axial stress component u 2 . Qualitatively, the u 2 stress component differs 
significantly between all the vortex studies. The results of Mehta and Bradshaw 
at 135 cm show that as the wall is approached along the plane of symmetry, the 
u 2 component monotonically increases. The results of Pauley and Eaton at 97 
cm show a peak in the near-wall region, but also show an additional peak on 
the plane of symmetry near the outer edge of the vortices. They also show a 
pair of peaks in the vortex core region symmetrically located about the plane 
of symmetry. The present results show a large peak on the plane of symmetry 
at the edge of the vortex pair. These differences can be traced to the velocity 
deficit (or lack of) in the vortex core region. It is likely that the vortices in the 
Pauley and Eaton study exhibited similar behavior to the present results at a 
location closer to the generators since a large velocity deficit is prese nt d irectly 
behind the generators. At Station 6, the present results show that the w 2 stress 
component exhibits two peaks symmetrically located about the y 0 midplane 
which were not observed at Station 5. Pauley and Eaton s data at 97 cm suggest 
that all the normal stresses will eventually exhibit this double peak behavior. 

Anisotropy between the axial stress component and the horizontal transverse 
component (u 2 -v 2 ) is shown in Fig. 5.21. At Station 5, the axial component is 
everywhere larger than the horizontal component. In the vicinity of the vortex 
core at Station 6, however, the transverse component exceeds the axial compo- 
nent by as much as 20%. This is in contrast to a two-dimensional boundary layer 
where the axial component is greater than either transverse components. Pauley 
& Eaton and Mehta k Bradshaw both report similar behavior of the normal 
stresses in the vicinity of the vortex core. Anisotropy between the axial and ver- 
tical normal stress components (u 2 — w 2 )) is shown in Fig. 5.22. Near the vortex 
core at Station 5, the stresses axe observed to be nearly equal. In the vicinity 
of y/R = 1.4, i:z/R = 0.35, the vertical stress exceeds the axial component. At 
Station 6, the vertical stress is everywhere less than the axial component, but 
they are still nearly equal in the vortex core region. Anisotropy between the ver- 
tical and horizontal transverse normal stress components (v 2 — w 2 ) is important 
in the generation of streamwise vorticity in non-circular ducts (see Section 5.4.1). 
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Contours of this quantity are plotted in Fig. 5.23. These plots show significant 
differences betweenJStations 5 and 6. It is_well known that in two-dimensional 
boundary layers, uj > uf > u\, where uj is the axial component, is the 
transverse component tangential to the surface, and u| is the component which 
acts normal to the surface (see, e.g., Fig. 18.5 of Ref. 54). For the transition 
duct flow, the above inequality is observed to hold true in the near-wall region at 
both Stations 5 and 6. In the outer region of the upper and lower wall boundary 
layers at Station 5, however, the component tangent to the surface exceeds the 
component which acts normal to the surface (negative anisotropy). In the region 
of the vortex pair at Station 5, a small pocket exists where the anisotropy is pos- 
itive. By Station 6, the region of positive anisotropy has grown and intensified. 
The results at Station 6 are qualitatively very similar to the results reported by 
Pauley and Eaton for a vortex pair at a development length of 97 cm (see Fig. 
4.39 of Ref. 31). 

The primary shear stress uv shown in Fig. 5.18 is observed to be positive 
everywhere except for a small region in the vicinity of the vortex cores. These 
regions of negative stress are a result of the positive mean rate-of-strain (dU jdy > 
0) associated with the aforementioned velocity ridge (see Section 5.4.3). Along 
the midplane z=0 where the axial mean velocity profile flattens, the primary 
mean rate-of-strain in both the y and z directions is nearly zero and, as a result, 
stress levels in this region are depressed. In the z-direction, the primary shear uw 
shown in Fig. 5.19 also changes sign as the summit of the velocity ridge is crossed. 
The secondary shear stress vw is another quantity that plays an important role 
in the production of streamwise vorticity (see Section 5.4.1). In the vicinity of 
the vortex cores, this stress is of the same order of magnitude as the primary 
shear stresses. 

The normal stress data were used to calculate the turbulence kinetic energy. 
These results are shown for Stations 5 and 6 in Fig. 5.24. The production of 
kinetic energy was also calculated. Neglecting streamwise derivatives, the pro- 
duction of kinetic energy is given by: 


P = 


dU dU 



uw— 

dz 


(5.3) 


Contours of this quantity are shown in Fig. 5.25. Generally, high levels of kinetic 
energy are associated with high production rates. A notable exception occurs 
on the duct upper and lower walls at Station 5. Here, the production in the 
near-wall region is significantly less than that observed at Station 6, which seems 
disproportionate when compared to the relatively small difference in kinetic en- 
ergy between these stations. Between Stations 5 and 6, the peak production 
is essentially the same, but the peak kinetic energy increases by approximately 
35%. 
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5.5.2 Turbulence profiles 

In this section, Reynolds stress profiles along the duct semi-major and semi- 
minor axes are presented and compared with the initial distributions measured 
at Station 1. In the following profile plots, the designation of the Reynolds 
stress components are relative to wall coordinates and not the z,y,z laboratory 
coordinates; that is, the v fluctuating velocity component is always directed along 
the wall coordinate of interest, either j /2 or 1/3 (see Fig. 3.2). All of the stresses 
are normalized by the bulk velocity and the wall coordinate is normalized by the 
local boundary layer thickness, as summarized in Table 5.2. 


Table 5.2. Normalized boundary layer thickness 
(S/R) at Stations 1,5 and 6 . 



Reb — 

88,000 

Rt\> — 

390,000 

Station 

1 / 2 -axis 

1 / 3 -axis 

y 2 -axis 

y 3 -axis 

1 

0.31 

0.31 

0.29 

0.29 

5 

- 

0.21 

0.58 

0.18 

6 

- 

- 

0.65 

0.25 


R = 10.214 cm 


Reference normal Reynolds stress distributions measured at Station 1 are shown 
in Fig. 5.26. Normal stress distributions measured at Stations 5 and 6 along 
the semi-major axis are shown in Fig. 5.27 and shear stress and kinetic energy 
profiles are shown in Fig. 5.28. It is readily apparent that the distributions along 
the semi-major axis deviate considerably from the initial profiles. The boundary 
layer on the duct sidewalls is subjected first to stabilizing convex curvature and 
stabilizing lateral convergence. In the second half of transition, the developing 
sidewall vortex pair creates stabilizing lateral convergence near the wall and 
destabilizing lateral divergence in the outer region of the boundary layer. In 
addition, the flow experiences destabilizing concave curvature. The common 
flow away from the wall creates a region of velocity deficit resulting in the double 
infection behavior of the mean velocity profiles shown in Fig. 5 . 12 . a. At Station 
5 , the axial normal and shear stress components show the largest deviation with 
strong attenuation in the near-wall region. Near the wall at Station 6 , these stress 
components have increased to a certain extent, but are still below the initial 
levels. The transverse normal stress components in the near-wall region decrease 
between Stations 5 and 6 . In wall-bounded shear layers, the Reynolds stresses 
are known to scale with the local friction velocity. Preston tube measurements 
were made in the duct mid-planes at Stations 1,5 and 6 from which the local 
friction velocity was deduced. These results are summarized in Table 5.3. 
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Table 5.S. Normalized friction velocity 
{Ur/Ub x 100) at mid-planes, JJe* = 390,000. 


Station 

O 

II 

s 

O 

II 

S 

1 

4.06 

4.06 

5 

2.70 

4.89 

6 

3.23 

4.38 


U h = 29.95 m/s 


The drop and subsequent rise in the axial Reynolds normal and shear stress 
components near the wall between Stations 1 and 6 correlates with the friction 
velocity behavior shown in Table 5.3. The transverse normal stresses appear 
to lag the development of the friction velocity, although it is presumed that 
all the stress components in the near-wall region will eventually increase as the 
vortex pair moves further away from the duct sidewalls. The decrease in the 
transverse normal stresses outweighs the increase in the axial normal stress so 
that a decrease in turbulence kinetic energy is observed between Stations 5 and 
6. Several factors contribute to the attenuation of the turbulence in the near- 
wall region. First, the primary rate-of-strain {dU /dy 2 ), and hence the wall shear 
stress, is reduced by the common outward flow of the vortex pair. Second, lateral 
convergence of the boundary layer acts to suppress turbulence. And finally, 
Fig. 5. 12. a shows that between Stations 4 and 6, the near-wall flow is subjected 
to streamwise acceleration which suppresses turbulence generation. The high 
turbulence levels observed in the mid-region of the boundary layer are a result 
of high primary rates-of-strain associated with the velocity deficit. 

The history of the flow along the y=0 mid-plane (semi-minor axis) is signifi- 
cantly different from the flow along the z=0 mid-plane (semi-major axis). In the 
first half of transition, the boundary layer is subjected to destabilizing concave 
curvature and destabilizing lateral divergence. In the second half of transition, 
the divergence decreases and the flow experiences stabilizing convex curvature. 
Normal stress distributions measured along the duct semi-minor axis at Stations 
5 and 6 are shown in Fig. 5.29 and shear stress and kinetic energy profiles are 
shown in Fig. 5.30. Wall proximity measurements at Stations 5 and 6 were lim- 
ited by the bevel gear arrangement necessary for probe rotation. This limitation 
was most serious along the upper and lower duct surfaces at Station 5 where 
the boundary layer is thinnest. Over the range of the boundary layer that was 
able to be measured, the net effect of the above flow conditioni ng is relatively 
small. The trends show an increase in the axial normal stress component and a 
decrease in the transverse normal stress components between Stations 1 and 6. 
The combined effect on the turbulence kinetic energy is very small. At Station 
5, the shear stress value nearest the wall is observed to decrease, although it 
is recognized that one point is not statistically significant. To investigate this 
further, turbulence measurements along the 3/3 axis were made at a bulk operat- 
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ing Reynolds number of 88,000 where the boundary layer is thicker and a larger 
region could be resolved. Shear stress and kinetic energy profiles for this case are 
shown in Fig. 5.31. These results confirm both the unchanged kinetic energy lev- 
els and the decreasing behavior of the shear stress. The decrease in shear stress 
is in contrast to the increased wall shear stress measured by means of the Preston 
tubes. This indicates that either the shear stress recovers in the near- wall region 
or that the flow is not in local equilibrium near the wall. Without the benefit of 
turbulence measurements at the intermediate Stations 3 and 4, it is difficult to 
speculate on the factors which lead to the flow condition at Station 5. Near-wall 
behavior will be examined more closely in Section 5.7. 

5.6 Turbulence Modelling Considerations 

The inviscid calculations by Burley et al [7,8] have shown that even for 
overall performance predictions, viscous effects cannot be neglected when calcu- 
lating transition duct flows. For accurate predictions the effects of turbulence 
must also be included. The level of turbulence modelling required, of course, 
depends on the information desired. From a computational standpoint, the most 
difficult aspect of the present configuration is the flowfield in the neighborhood of 
the vortex pair. Since the initial generation of streamwise vorticity via pressure 
gradient effects is essentially an inviscid process, even simple models will predict 
the presence of the vortex pairs. Accurate prediction of the diffusion of vortic- 
ity, however, relies on accurate modelling of the Reynolds stress components. 
Modelling efforts for the present configuration can be divided into two groups. 
The first group is concerned primarily with overall performance parameters and 
prediction of flow separation. The computational work of Liandrat ei al. [33] has 
shown that the primary features of the mean flowfield in the region of embedded 
vortices can be predicted reasonably well with simple mixing length models. The 
second group is concerned with the task of demonstrating the ability of particular 
turbulence models to predict detailed mean flow and Reynolds stress behavior. 
The present data set should be useful for this group. 

The simplest turbulence model is the algebraic (zero-order) eddy viscosity 
model based on the concept of Boussinesq. Here, the Reynolds stresses are 
assumed to behave like the molecular viscosity stresses and the Reynolds stress 
tensor for incompressible flow is written as: 


,dUi dUj. 2. , 

+ + 3 Siik 


(5.4) 


where u t is the turbulent eddy viscosity, 6 i} is the Kronecker delta function 
and k is the turbulence kinetic energy. Equation (5.4) is applicable to three- 
dimensional flows where only the Reynolds shear stresses (u,Uj, i ^ j) are 
important, inasmuch as the Reynolds normal stresses predicted by this equation 
are not in agreement with even simple flows. For the present configuration, 
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equation (5.4) may be adequate for prediction purposes through the end of the 
transition section, but will not be applicable if the flow is allowed to develop in the 
rectangular duct where the transverse Reynolds normal stresses are important 
in the generation and diffusion of corner-generated secondary flows. 

Since many design-oriented flow solvers are implemented with algebraic eddy 
viscosity turbulence models, it is worthwhile to examine the behavior of the 
eddy viscosity for the present flow. For the primary shear stresses, if streamwise 
gradients are neglected, equation (5.4) reduces to: 

- dU tK rx 

uv = —Vt — (5.5) 

oy 

— 9U 

u w = —vt~ 7 r~ (o.o) 

Oz 


Rearranging these equations to solve for the eddy viscosity yields: 


— K dV \ 

V '-’ = 

—n du \ 


(5.7) 

(5.8) 


where the y and z subscripts admit directional dependence (anisotropy) of the 
eddy viscosity. The behavior of the primary rates-of-strain (denominators in 
equations (5.7) and (5.8)) at Station 5 is shown in Fig. 5.32. Also indicated in 
these plots are horizontal traverses along which the component eddy viscosities 
were calculated. These results, plotted in terms of a viscosity ratio (u t /u), are 
shown for Station 5 in Fig. 5.33. Equivalent results for Station 6 are shown 
in Figs. 5.34 and 5.35. To avoid singular points where the strain-rates vanish 
(between solid and dashed lines in Figs. 5.32 and 5.34), the eddy viscosity is 
computed only where the denominator in equations (5.7) and (5.8) is greater 
than 5% of its maximum value in the cross plane. Although the eddy viscosity 
distributions shown in Figs. 5.33 and 5.35 show similar trends with respect to 
symmetry about the midplane z=0, the magnitude of the differences is enough so 
that the results should be considered only qualitative. In areas where the stress 
and strain in both transverse directions are nominally of the same magnitude, 
the eddy viscosity components axe observed to be nearly equal. Most of the 
regions where large deviations occur can be traced to either a large difference 
in the component strain rates or to inadequate resolution of a large velocity 
gradient, e.g., in the vicinity of the velocity ridge near the vortex core, this flow. 
This is in contrast to the results of Mehta and Bradshaw, who reported that the 
i / t y component in their embedded vortices was so ill-behaved so els to preclude 
plotting. 
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It was mentioned in Chapter 2 that boundary layer flows with extra rates-of- 
strain often exhibit more spectacular behavior than what is predicted by explicit 
terms which account for the extra strain rates. It is argued that the reason for 
this is that the extra rates-of-strain cause large changes in higher-order terms 
which appear in the Reynolds stress transport equations. These higher-order 
terms are modelled in terms of the Reynolds stresses and require the specification 
of empirical constants. Dimensionless turbulence structure parameters such as 
the cross-correlation coefficient are related to constants in turbulence models. 
Often the constants are determined based on the results of simple flows such as 
two-dimensional boundary layers, resulting in poor performance when applied 
to more complex flowfields. Some of these dimensionless structure parameters 
were computed for the present flow. The shear stress correlations Ruv — uv/u v , 
R uw = uw/u'tv' and R vw = vw/v'w' evaluated at Stations 5 and 6 are shown 
in Figs. 5.36, 5.37 and 5.38, respectively. The primary shear stress correlations 
R uv and R uw may be compared with the two-dimensional value of 0.45. Near 
the duct sidewalls and in the vicinity of the vortex cores, the R uv parameter is 
depressed due to strong lateral convergence of the boundary layer. At Station 
5, in the region where the peak turbulence intensities occur, a near constant 
value of approximately 0.55 is observed. By Station 6, the value has decreased 
and is close to the two-dimensional value. Near the upper and lower walls, the 
Ruw parameter generally is in agreement with the two-dimensional value. The 
secondary shear stress correlation R vw = vw/v'w' at both stations shows values 
in the range of 0.25 for most of the flowfield, but this parameter is also depressed 
in the region of the vortex core. Another dimensionless parameter of interest 
is the shear-energy ratio parameter al. This parameter is defined as the ratio 
of the resultant shear stress in a plane normal to a wall surface to twice the 
turbulence kinetic energy. In a two-dimensional boundary layer, this parameter 
has also been observed to be constant with a value of 0.15 (see, e.g ., Bradshaw 
ct al. [56]). For the present flow, two shear-energy parameters were computed: 

aly = Vut; 2 + vw 2 /2k (5-9) 

al* = \/uu: 2 + vw 2 /2k (5.10) 

which are applicable away from the corner region on the vertical and horizontal 
walls, respectively. The aly shear-ratio parameter at Stations 5 and 6 is shown 
in Fig. 5.39 and the al* parameter is shown in Fig. 5.40. In the regions where 
high turbulence levels are observed, the aly parameter is higher than the two- 
dimensional value and, like the shear stress correlation, is well below the initial 

value in the region of the vortex pair. Along the horizontal walls, the al* param- 

eter is elevated at Station 5, but returns to the initial value at Station 6. These 
results show that the transition duct produces a distortion of the turbulence 
structure, more so at Station 5 than at Station 6. The largest distortion occurs 
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in the vicinity of the vortex core region where lateral convergence suppresses 
turbulence. 

Further insight into the structure of the turbulence can be gained by exam- 
ining the data in terms of the invariants of the anisotropic stress tensor. This 
analysis is based on the concept of physical realizability limits of turbulence 
which has been used extensively to study the return to isotropy of homogeneous 
turbulence (see, e.g., Lumley [57]). The anisotropic stress tensor, first proposed 
by Rotta [58], is defined as: 


t u i u j ~ z^Sij 

li > = 2k 

This tensor must satisfy the Cayley-Hamilton theorem: 

A 3 - IX 2 + IIX - III = 0 
where IJI and III are the tensor invariants: 


(5.11) 


(5.12) 


Ill = 

The first invariant, 7, is identically zero by definition of the turbulence kinetic en- 
ergy. Invariant III defines the shape of the ellipsoid associated with the Reynolds 
stress tensor. A positive value of III indicates that there is only one principle 
component that is large, and a negative value indicates that two principle com- 
ponents are large. Limits on the anisotropic stress tensor can be defined by 
applying the condition that the Reynolds stress in any direction must go to zero 
as the strain rate in that direction goes to infinity. For example, as dU / dx goes 
to infinity, u 2 goes to zero and bn goes to -1/3. The largest level of turbulence 
that can occur in any one direction is 2 k, which occurs when the turbulence is 
one-dimensional. If all the turbulence is in the x direction, then u 2 = 2k and bn 
is equal to 2/3. Following Lumley [57], the limits on allowable turbulence are 
recast in terms of the tensor invariants II and III. These limits are illustrated 
in Fig. 5.41. The shaded area on this plot represents the region within which all 
physically realizable turbulence must lie. The nature of the turbulence at the 
boundaries is also labeled. The Reynolds stress tensor at Stations 1, 5 and 6 was 
recast in terms of the anisotropic stress tensor invariants. At Station 1, the tur- 
bulence is primarily contained in the axial component so it is expected that the 
III invariant will always be positive. The invariants at Station 1 are plotted in 
Fig. 5.42. a. Following Pauley and Eaton, the location of each data point on the 
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invariant map relative to its position in the flowfield is represented by a vector 
which has its origin at II = III = 0 and the tip at the data point location on the 
invariant map. This vector is plotted with its origin coincident with the physical 
location in the flowfield in Fig. 5.42.b. As expected, all the points correspond 
to positive III values. At the edge of the boundary layer the turbulence should 
approach an isotropic condition (zero vector length). The large vectors near the 
boundary layer edge indicate that the turbulence is approaching a more one- 
diroensional nature which is not realistic. This should not be taken to seriously 
though, since the magnitude of the vector is not related to the magnitude of the 
turbulence. That is, experimental uncertainty in regions of low-turbulence levels 
mn cause one component to be large relative to another resulting in an apparent 
high level of anisotropy. Invariant map plots of the Reynolds stress data at Sta- 
tions 5 and 6 are shown in Figs. 5.43 and 5.44, respectively. All of the data are 
observed to lie within the limits of realizability. In the region of the vortex core, 
the turbulence is nearly isotropic. In regions where the primary strain rates are 
large, the turbulence structure is much like that of the initial two-dimensional 
boundary layer at Station 1. Along the velocity ridge, where the primary strain 
rates are small, the III invariant is negative. Between Stations 5 and 6 the data 
cluster more toward and spread further out along the axisymmetric contraction 
limit indicating that the axial stress component is becoming more dominant. 

5.7 Wall Function Behavior 

The use of wall functions for predicting the present flow was analyzed at 
Stations 5 and 6. Reference coordinates for this analysis are shown in Fig. 5.45. 
Since flow variables in the near- wall region scale with the local friction velocity, 
peripheral skin friction coefficient distributions were obtained at Stations 5 and 
6. These distributions, which were deduced from Preston tube data measured 
along the periphery of the duct in Quadrant 1 (see Fig. 5.45), are shown in Fig. 
5.46. The results based on different diameter Preston tubes at Station 5 show 
systematic variations with a change in diameter, particularly for $/ s re f < 0.6, in- 
dicating that deviations from the law-of-the-wall exist at this station. Although 
only two Preston tubes were used at Station 6, the relatively good agreement 
suggests that the flow recovers to law-of-the-wall behavior downstream. Bound- 
ary layer profiles were measured along the y n traverses (n=l,2,5,6,7,8) shown in 
Fig. 5.45 with both pitot and hot-wire probes. Results obtained at Stations 5 
and 6, plotted in law-of-the-wall coordinates, are shown in Fig. 5.47. The fric- 
tion velocities used for these plots are based on the average distribution lines 
shown in Fig. 5.46. The non-dimensional outside diameters of the individual 
Preston tubes are superimposed on these plots. At Station 5, deviations from 
the law-of-the-wall are apparent for y + values greater than 80. In the near- wall 
region, however, all profiles tend toward law-of-the-wall behavior within the in- 
terval 30 < y + < 80. At Station 6, the law-of-the-wall is satisfied on all traverses 
for y + values between 30 and 200. This behavior implies that the law-of-the-wall 
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may be used for prediction purposes, recognizing that the limiting y + value on 
the first mesh line between the duct inlet (Station 1) and Station 5 may be less 
than the limiting value which applies at Station 5 (y + — 80). 

Wall function behavior for the turbulence kinetic energy ( k ) and the dissi- 
pation rate of turbulence kinetic energy (e) was also analyzed at Stations 5 and 
6. If only first-order effects are retained and terms containing secondary rates- 
of-strain are neglected in the reduced forms of the Reynolds stress transport 
equations, and the mean velocity profile along any normal to the duct surface 
is in accordance with the law-of-the-wall, then the following expressions (rela- 
tive to the x w ,y w ,z w wall coordinate system) apply for turbulent flow through 
non-circular ducts in local equilibrium [59]: 


where 
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with C p = 0.09 and « = 0.41, and where l p is Prandtl’s length scale. In the near- 
wall region along the transition duct horizontal and vertical walls at Stations 5 
and 6, ( dU/dy w ) > ( dU/dz w ), so that F p = 1.0. Further, it is reasonable to 
assume that along these walls, but excluding the comer region, Prandtl’s length 
scale can be specified as l p =Ky w . With these assumptions, the following alternate 
forms of equations (5.13)-(5.15) apply along the horizontal and vertical walls: 


(5.17) 

(5.18) 

(5.19) 

where e has been replaced by P (prod uction) in a ccordance with the local equilib- 
rium assumption. Distributions of Vuv 2 + uw 2 /U 2 , k/U 2 and Py w /Ul along 
several traverses parallel to the z axis at Stations 5 and 6 are shown in Figs. 5.48, 
5.49 and 5.50, respectively. Also indicated on these plots are the conventional 
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wall function values. Although the results are restricted to y values greater 
than 400 because of limitations imposed by the probe rotation mechanism, some 
conclusions can still be drawn. For example, it can be seen that distributions 
measured along the vertical midplane traverse y/R = 0 at Stations 5 and 6 gen- 
erally tend toward the conventional limits. This same comment applies for all 
other distributions measured at Station 6. At Station 5, however, the trends are 
more difficult to discern, especially along the y/R = 1.37 traverse, which lies 
within the sidewall boundary layer. 

The behavior of Vuv 2 + uw 2 /U^, k/Ul and Py/Uf along the semi-major 
traverse y 2 = 0 at Stations 5 and 6 is shown in Figs. 5.51, 5.52 and 5.53, respec- 
tively. The distributions measured at Station 6 approach the conventional limits. 
The deviation observed at Station 5 is probably attributable to the neglect of 
secondary rates-of-strain in the derivation of equations (5.13)-(5.15), inasmuch 
as Fig. 5.9 shows that secondary strain near the sidewall at Station 5 is a large 
percentage of the primary rate-of-strain, but that it decreases to a relatively 
small percentage by Station 6. The first data point for Station 5, however, is at 
y+ ^ 300 which is beyond the region where law-of-the-wall behavior is observed 
(refer to the yf distribution of Fig. 5.47.a). In reference to predictions, if y + 
on the first mesh line is restricted to an interval between 30 and 50, and if sec- 
ondary rates-of-strain are small within this interval, then wall function values on 
this mesh line may be close to the conventional limits at Station 5 and, indeed, 
at intermediate stations within the duct. Further measurements in the near-wall 
region will be required in order to verify this conjectured behavior. 
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a) Re b = 88, 000. 
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b) Re b = 390, 000. 

Fig. 5.1. Peripheral wall pressure coefficient distributions. 
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Fig. 5.4 
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c) Station 5. 
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Fig. 5.5. Total pressure contours ((Pt — Pj)/Qref)- 




.5 


1.0 


1.5 


y/R 

I Station 3. 
contours (U/Ub). 


(Continued . . . ) 











76 



0.0 0.5 1.0 1.5 

y/R 


a) Station 3. 

Fig. 5.7. Transverse mean velocity vectors (V r /Ub). (Continued . . . ) 
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b) Station 4. 

Fig. 5.7. Transverse mean velocity vectors (V r /Ub). (Continued 
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d) Station 6. 

Fig. 5.7. Transverse mean velocity vectors (V r /Ub) 
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b) Calculated. 

Fig. 5.8. Static pressure contours ((P — Pj)/Qref) at Station 5 
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Fig. 5.9. Axial distribution of divergence 
rate-of-strain parameter along midplanes y=0 and z=0. 





81 



a) y 2 traverse (see Fig. 3.2). 



b) y 3 traverse (see Fig. 3.2). 

Fig. 5.11. Total pressure profiles along the 
semi-major and semi-minor axes, Reb = 390, 000. 
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b) Station 4. 

Fig. 5.13. Axial mean vorticity contours (ft x R/Ub). (Continued . . . ) 






d) Station 6. 

Fig. 5.13. Axial mean vorticity contours (fl x R/Ub) 
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b) Station 6. 

Fig. 5.14. Axial turbulence intensity contours (u'/U x 10 2 ). 
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b) Station 6. 

Fig. 5.15. u^/Uj x 10 4 Reynolds stress contours. 

c-a 
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b) Station 6. 

Fig. 5.16. v 2 /U£ x 10 4 Reynolds stress contours. 
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b) Station 6. 

Fig. 5.20. vw/U£ x 10 4 Reynolds stress contours. 
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Fig. 5.23. Transverse normal stress anisotropy ^v 2 - w j )/ug x 10 4 ). 
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a) Station 5. 
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Fig. 5.24. Turbulence kinetic energy (k/Uj; x 10 4 ). 
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Fig. 5.26. Reynolds normal stress profiles at Station 1. 
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b) Station 6. 

Fig. 5.27. Normal stress profiles along semi-major axis at Stations 5 and 6 
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a) Shear stress. 
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b) Turbulence kinetic enegy. 

Fig. 5.28. Shear stress and kinetic energy profiles along semi-major axis. 
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a) Station 5. 
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b) Station 6. 

Fig. 5.29. Normal stress profiles along semi-minor axis at Stations 5 and 6 
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a) Shear stress. 
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b) Turbulence kinetic enegy. 

Fig. 5.30. Shear stress and t.k.e. profiles along semi-minor axis, Reb= 390, 000 
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a) Shear stress. 
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b) Turbulence kinetic enegy. 

Fig. 5.31. Shear stress and t.k.e. profiles along semi-minor axis, Reb= 88,000. 
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Fig. 5.32. Primary rates-of-strain at Station 5. 
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b) Quadrant 2. 

Fig. 5.33. Component eddy viscosity distributions at Station 5. 
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Fig. 5.34. Primary rates-of-strain at Station 6. 
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b) Quadrant 2. 

Fig. 5.35. Component eddy viscosity distributions at Station 6. 





b) Station 6. 

Fig. 536. Shear stress correlation R uv = uv/uV. 





b) Station 6. 

Fig. 5.37. Shear stress correlation R w = uw/uV. 
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b) Station 6. 

Fig. 5.38. Shear stress correlation R vw = vw/v'w'. 
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Fig. 5.40. Shear-energy parameter al z = \/iiw 2 + vw 2 /2k. 
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Fig. 5.41. Physical realizability limits of turbulence invariant map. 
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a) Invariant map. 
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b) Location of data points on invariant map. 

Fig. 5.42. Anisotropic stress tensor invariants at Station 1. 
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a) Invariant map. 



b) Location of data points on invariant map. 

Fig. 5.43. Anisotropic stress tensor invariants at Station 5. 
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Fig. 5.44. Anisotropic stress tensor invariants at Station 6. 
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Fig. 5.45. Reference coordinates for wall function behavior. 



Fig. 5.46. Periperal skin friction coefficient 
distributions at Stations 5 and 6, Reb= 390,000 . 
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b) Station 6. 

Fig. 5.47. Law-of-the-wall velocity profiles at Stations 5 and 6, Reb= 390,000 . 
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Fig. 5.50. Py w /Ur distributions parallel to semi-minor axis. 
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Fig. 5.52. k/U? distributions along semi-major axis. 
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Fig. 5.53. Py w /U? distributions along semi-major axis. 




CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 
6.1 Conclusions 

Incompressible, turbulent, swirl-free flow through a circular- to-rectangular 
transition duct has been studied experimentally. Mean flow measurements are 
presented at five axial stations and the complete Reynolds stress tensor is pre- 
sented at three axial stations for operating bulk Reynolds numbers of 88,000 and 
390,000. No significant Reynolds number dependence was observed and excellent 
symmetry prevailed through the entire test section. The results of the study may 
be summarized as follows: 

1) The inlet conditions to the duct, which correspond to partially developed 
turbulent pipe flow, proved to be an interesting case, in itself, considering 
the large differences between results of various studies on the subject. The 
present results showed that the boundary layer was significantly thinner 
than that reported by other investigators over nearly equivalent development 
lengths (x/D fa 9). Reynolds stress levels also differed. These discrepancies 
are generally attributable to inlet conditions at the pipe entrance. Excellent 
law-of- the- wall behavior was observed in the pipe and the skin friction was 
very near the fully- developed value. 

2) The results of the transition duct measurements show that streamwise vor- 
ticity generated in the first half of transition strengthens and concentrates 
along the duct sidewall in the second half of the transition resulting in a 
contra-rotating vortex pair (common flow away from surface) centered about 
the duct semi-major axis. The upwash action between the vortices creates 
a region of low wall shear stress. The vortex pair significantly distorts the 
mean flow and turbulence fields. 

3) In the transition section, transverse flow velocities of approximately 32% of 
the inlet bulk velocity were observed. In the constant area duct extension, 
the vortex pair causes a secondary flow to persist that peaks at about 12% 
of the inlet bulk velocity. 

4) In the vicinity of the sidewall vortices, a large velocity deficit exists. The 
large primary rates-of-strain associated with this velocity deficit dominate 
much of the vortex region. These large strain rates produce a region of 
relatively high turbulence intensity. 

5) In the vicinity of the vortex cores, nearly isotropic turbulence exists which 
is dramatically attenuated by lateral convergence effects. 
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6) Along the horizontal walls at the end of transition, a peculiar trend in the 
primary turbulent shear stress was observed which is a result of surface 
curvature and lateral divergence effects. This anomaly was manifested in 
a dramatic reduction in the shear stress midway across the boundary layer 
as the wall is approached. The measurements did not, however, resolve the 
near-wall region to determine if the shear stress recovered in the inner layer. 

7) At the end of transition, law-of-the-wall behavior indicates that the flow is 
in local equilibrium over only a very small region near the wall. Trends in 
turbulence-based wall function behavior at this station tend to support the 
local equilibrium assumption, but the present measurements do not fully 
resolve the near-wall region. At the end of the transition duct extension, 
local equilibrium is observed over a fairly large region near the wall. 

The present results can serve els an experimental database for CFD code calibra- 
tion/verification, inasmuch as the results demonstrate the appropriate symmetry 
and are based on well-defined inlet conditions. The results of this study should 
also provide a stringent test for turbulence models. To facilitate comparisons, 
the transition duct geometric vsiriables are available in Appendix A. Tabulated 
flow data in machine readable form are available from the author. 

6.2 Recommendations 

Between the time that the present study was initiated Eind completed, a 
number of studies dealing with circular-to-rectangular transition ducts appeared 
in the literature. Most notable Eire those of Reichert et al. [4] and Miau et al. 
[9]. Despite these additions, there are still some areas which have not been fully 
addressed. First among these is the case with inlet swirl. Inlet swirl may be 
beneficial from the standpoint of reduced noise and thermal plume. Although 
Reichert has provided mean flow data at several axial stations in a CR transition 
duct, a large region in the vicinity of the sidewsdl was not accessible. Supplemen- 
tary mean flow data in this region and turbulence measurements in the entire 
cross-plane will provide a more complete picture for the swirling case. Another 
area of interest is the case of laminar flow through a CR transition duct. This 
is of interest because laminar flow data ELre useful for code vsdidation purposes. 
Having both laminar and turbulent data allows the code developer to distinguish 
between deficiencies in the numerical scheme and in the turbulence model. Fi- 
nally, measurements in a different geometry transition duct may be useful. In 
particular, investigation of an identicEil duct to the one used in the present study, 
but without the area variation, may give further insight to the development of 
streamwise vortices, inasmuch as it has been observed that ducts which have 
no area expansion are devoid of vortices. Discrete vortices have been observed, 
however, in ducts which have an area expansion, for different aspect ratios and 
transition lengths. All of the above research extensions could be made in the 
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present facility with minimal hardware change. The swirling case requires only 
the addition of turning vanes in the inlet pipe or, preferably, a swirl generator 
similar to the one proposed by Reichert which very nearly approximates a solid- 
body rotation inlet condition. For the laminar flow case, in addition to removing 
the boundary layer trip, a means for reducing the flow rate through the facility 
would be required to ensure that the flow remained laminar through the entire 
test section. This could be accomplished by providing a controllable bypass flow 
downstream of the test section. 
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APPENDIX A 


TRANSITION DUCT GEOMETRY 

At any given streamwise station, the cross section of the transition duct can 
be described by the equation of a superellipse, namely: 

( y/a .)» + ( z/b)' = 1 (A.l) 

where a and b are the semi-major and semi-minor axes, respectively, of the su- 
perellipse. Special cases of the superellipse include a circle (a = 6, r) = 2), an 
ellipse (a / b, rj = 2) and a rectangle (a^b,rj = oo). The cross-sectional area at 
a given streamwise location is given by [60]: 

A„ = ^£(2/,)(4a&) (A2) 

where T refers to the “gamma function” , defined as: 

r oo 

r(r?)= / (17 > 0) (A3) 

Jo 

The design procedure for the transition duct is to specify the streamwise vari- 
ation of the semi-major and semi-minor axes (which define the side and upper 
(lower) walls of the superellipse) and the streamwise variation, if any, of the 
cross-sectional area. Once a(x), 6(x) and A cs (x) are known, equation (A. 2) is 
then solved numerically to determine the streamwise variation of rj. 

For the transition duct used in this study, the variation of the semi-major and 
semi-minor axes are specified by fifth-order polynomial functions of x/ R. These 
functions provide zero derivatives up to the second order at the beginning and end 
of transition. The aspect ratio (AR) and the length ratio ( L/D ) for the present 
duct are 3 and 1.5, respectively. The cross-sectional area begins to increase 
at the start of transition to a maximum of 15% above the inlet cross-sectional 
area, occurring approximately at the mid-point of transition. From there the 
area smoothly decreases until the end of transition where the cross-sectional 
area is equal to the inlet cross-sectional area. After solving equation (A. 2), the 
variation of r] is then approximated by a fifth-order polynomial function. Table 
A.l contains the polynomial coefficients for the variation of a/R, b/R and r? for 
the transition duct. Table A. 2 contains tabulated values of a/R , b/R, rj, A cs and 
the aspect ratio AR = a/b along the transition duct. Also indicated in Table A.2 
are the locations of the six data stations used in the present study. 
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Table A.l. Transition Duct Polynomial Coefficients 



a/R 

b/R 

n 

Co 

1.000000000 

x/R < 1.0 
1.000000000 

2.000000000 

C\ 

0.000000000 

0.000000000 

0.000000000 

C2 

0.000000000 

0.000000000 

0.000000000 

c$ 

0.000000000 

0.000000000 

0.000000000 

c 4 

0.000000000 

0.000000000 

0.000000000 

c 5 

0.000000000 

0.000000000 

0.000000000 

Co 

0.87568217 

1.0 < x/R < 1.96 
1.1102199 

8.4833465 

C\ 

0.38699868 

-0.34310201 

-21.329880 

c 2 

-0.38783144 

0.34380721 

26.148548 

Cz 

0.096543991 

-0.085532496 

-15.173079 

c 4 

0.041969148 

-0.037245857 

4.3614143 

C5 

-0.013362557 

0.011853259 

-0.49034940 

Co 

-0.37336560 

1.96 < x/R < 2.92 
2.2162779 

1.0971109 

C 1 

3.3510193 

-2.9683112 

2.3623579 

c 2 

-3.2796590 

2.9056348 

-2.8785242 

Cz 

1.5523217 

-1.3754753 

2.0620912 

C 4 

-0.33664930 

0.29831302 

-0.64136301 

c$ 

0.027262234 

-0.024157751 

0.080685615 

Co 

14.292266 

2.92 < x/R < 4.125 
-10.787227 

-1858.4412 

Cl 

-19.643522 

17.419488 

2828.7353 

C2 

11.204615 

-9.9360792 

-1704.8063 

c 3 

-3.0380988 

2.6941403 

508.59846 

c 4 

0.39689283 

-0.35195770 

-74.960370 

C5 

-0.020128978 

0.017849940 

4.3651652 

Co 

1.5463591 

x/R > 4.125 
0.51544988 

10.000000 

Cj 

0.0000000 

0.00000000 

0.0000000 

c 2 

0.0000000 

0.00000000 

0.0000000 

c 3 

0.0000000 

0.00000000 

0.0000000 

c 4 

0.0000000 

0.00000000 

0.0000000 

C5 

0.0000000 

0.00000000 

0.0000000 


R = 10.214 cm 
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Table A. 2. Transition Duct Geometry Variables 
Data 

Station x/R a/R b/R r] A R Act/Act^ 

1 - 1.00 1.0000 1.0000 2.0000 1.0000 1.0000 

H ^ ^ ^ ^ 

2 1.00 1.0000 1.0000 2.0000 1.0000 1.0000 

1.12 1.0003 0.9997 2.1047 1.0006 1.0192 

1.24 1.0025 0.9978 2.2154 1.0047 1.0376 

1.36 1.0078 0.9831 2.3331 1.0251 1.0555 

1.48 1.0174 0.9746 2.4588 1.0439 1.0729 

1.60 1.0316 0.9719 2.5934 1.0614 1.0897 

1.72 1.0509 0.9548 2.7383 1.1006 1.1056 

1.84 1.0752 0.9333 2.8947 1.1520 1.1198 

1.96 1.1041 0.9177 3.0645 1.2031 1.1319 

2.08 1.1371 0.8784 3.2497 1.2945 1.1410 

3 2.20 1.1734 0.8462 3.4528 1.3867 1.1464 

2.32 1.2123 0.8117 3.6770 1.4935 1.1475 

2.44 1.2527 0.7759 3.9261 1.6145 1.1441 

2.56 1.2936 0.7396 4.2050 1.7491 1.1362 

2.68 1.3341 0.7037 4.5199 1.8958 1.1241 

4 2.80 1.3729 0.6693 4.8791 2.0512 1.1087 

2.92 1.4093 0.6370 5.2932 2.2124 1.0911 

3.04 1.4423 0.6078 5.7775 2.3730 1.0725 

3.16 1.4712 0.5821 6.3524 2.5274 1.0544 

3.28 1.4954 0.5606 7.0484 2.6675 1.0381 

3.40 1.5147 0.5435 7.9106 2.7869 1.0248 

3.52 1.5290 0.5308 9.0131 2.8806 1.0152 

3.64 1.5385 0.5224 10.000 2.9451 1.0083 

3.76 1.5439 0.5176 10.000 2.9828 1.0026 

3.88 1.5460 0.5157 10.000 2.9979 1.0003 

5 4.00 1.5464 0.5154 10.000 3.0000 1.0000 

4J. JJ. JJ. -li 4 ^ 

8.00 1.5464 0.5154 10.000 3.0000 1.0000 
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APPENDIX B 


HOT-WIRE DATA ACQUISITION SYSTEM 

An all new integrated hot-wire data acquisition system was acquired for the 
present study. This system differs from previous ones used within the M.E. 
department in that, rather than operating directly on the analog anemometer 
output, data reduction operations are performed on a digital representation of 
the analog signal. This appendix contains a description of the system and an 
overview of some of the more important operating features. 


B.l Hardware 

The primary components of the system are a Thermal Systems Incorporated 
(TSI) IFA 100 Intelligent Flow Analyzer, a TSI IFA 200 analog-to- digital (A/D) 
converter and a Digital Equipment Corporation (DEC) VAXLab/VSII Scientific 
Workstation computer. A block schematic of this equipment is illustrated in Fig. 
B.l. 

The IFA 100 used in this study is a two-channel unit consisting of a Model 
158 slave cabinet housing two Model 150 constant temperature anemometers and 
two Model 157 signal conditioners. The signal conditioners contain a low-pass 
filter, a high-pass filter and an amplifier capable of applying a DC offset (when the 
high-pass filter is in the DC mode) and gain to the input (anemometer output) 
signal. The usable range of the Model 158 digital display and the input voltage 
range of the IFA200 digitizer is ±5 volts. Most wire and film sensors will have 
a voltage span much less than this. To utilize as much of this range as possible, 
the transducer signal must be offset and amplified. 

The IFA 200 provides high-speed simultaneous A/D conversions of one or 
more anemometer output signals. Selectable A/D conversions rates of up to 
50,000 Hz per channel provide the capacity of collecting data on rapidly changing 
signals. The computer controls the IFA 200 by sending ASCII characters through 
an RS-232 serial communication link. The parameters that are controlled through 
this link are the starting address (channel to be sampled first), the number of 
active channels and the sampling rate. For each sample, the A/D converter 
generates a 16-bit data word. The low 12-bits contain the binary representation 
of the transducer voltage. The high 4-bits contain the digitizer channel address 
which must be stripped from the 16-bit data word before converting the binary 
number to a hexadecimal voltage. The 16-bit data word is loaded into a 64- 
word-deep first-in, first-out (FIFO) memory. The data that enter the FIFO pass 
to the FIFO output and are sent via a DMA (Direct Memory Access) interface 
cable to a DRQ3B DMA controller card in the computer. 
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The DRQ3B Parallel DMA I/O module is a high-speed memory interface 
that allows real-time collection of 16-bit parallel digital data at transfer rates 
of up to 1.3 megaHz. The DRQ3B has a unique feature, referred to as double 
buffering , that can significantly increase data transfer rates as compared to con- 
ventional single buffer DMA transfers. With single buffer DMA, when a buffer 
fills, it is dequeued, but requires a software call to enqueue the next buffer. Dur- 
ing the time between dequeueing the previous buffer and enqueueing the next 
buffer, no data are transferred. This time between buffers must not exceed the 
time it takes to fill the 64-word-deep FIFO memory in the IFA200, or data will 
be lost. To improve this situation, the DRQ3B allows two or more buffers to be 
enqueued before any data are transferred. Now, when a buffer fills, the DRQ3B 
automatically begins filling the next buffer without any software intervention. 
The result is less down time between buffers, enabling higher transfer rates. 

B.2 Buffer-Swinging 

Transferring data at high sampling rates can quickly exceed the available 
memory in the computer. To avoid this, the digitized data Eire forwarded directly 
to the 71 megabyte hard disk in the computer by means of the buffer- swinging 
method of data acquisition. This feature is used to capture a continuous stream 
of data and operates in the following manner. Two or more memory buffers are 
enqueued to the DRQ3B prior to any data transfer. When the last buffer is 
enqueued, data transfer begins to the first buffer. When the first buffer fills, the 
DRQ3B automatically begins filling the second buffer. While the second buffer 
is filling, the data in the first buffer are written directly to the hard disk. When 
the second buffer fills and first buffer’s data have been transferred to disk, the 
process repeats as the second buffer’s data are sent to disk and the first buffer 
again receives data. 

B.3 Handshaking 

Synchronization of data transfer between the IFA 200 and the DRQ3B con- 
troller card is accomplished by an interlocked two-wire handshake. When prop- 
erly implemented, the handshake prevents data overrun and underrun conditions 
from occurring. The FUNCT OUT 0, STROBE and ACK pins on the DRQ3B 
input port are connected, via the DMA Interface Cable, to the READY, CYCLE 
REQUEST and BUSY pins, respectively, on the IFA 200 output port, f For data 

f At the time the data acquisition equipment was acquired, TSI did not sup- 
port DMA transfers from the IFA 200 to the DRQ3B controller card. To accom- 
plish this, the DMA cable had to be modified for compatibility with the DRQ3B. 
The pin-to-pin connections between the J1 and J2 output connectors on the IFA 
200 DMA Interface Cable and the J2 input port on the DRQ3B DMA card are 
listed in Table B.l. 
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transfer to begin, the computer asserts the FUNCT OUT 0 pin on the DRQ3B 
high. When the IFA 200 senses this on the READY line, it turns on the sampling 
strobe, begins filling its 64- word-deep FIFO memory and asserts the CYCLE RE- 
QUEST line low, indicating to the DRQ3B that valid data are available. Until 
the last buffer is enqueued on the DRQ3B by the computer, the DRQ3B asserts 
the ACK signal low, indicating that it is not ready to receive data. When the 
last buffer is enqueued, the DRQ3B releases the ACK line allowing it to go high, 
indicating that it is now ready to receive data. 

Once the digitizer has been activated and the DMA buffers have been en- 
queued, the following sequence of events occur for data transfer: 

1. When the DRQ3B receives STROBE low from the CYCLE REQUEST line 
of the IFA 200, it reads the data and places the value in its 512-word-deep 
FIFO memory. The DRQ3B then asserts the ACK signal low, indicating to 
the IFA 200 that it has received the data value. 

2. When the IFA 200 receives BUSY low from the ACK line of the DRQ3B, it 
sets the CYCLE REQUEST line high, indicating to the DRQ3B that valid 
data are not available, and prepares to transmit the next data word. 

3. When the DRQ3B receives STROBE high from the CYCLE REQUEST line 
of the IFA 200, it releases the ACK line allowing it to go high (unless the 
DRQ3B FIFO is full, in which case ACK is held low until a word is read out 
of FIFO and into memory). ACK high indicates that the DRQ3B is ready 
for the next word. 

4. When the IFA 200 receives BUSY high from the ACK line of the DRQ3B, 
it places valid data on the data lines and asserts CYCLE REQUEST low, 
starting the cycle over. 

B.4 Software 

Control of the DRQ3B is accomplished through the use of DEC’s VAXlab 
subroutine library. This software package provides useful programs for all aspects 
of the experimental study, including, real-time data acquisition, signal process- 
ing, plotting and data analysis. The subroutines are called from user supplied 
main programs written in the FORTRAN language. Control of the IFA200 is 
accomplished by sending ASCII characters, via an RS-232 serial communication 
link, by means of the FORTRAN WRITE statement. 

B.5 Digital Sampling 

The hot-wire techniques employed in the present study (see Appendix C) 
rely on an accurate evaluation of the mean and variance (F and e^) of the analog 
output signal from the hot-wire anemometer. The usual assumption that the 
velocity, and therefore the anemometer output voltage, is an ergodic random 
variable is made here. If the instantaneous voltage is decomposed into a mean and 
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fluctuating component, E = E + e, the true temporal average can be estimated 
by a finite average of discrete samples, provided that a sufficiently large number 
of samples N are taken: 

i rio+T 2 N 

E = lim - / E(t)it « -5> n (nA<) (B.l) 

T ^°° 1 Jt ° " n=l 

i r*o+T i ^ 

e 2 = lim - / (£(() - Efdi « - - Ef ( B.2 ) 

T -*°° 1 J* 0 A „=1 

where E n are the individual digitized samples, T is the total sampling time and 
At is the time increment between samples. For a fixed sampling time duration 
T, the number of samples is N — T/At. Since only single statistical quantities 
are required for the hot-wire measurements, it was not necessary to “capture” 
the waveform by digitizing at very high sampling rates. For the turbulence 
measurements, the sampling rate and duration time were set at 10 kHz and 15 
seconds, respectively, which yields a total of 150,000 samples for each data point. 
The convergence of the mean and mean-square values was monitored for each 
data point taken and it was observed that the above transfer rate and sampling 
time were more than adequate for all regions of the transition duct flowfield. 
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Table B.l. DMA Interface Cable Pin Connections 


IFA 200 DRQ3B 


Signal 

Name 

Port 

No. 

Pin 

No. 

Pin 

No. 

Port 

No. 

Signal 

Name 

CYCLE REQUEST 

J1 

02 

17 

32 

STROBE 

READY 

J1 

06 

22 


FUNCT OUT 0 

GND 

J2 

01 

42 


GND 

BUSY 


02 

18 


ACK 

GND 


03 

41 


GND 

GND 


04 

40 


GND 

GND 


07 

39 


GND 

GND 


11 

38 


GND 

GND 


12 

37 


GND 

GND 


13 

36 


GND 

GND 


15 

35 


GND 

GND 


17 

34 


GND 

GND 


19 

32 


GND 

GND 


20 

31 


GND 

GND 


21 

29 


GND 

GND 


22 

28 


GND 

GND 


23 

27 


GND 

GND 


24 

26 


GND 

DATA OUT 00 


39 

01 


DATA IN 00 

DATA OUT 01 


37 

02 


DATA IN 01 

DATA OUT 02 


35 

03 


DATA IN 02 

DATA OUT 03 


33 

04 


DATA IN 03 

DATA OUT 04 


31 

05 


DATA IN 04 

DATA OUT 05 


29 

06 


DATA IN 05 

DATA OUT 06 


27 

07 


DATA IN 06 

DATA OUT 07 


25 

08 


DATA IN 07 

DATA OUT 08 


26 

09 


DATA IN 08 

DATA OUT 09 


28 

10 


DATA IN 09 

DATA OUT 10 


30 

11 


DATA IN 10 

DATA OUT 11 


32 

12 


DATA IN 11 

DATA OUT 12 


34 

13 


DATA IN 12 

DATA OUT 13 


36 

14 


DATA IN 13 

DATA OUT 14 


38 

15 


DATA IN 14 

DATA OUT 15 

J2 

40 

16 

J2 

DATA IN 15 


Note: Unlisted pins are not connected. 


RS-232 Comnand 
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IFA 100 Anaaoneter 



Fig. B.l. Schematic of hot-wire data acquisition system. 




APPENDIX C 


HOT-WIRE TECHNIQUES 


C.l Introduction 

The derivation of the mean and turbulence hot-wire response equations for 
the Method A technique are presented in this appendix. The working form of the 
mean and turbulence response equations developed by Al-Beirutty ( Method B ) 
are presented in this chapter without derivation. For details of the development 
of these equations, the reader is referred to references [45] and [46]. 

The Method A single rotatable hot-wire technique which is proposed here 
is applicable for the measurement of the mean flow and Reynolds stress field 
in situations when the primary flow is nominally normal to the probe axis of 
rotation. This, in itself, is not unique. De Grande and Kool [61] proposed a hot- 
wire technique in which a single rotatable slant-wire is used to deduce the mean 
flow and turbulence field under similar conditions. Their technique is not based 
on an empirical cooling law relation, but rather relies on extensive pitch and yaw 
calibration. The time required for calibration becomes prohibitive when large 
amounts of data are being accumulated and frequent recalibration is necessary. 
Their technique has the advantage that commercially available probes can be 
used. The use a single slant-wire probe limits the range of sensitivity to the 
individual mean velocity and stress components. This is usually compensated 
for by obtaining data at a large number of rotational positions and using least- 
squares error minimization criteria to solve for the flow variables. The hot-wire 
technique proposed here relies on the sequential use of custom-made normal and 
slant-wire probes in order to maximize sensitivity to the flow variables. The 
response of the hot-wires to pitch and yaw flow angles is governed by an effective 
cooling velocity relation so that only a simple calibration procedure is required. 

Before continuing, it is important to note a restriction associated with 
Method A that is not present with Method B. Whereas with Method B, for zero- 
to-moderate skew angles, the probe body can be rotated 360 with little or no 
interference effects from the probe supports, there is a significant region of rota- 
tion associated with Method A where interference effects from the probe supports 
cannot be neglected. Within this region, no valid data can be expected; how- 
ever, the angular sweep of this region can be minimized and the angular position 
shifted by suitable design of the hot-wire probe itself. Two possible slant-wire 
probe configurations, along with an estimate of their regions of interference for 
use with Method A are illustrated in Fig. C.l. The impact of this restriction on 
Method A will be discussed as the response equations are developed. 
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C.2 Generalized Empirical Cooling Law Relation 

The starting point for the derivation of the mean and turbulence response 
equations for both Method A and Method B is the empirical cooling law relation. 
The form of this equation is dependent on whether a linearized or non-linearized 
hot-wire anemometer system is used. For a non-linearized system, the relation 
is given by: 

(E NL ) 2 -E 2 = B NL U c n (C.l.a) 

where E is the instantaneous bridge voltage and U e is the instantaneous cooling 
velocity; E 2 and B are the intercept and slope of the wire calibration curve when 
the calibration stream is normal to the wire. Conversely, when the anemometer 
output is linearized, the cooling relation is given by: 

E l = B L U t (C.l.b) 

These equations can be presented in a general form applicable to both linearized 
and non-linearized signals as: 

(E 2 - A) 2/m = B 2/m U 2 (C.l.c) 

where, for a linearized system: 

E 2 = (E l - Eo) 2 , A = 0, B = B 2 L , m = 2 

and for a non-linearized system: 

E 2 = E 2 nl , A = Eq, B = B nl , m = n 

In the following analysis, the mean and turbulence response equations are devel- 
oped for the non-linearized case, from which the linearized system equations are 
extracted as a special case. 

C.3 Mean Flow Equation Development 

The mean velocity vector can be evaluated in magnitude and direction when 
a hot-wire is rotated into four singular positions, such as Positions 1,2,3 and 4 
shown in Fig. C.2, and the corresponding mean bridge voltage is recorded. For 
Method A, this can be accomplished by rotating an inclined wire about the y- 
axis to Positions 1 and 2, and by rotating a normal wire about the same axis 
to Positions 3 and 4. Without knowing the mean velocity direction, however, 
there is the risk that the wire will be rotated into a position where interference 
from one of the probe supports will occur, rendering the data invalid. We will 
discuss later how this situation can be avoided, if beforehand, an estimate of the 
transverse flow angle in the x-z plane, 7, defined as shown in Fig. C.2, is known. 
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The reader should keep in mind that this is strictly a physical limitation and not 
a restriction inherent in the response equations themselves. 

For Method B, the four wire positions shown in Fig. C.2 can be achieved 
by rotating a single inclined wire about the x-axis. Unlike Method A , for zero- 
to-moderate skewness levels, the problem of probe support interference does not 
present itself with Method B, and, therefore, no knowledge of the mean flow 
direction is required apriori, other than that the flow is nominally in the x~ 
direction and that V and W are less than U . 

The primary difference between Method A and Method B is in the expressions 
for the effective cooling velocity term in equation (C.l.a). In the remainder of 
this section, the mean flow response equations for both Method A and Method B 
will be presented. Except where noted, the following discussion is applicable to 
both of these methods. 

When the time- averaged form of equation (C.l.a) is applied at each wire 
position shown in Fig. C.2, the following four response equations result: 


(£, 2 -£ 0 2 ) 2 / m = B 2/m (C e 2 )i (C.2. a) 

(E 2 2 -^) 2 / m = B 2/m (u7) 2 (C-2.6) 

(£ 3 2 - £ 0 2 ) 2/m = £ 2/m (u7) 3 (<?.2.c) 

(£ 4 2 -£ 0 2 ) 2/m = B 2/m {u7)< (C.2.d) 

Combining equations (C.2. a) & (C.2.b) and equations (C.2.c) & (C.2.d) yields. 

{E? - - £ 0 2 ) 2/m = W)i/Wh ( c - 3 - fl ) 

(E3 2 -4 2 ) 2/Tn /(^4 2 -£ 0 2 ) 2/Tn = ( C.Z.b ) 

Rearranging equations (C.3.a) and (C.3.b) gives: 

W)i - Si2(t?) 2 = 0 ( C-*.a ) 

(U?h ~ ^34<U?)< = 0 (C.4.6) 

where, for convenience, the following definitions are employed: 

E 12 = (E? - E 0 2 ) 2 /"*/(£ 2 2 - £ 0 2 ) 2/m (C.5.a) 

S 34 = (E 3 2 --&o 2 ) 2/m /( £: 4 2 --Eo 2 ) 2/m 


By eliminating the dependence on the wire calibration slope ( B ) in equations 
(C.4.a) and (C.4.b), we have placed the restriction that the same hot-wire must 
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be used in Positions 1 and 2 , and, likewise, in Positions 3 and 4. We do not 
require, however, that the hot-wire that is used in Positions 1 and 2 be the 
same as the one used in Positions 3 and 4, so long as the appropriate calibration 
intercept values are used in equations (C.5.a) and (C.5.b), 

In order to solve equations (C.4.a) and (C.4.b) for the mean velocity com- 
ponents, the effective cooling velocity ( U e ) must be expressed in terms of the 
mean velocity components (U, V, W), and the voltage terms (E 12 , E 34 ) must be 
expressed in terms of measurable mean (E) and mean-square (e 2 ) voltages. 

The voltage terms, and E 34 , in equations (C.4.a) and (C.4.b) can be 
expressed in terms of measurable voltages by substituting E = E + e into (E 2 — 
Eq) 2 /” 1 , expanding the resulting expression in a binomial series, and truncating 
at second order: 

(£ 2 - Elf " = (E 2 - Elf"\ 1 + (C.6.«) 

where: 

C = 2/m -|- (4/m)( ^ - 1)( ^ -=^ ) (C.6.b) 

Time-averaging this expression yields: 

_ _ 

(E* - Elf" = (£ 2 - Elf"[ 1 + JTTeI] (£•«•') 


Substituting equation (C.6.c) into equations (C. 5 .a) and (C.5.b) gives: 


E 12 


S1/S2 


(C.l.a) 


S34 = E3/E4 


where: 


E, = (£? - E 2 ) 2 '™[1 + 

E 2 = (£l-£ 0 ! ) J /"[l + ^lj] 

E s = (El - £ 0 2 ) 2 '’"[1 + 

E 4 = (El - £ 0 2 ) 2/m [l + 


(C.7.b) 


(C.7.c) 


The exact expression for the effective cooling velocity term, [ 7 e 2 , in equations 
(C.4.a) and (C.4.b) can be expressed in terms of the mean velocity components 
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{U, V, W ) and the Reynolds stresses (v~v~, ij = 1,2,3) at the four wire positions 
shown in Fig. C.2. as follows (recalling that different expressions are required for 
Method A and Method B ): 


Method A 


(see Appendix D for derivation) 


{U 2 ) i + £] 

+ K3,ih 2 BN [(—) 2 + jp] + 2/^2, i[(—) + jp]} 

(U?h =Ko,2U 2 {K lt2 [l + ^) + [{jjf + 

+ ^3,2h 2 BN [( — ) 2 + jp] - 2K 2 ,7[(^)+ Jp]} 


(C.S.a) 


(C.&.b) 


2 

(^e 2 )3 = ATo,3^ 2 { 2(^.3 + h 2 BN I< 3,3)[1 + Jp] 

+ t(^) 2 + ^ 72 ] + 2^ 1,3 + + Jp ] 

W HIV 

— (•K’l.a — ^bjv^.s)^) + fjyll 


(C.8.c) 


(U 2 )i =Ko,4U 2 {-(Ki t 4 + hg N A' 3 i 4 )[l + y 2 ] 


, r /^2 , v2 i , W V (C.8.d) 

+ [(??) + 777 ] + o ( Al > 4 + h BN K 3,4)[( jj ) + rr 2 J 


•If 


t/ 2 


u 


u 2 


w. 


uw. 


+ (^1,4 — ^BN^3,4)[( y ) + ^ 2 ]} 


Method B (see Appendix A, reference [47] for derivation) 


{U 2 ) i =/v 0 , 1 L/ 2 {[l + ^] + -K’l.i K 77) 2 + 777] 


■If 1 


tf 2J 


+ ^3, 1 fcU(f) 2 + ^] + 2^,,[(^)+p 


(C.9.a) 
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(U c 2 ) 2 =Ko }2 U 2 {[\ + — ] + K lt2 [(—) 2 + jp] 


,W 


+ K 3i 2h 2 BN [( — ) 2 + Trj] “ 2K2,2[(jj) + 7 - 72 I) 


ttt; , 


U 


U 2 


u ’ U 21 


(C.9.6) 


(Ue 2 )z =Ko, 3U 2 {[1 + jp] + K 3j 3 h 2 BN [( — ) 2 + Jjj] 


w 


+ ^- 3[(^) 2 + ^]+ 2 ^ 2 , 3 [(^)+^]} 


W „ uw. 


u 


u 2 


U J U 2 


(C.9.c) 


(U e 2 ) 4 =/^ 0> 4^ 2 {[l + Jp] ■*" K 3 ,4h 2 B N[(jj) 2 + Jp] 


w. 


IV ‘ 


V, 

u 

. w , 


uu>. 


+ Ki A [{^j) 2 + jp\ - 2 K 2 ,4[{jj) + jp]} 


(C.9.d) 


where Ao, A’i, AT 2 and AT 3 (the second subscript denotes wire position) are func- 
tions of the wire inclination angle (a), the normal cooling coefficient (hfj) and 
the tangential cooling coefficient (k ) ; and h B fj is the binomial cooling coeffi- 
cient. The functioned forms for these variables are given in Appendix D. When 
equations (C.8.a) through (C.8.d) or equations (C.9.a) through (C.9.d) are sub- 
stituted into equations (C.4.a) and (C.4.b), the following non-linear algebraic 
equations, applicable to both Method A and Method B, result: 


Kq v + K U (V/U) + A' 2 „(V nvf + K, v (W/V) 2 +C, = 0 (C.lO.a) 
K 0w + K lv ,(\V/U) + KiJW/Uf + KiJV/Vf + C„ = 0 (C.10.6) 

where 

C„ = K„ v (^) + A',„(p) + K u (£) + (C.ll.a) 


7/2 7/7/1 

C w = K 0w (jp) + Kiv(jp) + K 2w {—) + K 3w (jp) (C.11.6) 

and where, for convenience, the following definitions have been employed: 
Method A 

Ko v =Ao ) iAT] i i — A 0l2 A li2 E 12 
Kiv =27<"o,iA' 2) i + 2Ao, 2 A' 2i2 Ei 2 
I\ 2v =A'o,i — Ao, 2 Ei 2 
A^t, =h 2 BN (K 0 ,\ 7v 3> i — A r o >2 A r 3|2 Ei 2 ) 


(C.12.a) 
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Kow =(-^)[Ko t 3{Ki t 3 + h 2 BN K3 t 3) 

— #0, 4(^1, 4 + h 2 BN K3 t 4^34] 

K lw = — [i^o,3(J^i,3 — h 2 BN K. 3 , 3 ) 

+ Ko,4(Kl,4 — h 2 BN K3 t 4'E34] 

1 ^ 2 ui —Kqw 
K 3 W = Kq,3 — -^0,4^34 

Method B 

Kov =I<o,i — Ko, 2^12 
K\ v =2 Kq,\ K 2t i 2Ko, 2 ^ 2 , 2^12 
K 2v =Ko,lKl,l — ^0,2^1,2^12 

&3V =h 2 BN (Ko t \K 3 t i — K\ 2 K3,2^2) 

and, 

Kq w =K 0t 3 — Ko,4^34 
Klw =2/vo,3-^2,3 + 2Ko > 4 K 2 , 4 S 34 
K 2w =A'o, 3 /vi, 3 — ATo ,4 K 1,4^4 
K 2w =h 2 BN (K 0 ,3K 3 ,3 - Kq t 4 K 3 f 4 S 34 ) 

The terms C v and C w in equations (C.lO.a) and (C.lO.b) account for the errors 
which would result from neglecting the fluctuating velocity components relative 
to mean components in the expression for the effective cooling velocity. 

In a low turbulence intensity field (< 10%), the correction terms C v and C w 
in equations (C.lO.a) and (C.lO.b) can be neglected without introducing serious 
errors in the computed mean velocity ratios ( V/U ) and ( W/U ). Under these 
conditions, the mean response equations reduce to: 

Ko, + KUV/U) + K»<y/ur + K^W/U ) 1 = 0 (C.14.0) 

Jvou> + K lw {W!U) + Kijw/vy + K 3w (V/U) 2 = 0 (C.14.6) 

Equations (C.14.a) and (C.14.b) represent two coupled non-linear algebraic equa- 
tions which can be solved simultaneously for the relative mean velocity compo- 
nents (V/U) and (W/U). When second-order mean velocity terms in equations 
(C.14.a) and (C.14.b) are neglected, solutions for (V/U) and (W/U) which are 
first-order accurate may be obtained directly as: 

(V/U ) 0 = —Kov/Kiv 

(W/U) 0 = -Kow/Ki w 


(C.12.6) 


(C.13.a) 


(C.13.a) 


(C.15.a) 

(C.15.6) 
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which axe appropriate for evaluating ( V/U ) and ( W/U ) in slightly skewed flows. 
For the more general case, however, equations (C.lO.a) and (C.lO.b) represent two 
non-linear algebraic equations in seven unknowns (V/U, W/U, u 2 /U 2 , v 2 /U 2 , 
w 2 /U 2 , uv/U 2 , uw/U 2 ). Therefore, five more independent equations are required 
for closure, which, in turn, requires rotation of the slant and/or normal-wire to 
at least five different positions. If this is done arbitrarily, then the coefficient ma- 
trix may become ill-conditioned. In order to overcome this difficulty, equations 
(C.14.a) and (C.14.b) are first solved for the mean component ratios (V/U) and 
(W/U) after initial estimates are made using equations (C.15.a) and (C.15.b). 
The calculated results are then used to solve the turbulence response equations 
(to be developed shortly) for the complete Reynolds stress tensor. The calculated 
stresses can then be used to correct the mean velocity components by simulta- 
neously solving equations (C.lO.a) and (C.lO.b). The improved mean velocity 
values axe subsequently used to update the Reynolds stress tensor. This iteration 
procedure is repeated until convergence is achieved. 

Once the mean velocity ratios (V/U) and (W/U) axe known, the axial mean 
velocity (U) can be determined by averaging its calculated value with the wire 
at each of the four positions shown in Fig. C.2. When the expression for the 
effective cooling velocity at Position 1 (equation (C.8.a) or (C.9.a)) is substituted 
into equation (C.2.a), the following equation for the axial velocity with the wire 
at Position 1 is obtained, after some rearrangement: 

Method A 


17, -{(£? - + CV,7(£? - £ 0 2 )]/(Bo)^") 1/! 


u 


V 


W, 


w 1 


{#1,1 [1 + jp \ + E( 7f) 2 + 777] + ^ 3 , ih BN [(—) + 777] 


U 


U 2 


u 


u 2 


+ 2 A - 2 ,.[(^)+^]}-‘ /2 


(C.16) 


Method B 


17 , ={(£? - E„ 2 ) 2/m [ 1 + C,e, 2 /(£? - B 2 )]/^) 2 '"} 1 ' 2 




. V , 


w. 


w i 


{[1 + Jp \ + # 1,1 [({/) 2 + \ + K 3 ,\ h 2 BN [( — ) 2 + 

+ 2 A - ! ,,[(£)+^!}-/ 2 


(C. 17) 


where B c , the slope of the calibration curve, is either: 


B C = B 


(C.lS.a) 
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when the calibration stream is normal to the wire, or: 

B c = B(K 0 ) m ' 2 (C.18.6) 

when the calibration stream is aligned with the probe body (as in Method B ), or: 

B c = B(K 0 K !) m/2 (C.18.c) 

when the calibration stream is aligned normal to the probe body (as in Method A). 
Similar expressions for the axial mean velocity can be derived for wire Positions 
2, 3 and 4. 

The total mean velocity (U 0 ) is related to the axial mean velocity (17) by: 

Uo = Uj cos 4> (C.19.a) 

whfTP 

<f> = tan“ 1 (V , r /17) (C.19.6) 

VrfU = [(V/U ) 2 + (W/U ) 2 ] i/2 (C.19.C) 

This section will be concluded with an explanation of how the possible probe 
support interference associated with Method A can be avoided if data are taken 
by means of probe rotation about the y-axis and an estimate of the transverse 
flow angle 7 in the x-z plane is known. 

Interference can be avoided if, rather than recording data relative to the 
primary coordinate system, data axe collected relative to a another coordinate 
system which is rotated about the y-axis through an estimated yaw flow angle 
70, which is close to the actual flow angle 7, such that the transverse velocity 
component in the x-z plane is small. This allows the use of a hot-wire probe 
body designed to tolerate small amounts of unknown yaw skewness. The mean 
velocity components tire computed by the method described above, relative to 
the secondary coordinate system, and then they are transformed to the primary 
coordinate system. Estimates of 70 can be established, at least to first-order 
accuracy, by pressure nulling of a two-tube Conrad probe, or by means of the 
hot-wire nulling technique described by Page [62]. 

C.4 Turbulence Equation Development 

As was the case for the development of the mean response equations, the 
starting point for the turbulence response equations is the instantaneous empir- 
ical cooling relation: 

(E 2 - El) 2/m = B 2/m U 2 


(C.l.c) 
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Also, as before, the left-hand side is expanded by means of the binomial theorem 
and truncated at second-order: 


(E 2 - = (£ 2 - 4 2 ) 2 '"[ 1 + WmJ-sA* + 


EP-E 2 E 2 — Eq 


(C.6.a) 


where, 

C = 2/m + (4/m)(£-l)(j 5 ^) (C.6.4) 

In contrast to the mean velocity development, the effective cooling velocity term 
in equation (C.l.c), is derived relative to another orthogonal coordinate system. 
For Method A, this coordinate system, x',y,z', is generated by a rotation about 
the y-axis such that the x'-axis coincides with the resultant mean yaw velocity 
vector, U r (refer to Fig. C.3.a). For Method B, a coordinate system, x,y',z', is 
generated by a rotation about the x-axis such that the z' = 0 plane coincides with 
the resultant mean transverse velocity vector, V r (refer to Fig. C.3.b). These new 
coordinates axe employed here because the Reynolds stress components, when 
calculated relative to these new coordinates, do not involve the solution of a 
full 6x6 system of equations under certain conditions. The exact forms of the 
instantaneous effective cooling velocity relative to these new coordinates are: 

(see Appendix E) 


Method A 


U t 2 = KoUSifm + /»(«' /U r ) + Mv/Ur) ~ f„{w' /U r ) 

+ K 12 (u ,2 /U r 2 ) + ( V 2 /U r 2 ) + K U (tV ,2 /U r 2 ) 
+ 2Kio(u'v/U r 2 ) - 2Ki 3 (u'w'/Ur 2 ) 

- 2K u (vw' /Ur 2 )} 


(C. 20) 


where 

f m = K u + 2K 10 (V/U r ) + (V/Ur) 2 
f u = 2K\2 + 2K 10 (V/Ur) 
f v = 2A'io + 2(V/U r ) 

f w =2K u + 2K u (V/Ur) 


(C. 21) 


Method B 


(see Appendix B, reference [45]) 


U 2 = IUU 2 {f m + f»(u/U) + f v (v'/U) + f w {w'!U) + ( u 2 /U 2 ) 

+ K 12 (v' 2 fU 2 ) + Ku(w' 2 /U 2 ) + 2K 10 (uv'/U 2 ) 
+ 2K u (uw' /U 2 ) + 2K 13 (v'w'/U 2 )} 


(C. 22) 
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where, 


fm = 1 + 2 K 10 (Vr/U) + K 12 {Vr/U ) 2 


/„ = 2 + 2K 10 (Vr/U) (C.23) 

f v = 2 K\q + 2Ki2(V r /U) 
f w = 2K n +2K u (V r fU) 

where K 10 , K u , K u , K 13 , K u are functions of the probe rotation angle, 0, and 
of Ki , K 2 and K 3 , which, in turn, are functions of h N , k and a (refer to Appendix 
D and E for explicit forms), and u\ v' and u;' are velocity fluctuations along the 
x\ y' and w' axes, respectively, as shown in Figs. C.3.a and C.3.b. If equation 
(C.6.a) and either equation (C.20) or (C.22) is substituted into equation (C.l.c) 
and subsequently time- averaged, then the resulting expression would not take 
into account first-order effects in the fluctuations and would, therefore, not be 
appropriate for analyzing low intensity flows. Also, the resulting set of equations 
would form a 6 x 6 system of equations whose solution can lead to inaccuracies in 
the calculated Reynolds stresses. In order to avoid these shortcomings, consider 
equation (C.l.c) which is first time-averaged to yield: 


(E 2 - E\ l) 2 !™ = B 7 f m U e 2 (C. 24) 

and then consider the difference between equations (C.l.c) and (C.24) in a nor- 
malized form which eliminates the contribution of the slope factor, B 2 m , namely. 


(E 2 - £ 0 2 ) 2/m - (E 2 - E 2 ) 2 /™ _ U 2 ~ Ue 2 
(£ 2 -£ 2 ) 2/m U e 2 


(C. 25) 


On the basis of equations (C.6.a) and (C.20) or (C.22), the time-averaged terms 
in equation (C.25) can be written as: 


(E 2 - El ) 2 = (E 2 


^> 2/ "[ 1 + i£W 


(C.26) 


Method A 

U? = KoUfifm + K 12 (^/Ur 2 ) + &/Ur 2 ) + K U (^ /U r 2 ) 

+2K,o(^/Ur 2 ) - 2K u (u'w l /Ur 2 ) ( C - 27 ) 

- 2K\\{yw' !U r 2 )} 


Method B 

Uf = K 0 U 2 {fm + &/U 2 ) + Kni^/U 2 ) + K 14 (w' 2 /U 2 ) 

+2 K 10 (^/U 2 ) + 2Kii(uui' /U 2 ) + 2K i3 (v r w' /U 2 )} 


(C. 28) 
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If the appropriate terms developed thus fax are substituted into equation (C.25) 
an expression which has the following functional form results: 

y(e,e , e^, Ui, tijUj, Ujtij) = 0, i,j = 1,2,3 (C.29) 

since all mean flow field variables are known from the previous mean flow analysis. 
Equation (C.29) implies that in order to retain first-order fluctuation effects, 
equation (C.25) should be squared and time-averaged: 


[(E 2 - Eg) 2 /™ - (E 2 - E 0 2 ) 2 /™] 2 _ \U 2 - U e 2 ] 2 
[(E 2 - Ej,) 2 / m ] 2 ” [U?] 2 


(C.30) 


Now, if the appropriate terms are substituted into equation (C.30), and the 
squaring and time-averaging performed, then, retaining only second-order terms, 
the following generalized turbulence response equations result: 

(see Appendix E) 


Method A 


F uu (u' 2 /U 2 ) + F vv (v 2 fU r 2 ) + F ww (w' 2 /U r 2 ) 

+ F uv (^/Ur 2 ) - F uw (^/Ur 2 ) - F vw (vw r /U r 2 ) = F mm 


(C- 31) 


where, 


F uu = A 2 — 2/ m nE 12 
F vv = f 2 - 2/ m n 


F W w = fw* — 2/ m nE 14 

F uv = 2 (/./. - 2/ m nE 10 ) (C.32) 

F uw = 2(/ u /u, — 2/ m nA' 13 ) 

E vl „ = 2(f v f w — 2f m UKn) 

F m m = /m 2 n 


Method B\ ( 


see 


Appendix B, reference [45]) 


F U u(u 2 /U 2 ) + F vv (v' 2 /U 2 ) + F ww (w' 2 /U 2 ) 

+ F uv (uv'/U 2 ) + F uw (uw' /U 2 ) + F vw (v'w'/U 2 ) = E mm 


where, 



= A 2 - 2/ m n 

= /t, 2 — 2/ m IIA , 'i2 
= U 2 -2f m UK 14 
= 2 (f u fv - 2f m UK 10 ) 

= 2(/ u /u, — 2/ m nA n ) 
= 2(f v f w - 2/ TO ILA' 13 ) 
= /m 2 n 


(C.33) 


(C.34) 
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with (applicable to both Method A and Method B ), 


II = (16/m 2 ) 


\e*£ 7 /(E 2 ~ EZf\ 

[l + 2Ce 2 /(F 2 -F 2 ) 2 ] 


(C.35) 


and fjn, fu, fv and f w are defined by equations (C.21) or (C.23). The compo- 
nent e 2 in the numerator of equation (C.35) is a consequence of the first-order 
fluctuation (e) in equation (C.6.a), while e 2 in the denominator is a consequence 
of the second-order fluctuation (e 2 ) in this equation. 


C.4.1 Low Turbulence Intensity Flow 


For flows of low turbulence intensity, the second-order terms in equations 
(C.6.a) and (C.20) or (C.22) can be neglected in comparison to the first-order 
terms. Substitution of the first- order forms of equations (C.6.a) and (C.20) or 
(C.22) into equation (C.30) yields the following turbulence response equations 
applicable to low intensity flows: 


Method A 


(see Appendix E) 


fu\W 2 /U r 2 ) + f v \v 2 /V 2 ) + f w 2 (w' 2 /U r 2 ) 

+ 2f u fv(yv/U r 7 ) - 2fJ w (M/Ur 2 ) - 2f v f w (vw' fU r 2 ) = 


(CM) 


(see Appendix B, reference [45]) 


fj(uyu 2 ) + f„\v’ 2 /v 2 ) + fj(w 2 /u 2 ) (C37) 

+ 2f u f v (uv i /U 2 ) + 2f u f w (uw'/U 2 ) + 2f v f w (v'w'/U 2 ) = it f m 2 

with (applicable to both Method A and Method B ), 

tt = (16 /m 2 )[SE 2 /(E 2 - E 2 ) 2 } (C. 38) 

In principle, the six components of the Reynolds stress tensor relative to the 
x',y,z' coordinates shown in Fig. C.3.a or to the x,y',z' coordinates shown in 
Fig. C.3.b could be determined by applying data taken at six fixed probe po- 
sitions to equation (C.36) or (C.37), respectively, for low intensity flows. This 
is possible, since, in general, the wire will be sensitive to all three fluctuating 
velocity components at a given rotation angle. The resulting 6x6 system of 
equations may, however, become ill-conditioned if these fixed probe positions are 
specified in an arbitrary manner. 

To increase the accuracy of his calculations, Al-Beirutty employed a sensi- 
tivity minimization criteria to Method B which significantly improves the condi- 
tion number of the coefficient matrix. He accomplished this by observing that 


Method B 
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at certain rotation angles (“optimum planes”), the wire is either insensitive to 
a particular transverse velocity fluctuation, v' or w\ or is equally sensitive to 
both. As will be shown shortly, when data are taken at these optimum planes, 
the resulting 6x6 system of equations becomes partially decoupled, and, as a 
result, the problem is reduced to the solution of several lower-order systems of 
equations. The behavior of the sensitivity coefficients in equation (C.30) will be 
examined to determine if similar criteria can be applied to Method A as well. 

The sensitivity coefficients /„, f v and f w appearing in equations (C.36) 
and (C.37) can be written in terms of the probe rotation angle 0 when the 
appropriate expressions for the K\ o through Ku coefficients (given by equation 
(E.7) in Appendix E) are substituted into equations (C.21) and (C.23): 

[Method ~A\ 


f u = 2{K\ cos 2 9 + h 2 BN K 3 sin 2 0) + 2K 2 cos 0( V/U r ) (C.39.a) 

/„ = 2 K 2 cos 0 + 2(V/Ur) (C.39.6) 

f w = 2(K\ — h 2 BN K 3 ) sin0cos0 + 2K 2 sm0{V jU r ) (C.39.c) 


Method B 

f u = 2 + 2 K 2 cos 0(V r /U) (CAO.a) 

f v = 2I\ 2 cos 0 + 2 (Ki cos 2 0 + h 2 BN K 3 sin 2 0)(V r /U) (C.40.6) 

f w = 2K 2 sin 0 *f- 2(A i — h 2 BN K 3 )sm0cos0(V r /U) (C.40.c) 


In order to examine the wire’s sensitivity to the fluctuating velocity components, 
it is expedient to plot the Reynolds stress coefficients appearing in equations 
(C.36) and (C.37) as a function of probe rotation angle ( 0 ) for various levels 
of skewness. The normal and slant-wire cases for Method A are shown in Figs. 
C.4 and C.5, respectively, and for Method B in Figs. C.6 and C.7, respectively, 
utilizing typical values for the cooling coefficients: k = 0.15, hw = 1.0 and hsN 
= 0.9. Note that for Method B, only positive skewness angles need be considered 
since negative singles can be transformed to positive angles by a 180° rotation of 
the probe. 

Method A 

In order for the hot-wire to be insensitive to the w' fluctuating velocity 
component, the coefficient f w must be equal to zero. An examination of Fig. C.5 
reveals that f w is zero for a slant-wire at any skewness level when the wire is in 
the 0 and 180° positions. These locations will be denoted as Positions A and B, 
respectively. The offset prong slant-wire configuration (refer to Fig. C.l.b) was 
designed to operate in these positions with minimal probe support interference. 
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For a normal-wire (see Fig. C.4), f w is zero at 8 = 0° and 270°, although the 
former position is not admissible due to interference effects. The latter position 
is designated as position N a ■ The above conditions allow three of t he R eynolds 
stresses to be calculated independent of the remaining three. The u' 2 , v 2 and 
u'v stress components can be obtained by applying normal-wire data at position 
N a and slant-wire data at Positions A and B to equation (C.36) and solving the 
resulting 3x3 system of equations. 

With three stress components known, attention is turned to the condition 
when the wire is insensitive to the v fluctuating component, which occurs when 
f v is zero. For Method A, this condition proves of little use, inasmuch as for the 
limiting case when the skewness is zero, the slant-wire positions where /„ = 0 
correspond to positions where f w = 0, namely; 6 = 90 and 270°. Also, for 
the normal-wire, in general, there is no location where /„ = 0. The condition 
where the slant-wire is equally sensitive to both transverse fluctuating velocity 
components is also of little use. The reason for this is because the angular 
locations where these conditions occur do not always coincide with, or are even 
near, the maximum sensitivities of the unknown stress components. Therefore, 
for Method A , it appears that fixed probe positions are necessary. 

With reference to Fig. C.5.b, there are four fixed positions, 6 = 45, 135,225 
and 315°, designated as Positions E a ,F a ,G a and H A , where the slant- wire is 
always strongly sensitive to the vw' stress component, at least over the skewness 
range considered. The straight prong slant-wire configuration (refer to Fig. C.l.c) 
was designed to operate in each of these positions with minimal probe support 
interference. Slant- wire data taken at these four positions can be used to evaluate 
the vw' stress component directly. An algebraic expression (derived in Appendix 
E) for this component is given as: 

vw' —{(fufw)F A ( 7t E A ~^H a )Ujt^)e a ~ (fufw)E A ( V FA ~ ) (/m )ful 

V? = (2f u f w ) FA (2f v f w ) EA - (2fuU)E A (2fvfu,)F A 

(G.41 j 

At this point, expressions for four of the six Reynolds stresses have been devel- 
oped and there are no more wire positions available that will further decouple 
the response equations. Therefore, in order to get the remaining two Reynolds 
stresses {w ' 2 and u'w'), at least two additional probe positions are required, and, 
ideally, these positions should be at rotation angles where the sensitivity to the 
unknown stresses is at or near a maximum. With reference to Fig. C.4.b, for the 
normal-wire, two fixed positions exist where the wire is nearly equally sensitive 
to both u' and w' fluctuations, namely; 6 = 225 and 315°, which are denoted as 
Positions C A and D A , respectively. When normal- wire data at Positions C A and 
D a and the four known Reynolds stresses are applied to equation (C.36), then 

a 2 x 2 system of equations results, which, when solved, yields the w ' 2 and u'w' 
Reynolds stresses. 
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Method B 

With reference to Fig. C.7, two slant-wire probe positions result when the 
condition f w = 0 is applied to Method B: 9 = 0 and 180°, which are denoted as 
Positions A and B , respectively. The normal- wire also has two positions where 
f w = 0: 9 = 0 and 90° (refer to Fig. C.5). For the normal- wire, 0 = 0 has the 
additional benefit that /„ = 0 as well. This position is denoted as position Ng. 
When normal-wire data at position Ng are applied to equati on (C.37) a single 
algebraic equation for the u 2 stress component results. The v' 2 and uv' stress 
components can be obtained by applying slant-wire data at Positions A and B 
to equation (C.37) and solving the resulting 2x2 matrix. 

When the condition /„ = 0 is applied to Method B, two discrete slant- 
wire positions result, designated as Cg and Dg in Fig. C.7.a. At these positions, 
which are not fixed, but rather vary with skewness level, the slant- wire is sensitive 
only to u and w' fluctuations. If the u 2 stress component, which is known from 
previous calculations, and slant-wire data at Positions Cg and Dg are app lied 
to equation (C.37), then a 2 x 2 matrix results which may be solved for the w ' 2 
and uw' stress components. The remaining unknown stress component, v'w 1 , is 
determined by solving an algebraic expression which utilizes slant-wire data at 
Positions Eg, Fg, Gg and Hg shown in Fig. C.7.b. These four wire positions 
correspond to locations where the wire is equally sensitive to both transverse 
fluctuating velocity components, or, in other words, where f 2 = f w 2 . The 
expression for v'w' (derived in Appendix C of reference [45]) is given as: 

vW_ _ (/u/uQfb^Eb - 7r//fl)(/m 2 )E B ~ (/u/u>)e b (^Fb ~ '^Gb )]JvS)Fb 

U> ~ (2/ u /u,)fb(2/u/u))fb ~ (2/u/u))eb(2/ v /u>)f b 

(C.42) 

For all slant-wire measurements referred to Method B , the straight prong probe 
configuration shown in Fig. C.l.c is adequate, since probe interference effects do 
not occur with the wire in Positions A through Ng , provided that the skewness 
level does not exceed 30° in a low intensity flow and 20° in a moderate intensity 
flow [45]. 

C.4.2 Moderate Turbulence Intensity Flow 

The preceding analysis is applicable only to low turbulence intensity flows. 
For the more general case, equation (C.31) or (C.32), as applicable, must be 
solved. Second-order velocity fluctuations which appear in the expressions for 
the effective cooling velocity (equation (C.20) or (C.21)) prevent the decoupling 
of the stress components, so that a full 6x6 system of linear equations must be 
solved. 
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Method A 

For Method A, slant-wire data at Positions A,B,Ca and Da, in addition 
to normal- wire data at position N a, result in a system of five coupled linear 
equations. The sixth equation is obtained by utilizing slant- wire data obtained 
at Positions E a , Fa, G a and Ha- The generalized response equation (derived in 
Appendix E) is given as: 


(SiK 12 , Ea “ S 2 K UiFA )(u' 2 /U r 2 ) + (Sl - S 2 )(vyU r 2 )+ 

(siKu.Ea - S 2 Ki4,F A ){w ' 2 /Ur 2 ) + 2(si ATiO, Ea - S2Kio,F A )( u ' v /Ur 2 )- 


(S lvw - S 2 vw ){vw' /Ur 2 ) = -{sJ^Ea - s2 fm,F A )/ 2 


(C.33) 


where 


51 = — 8 [(/u/u>)Fa — (fmHi 3 )F A {^FA + He* )](n£7 B H H B )(fm)E A 

5 2 = —S[(fufw)E A ~ (fmKi3)E A (J^EA d^Fa^C^F* — n<?y* )(/»» i)Fa 


Sivw = l6[(fufw)F A ~ A'i3,F^(/m)F^(nF >l + n G ,«)]x 
[{fvfw)E A ~ Fu,E A (fm)EA(^EA + n H-»)l 
S 2 vw = lG{(fufu,)E A - K 13 ,EA(fm)EA(nE A +KHa)]x 
[( fvfw)F A ~ Fu,F a (fm)F A (^Fa + n Gx)l 

and where II is given by equation (C.35) and K\o through A 14 are given by 
equation (E.7). 

Method B 

For Method B, slant-wire data at Positions A,B,C b and D B , in addition 
to normal- wire data at position N B , result in a system of five coupled linear 
equations. The sixth equation is obtained by utilizing slant-wire data obtained 
at Positions E B , F B , G B and H B . The generalized response equation (derived in 
Appendix C of reference [45]) is given as: 

(s 1 - s 2 )(u 7 /U 2 ) + (siK 12iEb - s 2 Ki 2 t F B )(v' 2 /U 2 )+ 

(siA'i4 i f b — s 2 Ku t F B ){ W ' 2 /U 2 ) + 2(si A io.Eb — s 2 K\o,f b )( uw ' /U )+ (^-^) 
('S’lwui — S 2 vw )(v'w'/U 2 ) = —{sifm,E B ~ s 2 fm,F B )/ 2 


where 

Sj = -8 [(/u/uOFb - A' ni FB(/m)FB( n FB +n G fl)](n Efl - U HB )(fm)E B 
s 2 = -8[(/ u /u,)fb - K u ,E B (fm)E B ( n^B +n// B )](n F B -n G B)(/-)FB 
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and, 

S\vw ~ 1 6[(/«/u>)f b - Kn,F B (fm)F B (^FB +n GB )]x 

[(/i/ui)Eb ^13, Eb (.fm )f B (^Eb d" IIi/ B )] 

$2vw = 16[(/u/ur)£; B d" H// B )] X 

[{fvfw^Fa -^13, F® (/n»)F B (Hf b d" n Ga )] 

and where II is given by equation (C.35) and K 10 through Ku are given by 
equation (E.7). 

C.5 Solution of Mean and Turbulence Equations 
C.5.1 Low Turbulence Intensity Flow 

For low intensity flows the velocity ratios V/U and W/U and, in turn, V/U r 
{Method A ) or V r /U ( Method B ), are found by applying the Newton- Raphson 
procedure to the simultaneous equations (C.14.a) and (C.14.b). First-order ap- 
proximations given by equations (C.15.a) and (C.15.b) are used as initial esti- 
mates in the calculation procedure. Once V/U r or V r /U (as applicable) is known, 
the optimum planes locations are determined from equations (C.39) or (C.40), 
as applicable. Data can then be taken at these locations or interpolated to these 
positions from data taken in fixed increments. These data are applied to equa- 
tions (C.36) and (C.41), or (C.37) and (C.42) from which the Reynolds stress 
components are computed relative to the secondary turbulence coordinate sys- 
tem. The Reynolds stress tensor is then transformed to the primary coordinate 
system. 

C.5. 2 Moderate Turbulence Intensity Flow 

For moderate intensity flows, the velocity ratios V/U and W/U axe com- 
puted by simultaneously solving equations (C.lO.a) and (C.lO.b). These equa- 
tions, however, require knowledge of the Reynolds stress components. To allevi- 
ate this problem, V/U and W/U are first estimated by solving equations (C.14.a) 
and (C.14.b). From these results, optimum planes are established and Reynolds 
stresses relative to the secondary coordinate system are computed using either 
equations (C.31) and (C.43), or (C.33) and (C.44), as applicable. The stress 
tensor is then transformed to the primary coordinates and substituted into equa- 
tions (C.lO.a) and (C.lO.b) from which improved estimates of V/U and W/U are 
calculated. These improved estimates are, in turn, used to update the Reynolds 
stresses and so on until convergence is achieved. 

A disadvantage of the above procedure is that two slant-wire probe configu- 
rations are required for Method A in order to obtain data at Positions A, B and 
Ea through Ha ■ An alternative procedure would be to obtain slant-wire data at 
three (or more) fixed positions; 6 = 45,225 and 315°, all of which are attainable 
with the offset prong slant-wire probe configuration shown in Fig. C.l.b. By 
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applying slant-wire data at these positions and norm al-w ire data at Positions 
Ca and Da to equation (C.36) the stress components, uj' 2 , u*w 9 and vw* can be 
solved so as to minimize the error (in a least-squares sense) associated with the 
over-specified conditions. It should also be noted that in many practical flows, 
the secondary shear stress vw is at least an order of magnitude less than the 
other stress components and can be neglected. Under these circumstances there 
is no need for the straight prong slant-wire probe shown in Fig. C.l.c. 
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a) Normal-wire. 



b) Offset prong slant-wire. 



c) Straight prong slant-wire. 

Fig. C.l. Hot-wire probe configurations for Method A. 
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WIRE IN THE x-y PLANE 
(LOCATED AT O) 



WIRE IN THE x-z PLANE 
(LOCATED AT O ) 


Fig. C.2. Hot-wire positions for mean flow determination. 
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a) Method A. 



ROTATABLE HOT 
WIRE PROBE 



b) Method B. 

Fig. C.3. Reference coordinates for Reynolds stress determination. 


(-2f u f w ) (-2f v f.) 
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sjU«!3!i^aoQ ssajjs spjouAay 


a) f^, fj and 2f u f v coefficients. b) f£, -2f u f w and -2f v f w coefficients. 

Fig. C.4. Behavior of Reynolds stress coefficients, normal-wire. Method A. 
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o 

(N 


sjuopuj»oo 8S»J*S spiou/ey 


a) f 2 ,f v 2 and 2f u f v coefficients. b) f£, -2f u f w and -2f v f w coefficients. 

Fig. C.5. Behavior of Reynolds stress coefficients, slant-wire, Method A. 
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s|uei 3 [jj »03 ss»J|s sp|ou/*y 


a) f„,Cj and 2f u f v coefficients. b) f^,2f u f w and 2f v f w coefficients. 

Fig. C.6. Behavior of Reynolds stress coefficients, normal-wire. Method B. 






170 














APPENDIX D 


SUPPLEMENTARY ON 
MEAN FLOW EQUATION DEVELOPMENT 


This appendix contains supplementary information on the derivation of the 
mean flow response equations for Method A. The corresponding information for 
Method B can be found in Appendix A of reference [45]. 


D.l Effective Cooling Velocity 


In order to develop an expression for the effective cooling velocity, the in- 
stantaneous total velocity vector must be decomposed into components which 
act normal and tangential to the wire. Fig. D.l shows these decomposed com- 
ponents acting on a wire positioned in the x"-y plane. In Fig. D.l, the U\, U 2 
and U 3 velocity components are directed along the x ,y ,2 coordinate axes, re- 
spectively, which rotate with the wire about the y-axis and are related to the 
laboratory (primary) x, y , z coordinate system through the rotation angle A. The 
normal, parallel and binormal effective cooling velocity components acting on a 
wire which is rotated through an arbitrary angle A about the y-axis relative to 
the x,y,z primary coordinate system are: 


(U e )N = h^j{[(U + u) cos A — (W + w ) sin A] sin a + (V + v) cos a} (D.l) 
(Ue)p = k{[(U + u) cos A — (W + ttf) sin A] cos a — (V + t>)sina} (D.2) 
( U c )bn = hBN{(U + u)sin A + (W + u>)cosA) (£-3) 


where hjv, h, and h q y represent weighting factors, in the normal, parallel and 
binormal directions, respectively. The total effective cooling velocity is the Eu- 
clidian sum of the individual components: 


Ue 2 = ( U e f N + ( U t fp + {U e f BN (DA) 


Substitution of equations (D.l), (D.2) and (D.3) into (D.4) yields the exact ex- 
pression for the effective cooling velocity for a hot-wire rotated an arbitrary angle 
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A about the y- axis: 


U e 2 = KqU 2 {K\i[\ + 2 (— ) + (^ 7 )] 


V 2 v V v 2 
+ [(^7) + 2 <^%) + (^)] 

uj 

+ Ah[(^-) + 2(^)(^) + (^)] 

+ 2iC 1 o[(^) + (jj) + (|y)(^) + (^)] 

-2 *,,[(£) + (E) + (£)(£) + ( -)] 

- 2 *„|<£)(£) + (£)(£) + (£)(£) + Q] 


(0.5) 


where, 


and, 


A”io — AT 2 cos A 
An = AT 2 sin A 

K \2 = A] cos 2 A h 2 BN Kz sin 2 A 
A'13 = (A'i — sin A cos A 

A14 = I\ 1 sin 2 A *f- figfyKz cos 2 A 


K 0 = (h 2 jv cos 2 a + k 2 sin 2 a) 

K\ = ( h 2 N sin 2 a + A: 2 cos 2 a)/A 0 
J\ 2 = ( h 2 N — A: 2 ) sin a cos a/ Ao 
A 3 = 1/A'o 

Time-averaging equation (D. 5 ) yields: 


(DA) 


(D. 7) 


^ 2 = A ' oC / 2 { A 12 [1 + ^] + [(^) 2 + ^] 


t/ 2 ' 


C/' 


t / 2 


,w. 


+ Ai 4 [(^) 2 + ^] + 2 Aio [(^)+^] 

- 2 if„K 2 ) + p ] _a riI [(^) + P)} 


l/V . 


(D.8) 


The expressions for the effective cooling velocity at wire positions 1 , 2,3 and 4 
in Fig. 4.1 are obtained by substituting the corresponding rotation angles (A = 
0 , 180 , 45 and 315 °) into equation (D.8). At wire positions 1 and 2, where a 
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slant- wire (a = 45°) is utilized, the following expressions result: 


(U t 2 )i =K 0 ,iU 2 {K hl [l + + [(^) 2 + jp\ 


+ Kz,\h 2 BN [{ — ) 2 + y 2 ] + 2K2,l[(jj) + y 2 ]} 


P-9) 


(v,*h + £l + l(^) 2 + £] 10) 

+ K„h\ N \(^f + j£] - 2K,^) + j£]} 

A normal- wire (o = 0°) is utilized at positions 3 and 4, for which the coefficient 
K 2 = 0. The expressions for the effective cooling velocity at these positions are: 


(U 2 ) 3 =Kq^U 2 {-{K \ i 2 + h 2 BN K3j){l + y 2 ] 

+ [( Jjf + ^ 2 ] + + h BN K 3,3)[(jj) 2 + Jp] 

W uw 

— Pi, 3 — I) 


P-11) 



D.2 Cooling Coefficient Modeling 

In general, the cooling coefficients h/v, k and h B ^ are functions of the wire 
rotation angle. The tangential cooling coefficient, k, is modeled by a correla- 
tion developed by Friehe and Schwarz [63]. This correlation accounts for the 
dependency of k on the wire orientation angle (ip): 

k(xp) = {[1 - i> r (l - cos 1/2 ip)} 4 — cos 2 V’} 1/2 / sin (D.13) 


1 — [fc r 2 sin 2 a -I- cos 2 a] 1 / 4 
(1 — cos 1 / 2 a) 


where, 
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and where xp is the angle of incidence between the total velocity vector and the 
normal to the axis of the wire. The functional dependency of xp on probe rotation 
angle for Method A is given by: 


The reference value of k r is determined from a calibration procedure in fully- 
developed pipe flow. This procedure is described in detail by Al-Beirutty in 
Appendix E of reference [45]. 

In order to demonstrate the behavior of k, it is expedient to examine its 
behavior relative to the x ' ,y,z' coordinate system shown in Fig. C.3.a. In this 
coordinate system, the mean transverse velocity component directed along the 
w'~ axis has been nulled out. The probe rotation angle for the x' , y, z' coordinate 
system is specified as 0 and is related to the rotation angle in the x, y, z coordinate 
system by: 

0 = A + tan - X (W/U) (D. 14) 

Relative to this new coordinate system, the dependence of the wire orientation 
angle {xp) on probe rotation angle ( 0 ) is given by: 

ip = 7t/ 2 — cos -1 [cos £ cos a cos 0 — sin £ sin a] 
where £ is the pitch flow angle defined as: 

£ = tan -\VlU r ) 

A plot of k versus probe rotation angle 0 for various values of £ is shown in Fig. 
D.2.a for a normal-wire, and in Fig. D.2.b for a slant-wire. 

For Method B, Arterberry [47] determined th at the normal and binormal 
cooling coefficients could be specified as constants; (h# = 1.0, Hen — 0-9)- This 
is a reasonable assumption since for low-to-moderate skew angles, the normal 
velocity component is much larger then the binormal component, regardless of 
probe rotation angle. For Method A, however, equations (D.l) and (D.3) show 
that the normal and binormal velocity components for the slant-wire are strongly 
dependent on probe rotation angle, even under unskewed conditions. For exam- 
ple, consider a slant-wire in unskewed flow conditions (V=W=0) in the absence 
of turbulence fluctuations. In the limit as A approaches 0 or 180° (refer to Fig. 
D.l and equations (D.l) and (D.3)), the effective normal velocity component 
(U e )w approaches hfjU sin a and the effective binormal component is zero. For 
this case, 1.0 and hBN=0.9 , since these positions axe equivalent to positions 
associated with Method B. In the limit as A approaches ±90°, the normal veloc- 
ity component approaches zero and the effective binormal component approaches 


xp = tt/ 2 — cos" 


_i [ cos a cos A — ( V/U ) sin a — ( W/U ) cos a sin A 

[ (i -i- ( VJuy + Jw/vyyf* 
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the total velocity i.e., (U e ) B N=h B fjU=Uo- For this case, h^— 0.9 and h B s 1*0* 
These limiting cases for the slant-wire were accounted for in the present study by 
specifying h N = 1.0 and h BN =0.9 when \U N \ > \Ubn\, and h N = 0.9 and h B N = 10 
when |17 bn| > \Un\- For the normal-wire, over the useful range of rotation, and 
for zero-to-moderate skewness levels, the binomial velocity component is always 
significantly greater than the normal component, so that the h N and h BN can 
be specified as 0.9 and 1.0, respectively. 
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APPENDIX E 


SUPPLEMENTARY ON 
TURBULENCE EQUATION DEVELOPMENT 

This appendix contains supplementary information on the derivation of the 
turbulence response equations for Method A. The corresponding information for 
Method B can be found in Appendices B and C of reference [45]. 

E.l Effective Cooling Velocity 

In order to develop an expression for the effective cooling velocity relative 
to the secondary turbulence coordinate system, the instantaneous total velocity 
vector must be decomposed into components which act normal and tangential 
to the wire. Fig. E.l shows these decomposed components acting on a wire po- 
sitioned in the x"-y plane. In Fig. E.l, the U u U 2 and U 3 velocity components 
are directed along the x",y,z" coordinate axes, respectively, which rotate with 
the wire about the y-axis and are related to the turbulence (secondary) x,y,z 
coordinate system through the rotation angle 0. The normal, parallel and bi- 
normal effective cooling velocity components acting on a wire which is rotated 
through an arbitrary angle 0 about the y-axis relative to the x',y,z' secondary 
coordinate system are: 

(U e ) N = h N {[(V r + u')cos8 - w' sm6]sina + (V + v)cosa) (E.l) 
(U e )p = k{[(Ur + u')cos0 - w' sin0]cosa - (V + r)sino) (E. 2) 
(U e ) B N = h B N{(U r + u')sin0 + u/cos0} (£-3) 

where h N , k, and h B x represent weighting factors in the normal, parallel and 
binormal directions, respectively. The total effective cooling velocity is the Eu- 
clidian sum of the individual components: 

U t 2 =(U e ) 2 N + (U e ) 2 P + (U e ) 2 BN (EA) 

Substitution of equations (E.l), (E.2) and (E.3) into (E.4) yields the exact ex- 
pression for the effective cooling velocity for a hot-wire rotated an arbitrary angle 
0 about the y-axis: 

U 2 = KoUSifm + fu(u'/Ur) + fv(v/Ur) - f w (™'/Ur) 

+ Ku(u' 2 /U r 2 ) + ( v 2 /U r 2 ) + Kn(w ,2 /U r 2 ) 5) 

+ 2/vio(u 'u/t/r 2 ) - 2K\z(u'w' jU 2 ) 

- 2K u (vw' /U r 2 )} 


pAot_12i- mTtNT!0NftLLY 


FRECEDING PAGE BLANK NOT FILMED 
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where, 


fm = K l2 + 2K l0 (V/U r ) + (V/U r ) 2 
f u = 2K u + 2K 10 (V/U r ) 
f v = 2K 10 + 2(V/U r ) 
f w = 2K u + 2K u (V/U r ) 


(E.6) 


and, 


ivj 0 = K 2 cos 9 
K\\ — K 2 sin# 

K 12 = K 1 cos 2 6 + h 2 BN K 3 sm 2 6 (£.7) 

K\ 3 — ( K\ — h 2 BN K 3 )smBcos9 
Ku ss ATi sin 2 9 + h 2 BN K 3 cos 2 9 

and where the coefficients Ko through K 3 are given by equation (D.7). 

E.2 Turbulence Response Equation 

The working expression for the mean-square voltage response of a hot-wire 
transducer in terms of the effective cooling velocity is given by equation (C.30) 
as: 


j(£ 2 - E 2 ) 2 /™ - (£ 2 - £ 2 ) 2 /™] 2 [f/ e 2 - t/e 2 ] 2 


[(£ 2 - Elf !'"] 1 


\V,‘\ 


*12 


(£. 8 ) 


Substitution of equation (E.5) and its time-averaged form into the RHS numer- 
ator of equation (E.8) yields: 


[1 U e 2 - U 2 ) 2 =K 0 2 U r 4 {fJ(u' 2 /U r 2 ) + fA^/Ur 2 ) 

+ fJi^/Ur 2 ) + 2f u f v (t?v/U r 2 ) 

- 2f u f w (y^/U r 2 ) - 2 f v f w (v^/Ur 2 )} 
+ O[2f u K 12 K 0 2 U r 4 (^/U r 3 )\ 
and the RHS denominator of equation (E.8) becomes: 


(£.9) 


[U e 2 ] 2 =K 0 2 U r 4 {f m 2 + 2/ rn [A' 12 (u' 2 /£ r 2 ) + ( v 2 /U r 2 ) 

+ Kui^/Ur 2 ) + 2K 10 ffi/U r 2 ) - 2K 13 (M/U r 7 ) 


(£. 10 ) 


- 2K n {vw' /U 2 )]} +0[KoK l2 U 2 {u' 2 fU 2 )] 2 

where the last terms on the RHS of equations (E.9) and (E.10) represent the 
order of the truncated terms. The first term in the LHS numerator of equation 
(E.8) is given by equation (C.6) as: 

(E 2 - Elf' m = (E 2 - Eifl-l 1 + + fTTEjl 
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where. 


2 

C = 2/m + (4/m)(- - 1)( & _ - g ?) 


(£. 12 ) 


Substitution of equation (E.ll) and its time- averaged form into the LHS numer- 
ator of equation (E.8) yields: 

e*E 2 


[(E 2 - E 2 ) 2 / m - (E 2 - E 2 ) 2 /™) 2 = (E 2 - £o) 4/n, [(4/m) 2 ^2— - 2 ^] (£.13) 

and the LHS denominator of equation (E.8) becomes. 

o Ce? 

[(£2 _ Ej)»/-f = (£ J - £ 0 2 ) 1/m [l + (jnrgjjl (E ' 14) 

where third and higher-order voltage terms have been neglected in the develop- 
ment. 

Substituting equations (E.9), (E.10), (E.13) and (E.14) into equation (E.8), 
and retaining only second-order terms, the following generalized turbulence re- 
sponse equation results: 


where, 


F uu (u' 2 /U r 2 ) + F vv (v 2 /U r 2 ) + F ww (w' 2 /Ur 2 ) 

+ FuvilFv/Ur 2 ) ~ F uw (u'w' /U 2 ) - F vw (vw' /Ur 2 ) = Fmi 

F uu = f 2 - 2/ m nA r 12 
f vv = f 2 - 2/ m n 
F ww = fj - 2/ m nA 14 
F uv = 2 (/./„ - 2/ m nA 10 ) 

F uw = 2 (/„/«, - 2/ m nA 13 ) 

Avu, = 2(f v f w -2f m UK u ) 

F m m ~ fm II 


(£.15) 


(£.16) 


and, 


n = (16/m 2 ) 




(£.17) 


(1 + 2 CeV(£2 - £ 0 2 ) 2 ) 

and f m , /„, /„ and / u , are defined by equation (E.6), and where K 10 through 
A'h are given by equation (E.7). 
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E.3 Algebraic Expression for vw' l Ur 2 

When slant-wire data taken at positions Ea, Fa,Ga and Ha are applied to 
equation (E.15), the following four equations result: 

(F.uhA^/vS) + (f vv )e a <?/uS) 

+ (.f^eA^/vS) + (F.AeJ'Jv/US) ( e - 18 ) 

- (F uw ) Ea (u’w' /UA) - (F vw )E A {yw’ /UA) = (F mm )E A 

(F..)fA^/UA) + (ttW A) 

+ (^(^’1 + (F..) f ,(S/t;, ! ) ( £19 ) 

(-Fuu;)F y4 {u’lV 9 /U r ^) “ (F*vw)f a (viV* fU r ) = (-Fmm)fA 

(F.uhA^/UA) + (F ve )c A (S/UA) 

' +(F»»)G^^/^ S ) + (F.„) C .(^/^ 2 ) (£- 20 > 

- (Fu^gJu'w'/UA) - (F, w )aA vw '/Fr 2 ) - (Fm m )a A 

(F.Ah^/VA) + (f,.KW ! ) 

+ (FM^JVA) + (F..)h a (?v/UA) (E.21) 

- (F.AhAu'W/UA) - (F„)hAvw‘ /U r 1 ) = (f.mlifj 

where the expressions for F UU ,...F VW are given by equation (E.16). Subtraction of 
equation (E.21) from (E.18) and (E.20) from (E.19) yields the following combined 
equations: 

[{F„)e a - (F UU ) HA )(^/U r 2 ) + [(F vv ) Ea - (F VV ) HA ](^/Ur 2 ) 

+ [(E^u,)^ — (F WU! )h a ](w' 2 /U r 2 ) + [(F uv )e a — (F uv )H A ](u'v/U r 2 ) (e.22) 

— [(-Futu)^ ~ (Fuw)h a ]( u 'w' /U r^) — [(-Pvu,)^ — (F vw )h a ]{vw' jU T ) 

= [(-E mm )e a — (F mm )H A ] 

1(F uu )f a - (F„)c a ](^/UA) + {(F^)f a - (F„)o a \(^IVA) 

+ [(-Etuu>)F,i — {F ww )c A \{w >2 !U r 2 ) + [(•F , „„)j , x — {F uv )g a \{ u>v /F t ) (£.23) 

- [(F uw )f a - (F UW )G A }(^/Ur 2 ) - [(F vw )f a - (F VW ) GA )(^/U r 2 ) 

= [(f'mm)F^ — (F mm )G A ] 

Multiplication of equation (E.22) by [(Fuu,)^ — (-Futi^G*] an< l equation (E.23) 
by \(F uw ) Ea - (F uu ,)//^] and subtraction of the resulting equations yields: 

(S llItI - S 2u u)(^/U r 2 ) + (S lvv - S 2vv )(^/U r 2 ) 

+ (Siww - S 2w u,)(^/Ur 2 ) + (S luv - S 2uv )(u 7 v/U r 2 ) 

— ( 5" i mi- S 2vw ')(vw > /U T ) = (S\ mm S 2mm ) 


(£. 24 ) 
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where, 


S\ij [(-FuU))Fy| (•^UU')G y i][(-f 1 »i)Fyi 

$2ij = ~ (•Fuu))wx][(-^ , »i) f U — (A»j)ga] 


(i,j = u,v,w,m) 


By observing the behavior of the coefficients f u ,f v ,fw,fm,K i0 ,K u ,Ki 2 ,Ki 3 
and Ku at positions Ea,Fa,Ga and Ha, it can be shown that: 


and, 


(U 2 )e a = (fu 2 )H A 

Uv 2 )e a = Uv 2 )h a 

Uw 2 )e a = (fw 2 )H A 
( fm 2 )E A = ( /m )h a 

( fufv)E A = (fufv)H A 
(fuf w)e a = ( fufw)li A 

{fvfw)E A ~ ~{fvfw)H A 


(I<io)e a = (Kio)h a 
(K n )e a = -(Ku)h a 
(A'i 2 )e a = (Ku)h a 
(K\ 3 )e a = -( Ku)h a 
{Ku)e a = (Ku)h a 


; (/» 2 )f a = U 2 )g a 

; Uv 2 )f a = Uv 2 )g a 
; {Iw 2 )f a = {fw )gji 

; (fm 2 )F A = (fm 2 )G A 

j (/u/v)f.a = (fuf v)g a 
; (/u/u>)f4 = ~{fufw)Gj 

5 ( fvf\v)F A ~ (/v/u>)g> 

(Kiq)f a = (A'io)Ga 
(Ku)f a = -(Ku)g a 
(K 12 )f a = (A'i 2 )g a 
(Kis)f a = -(ATi 3 )g a 


(£.25. a) 


(£.25.6) 


(A' 14 )fa = (A"h)ga 

By imposing conditions (E.25.a) and (E.25.b) on equation (E.24), the following 
final form is obtained: 

(siK 12l E A ~ S 2 I<l2 i F A )(u t2 /U r 2 ) + (Si ~ S 2 )(v 2 /U r 2 ) + 

(siKi4,E A 5 2^14 ,F A ){ wt2 /Ur ) + 2(^i K\$ y E A “ 5 2^ v 10 ,F A ){ UV /Ur ) 

(‘S'lvui S 2 vw){.VW /Ur ) (^1 fm } E A $2 /iti,Fa)/2 

where 

= -8 [(/ u /u,)fa -(/mA 13 )FA(nFA +n G>l )](n Efl -n HB )(/m)EA 

52 = -8[(/ u / u ,)Fa -(/n>A'i3)EA(n£A +nFA)](nFA -nGA)(/m)FA 


Sivw = 16[(f u f w )F A - Kl 3 ,F A (frn)F A (n FA +Kg a )}x 

\(fvU)E A - I<n,E A {fm)E A (TlE A + A )} 

S 2 VW = 16 [(/u/u>)Ea - F i 3 t E A {fm)E A ( 1 ^E A + Hwa)] X 
[{fvfw)F A “A 'n,F A {fm)F A {^F A + n GA)l 
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and where the expression for II is given by equation (E.17). An algebraic ex- 
pression for vw' , applicable to low turbulence intensity flows, can be derived 
from equation (E.24) in a similar manner; however, it can also be extracted as 
a special case from equation (E.24) by substituting Fa = /,/,•; ( i = u, v, w), 

Fij = 2 fifj\ ( i,j = u,v,w ; t ^ j ) and F mm = f m 2 ir into equation (E.24), 

which yields: 


Sluu — Slw — — 0 

^2ttu = S>2vv = $2ww = Siuv = 0 


and,. 

Si vw = [2(fufw)F A *^(fufw)G A ][^(fvfw)E A ~ 2{fvfw)H A ] 

= (4fufxv)F A (4/ v /uj)Ea 
$2vw ~ i^fufw )e a {^fvfw^F A 

$1 mm ^ {^fufw)F A (/m )e a {^E a 

$2 mm = (4f u fw)E A (/m )f>* (^F* ) 

Substitution of these coefficients back into (E.24) yields, after some rearrange- 
ment: 


vw * _ \{fufw)F A {' K E A )e a (fufw)E A ( 7r F A ^G^jfm ) f ^] 

Ur 2 i^fufxv)F A i^fvfw)E A ""( 2 fufw)E A i^fvfw)F A 


(E. 27) 



U\ = ( U v + u r ) cos 9 — w* sin 6 
U 2 = {V + v) 

U3 = ( U r + n') sin 0 + w' cos 6 


Uff = U\ sin a + U 2 cos a 
Up = U\ cos a — U 2 sin a 



Fig. E.l. Effective cooling velocity for turbulence flow calculations. 



APPENDIX F 


VERIFICATION OF THE HOT-WIRE TECHNIQUE 

The low-intensity form of the Method A hot-wire technique proposed in 
the Appendix C has been partially verified by means of measurements in fully- 
developed turbulent pipe flow under simulated skewed flow conditions. In addi- 
tion, the implementation of the Method B technique into the present study was 
verified in pipe flow for the unskewed case. The pipe flow facility, as well as the 
full verification procedure for Method B , is described in detail in references [45] 
and [47]. 

F.l Fully-developed pipe flow 

A schematic of the pipe flow verification setup for Method A is shown in 
Fig. F.l. The probes are inserted into the flow through a slot (6.5 x 13 mm) 
located one diameter upstream from the pipe exit. The outside edges of the 
slot were chamfered to allow the probe to be skewed ±20 degrees in the pitch 
direction. The verification measurements were made on the far side of the pipe to 
minimize disturbances due to the slot. The diameter of the pipe is 10.16 cm and 
all the verification measurements were made at an operating Reynolds number of 
125,000 based on the centerline velocity and the pipe diameter. The development 
length for the pipe flow is L/D = 136. 

Along a radial traverse of turbulent pipe flow, azimuthal symmetry imposes 
the condition that the uw and the vw Reynolds stress components (relative to 
pipe flow coordinates) are identically zero. When the flow is fully- developed, the 
uv shear stress component, when normalized by the friction velocity, is a linear 
function of the radial location if the very near-wall region is excluded i.e.,: 



where the friction velocity is given as: 

u r /u b = vTTs (f.2) 

and U b (bulk velocity) and / (friction factor) are determined from: 

l/V7 = 2.01o glo (Re tN /7)-0.8 


with, 


Ub/Ud — RebjRtd 

Re cl = U c iD/v 



INTENTIONALLY KAMI 
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and where the Reynolds number ratio is given by a correlation proposed by 
Eppich [40]: 

Re b /Re c , = 0.045 log 10 (Re c/ ) + 0.600 

For the present study, the friction velocity was calculated from equation (F.2) 
and also measured independently with a Preston tube using the calibration given 
by equation (3.3). The results from the two methods were virtually the same: 

Ur/U b = 0.04726 @ Re ct = 125,000 


F.2 Verification approach 

Since, for low intensity flows, the yaw angle 7 (refer to Fig. C.2.a) can be 
determined with sufficient accuracy by nulling of a two-tube Conrad probe or by 
the normal hot-wire nulling technique described by Page [62], the low intensity 
verification procedure for Method A was initially performed for unknown pitch 
angles ( only. For this case, the transverse velocity component W is zero and 
the primary (mean flow) coordinates and secondary (turbulence) coordinates 
coincide (see Fig. C.3.a). Wire Positions 1,2,3 and 4 in the primary coordinate 
system (see Fig. C.2), then, are equivalent to wire Positions A , B, Ca and Da in 
the secondary coordinate system, respectively. For the following discussion, the 
probe rotation angle will be denoted by 9 and the prime (') superscript associated 
with the secondary (turbulence) coordinate system will be dropped. The probe 
rotation angle 9 is defined as shown in Fig. F.2. 

The proposed procedure for determining the mean velocity components and 
Reynolds stress components is presented in detail in Appendix C and is summa- 
rized here. The V/U and W/U mean velocity ratios are determined by applying 
data obtained with the offset prong slant-wire probe (Fig. C.l.b) at positions 
9 = 0 and ISO 0 and normal-wire data obtained at 9 = 225 and 315° to equation 
(C.14). Once the mean velocity ratios are known, the axial velocity component 
U is evaluated at wire Position 1 by means of equation (C.16), and it’s equiv- 
alent for wire Positions 2,3 and 4, and then the four values are averaged. The 
procedure for determining the components of the Reynolds stress tensor is as 
follows: 

1) u 2 ,v 2 ,uv - Data obtained with a normal-wire at 9 = 270° and with the offset 
prong slant-wire configuration (Fig. C.l.b) at 9 = 0 and 180°, when applied 
to equation_(C.36), results in a 3 x 3 system of equations which, when solved, 
yields the u 2 ,v 2 and uv Reynolds stress components. 

2) vw - Data obtained with the straight prong slant-wire probe (Fig. C.l.c) 
at 9 = 45,135,225 and 315°, when applied to equation (C.41), yields an 
algebraic expression for the vw stress component. 
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3) w 2 ,uw - Normal- wire data obtained at 0=225 and 315° and the stress com- 
ponents from 1) and 2) above, when applied to equatio n (C .36), results in a 
2x2 system of equations which, when solved, yields the w 2 and uw Reynolds 
stress components. 

By considering pitch skewness only in the fully-developed pipe flow verification, 
the uw and vw stress components are identically zero regardless of the pitch angle. 
Further, the w 2 stress component is invariant as the pitch angle is changed. 

F.3 Hot-wire probe design and construction 

The proposed Method A hot-wire technique requires the sequential use of 
two and possibly three probe configurations which are illustrated in Fig. C.l. 
All the probes for the Method A technique were constructed from scratch with 
the basic dimensions and materials as shown in Fig. F.3. Detail dimensions of 
the probe prongs for the three probe configurations are shown in Fig. F.4. One 
of the constraints imposed on the probe design was that it had to fit through a 
4.57 mm (0.180 in) access hole which corresponds to the diameter of a standard 
Thermal Systems Inc. (TSI) probe holder. This constraint places a limitation on 
the amount of offset that can be applied to a probe such as the slant- wire probe 
shown in Fig. F.4.b. Given a diameter that the probe must pass through without 
interference, a compromise must be made between probe offset and the distance 
between the probe prong tips (hot-wire length). 


F.4 Calibration 

The hot-wire probes were calibrated in the core flow at Data Station 1 of 
the transition duct flow facility (see Fig. 3.2). A Pitot tube placed adjacent 
to the hot-wire probe and the local wall static pressure provided a reference 
calibration velocity. The normal-wire was calibrated in the 0 = 270° position 
and the straight prong slant-wire was calibrated in the 0 = 0° position. The 
offset prong slant- wire was calibrated in both the 0 = 0 and 180° positions. A 
typical plot of an offset prong slant-wire probe calibration is shown in Fig. F.5. 
Clearly, the calibrations are not equivalent. There are two possible reasons for 
the differences. First, unlike Method B where, for unskewed flow, the stagnation 
line for the normal velocity component on the hot-wire is unchanged regardless of 
rotational position, for Method A the stagnation line rotates around the hot-wire 
as the probe is rotated. Waviness or other imperfections in the wire may cause it 
to respond differently when the flow is exposed to the front or back “side” of the 
wire. The second cause is the influence of the upstream probe prong when the 
wire is in the 0 = 180° position. There are indications that the latter cause is 
responsible for the discrepancy. If wire imperfections were responsible, then the 
calibration curves should vary randomly with each different wire. Calibration 
with separate but similar probes and several different wires has shown that the 
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calibration curve for Position B (6 — 180°) consistently rides slightly higher then 
the curve for Position A. If it is assumed that the wire is exposed to the true 
c ali bration velocity when in Position A and that the wire itself is perfect, that 
is, it will respond the same in any orientation so long as the effective cooling 
velocity remains constant, then Fig. F.5 indicates that the wire is exposed to a 
higher velocity when in Position B than when in Position A. It is plausible that 
the higher velocity in Position B is a result of the local acceleration of the flow 
as it is displaced about the upstream probe prong. This behavior is illustrated 
in Fig. F.6. 

The question arises as to whether the prong influence is significant enough to 
preclude accurate flow measurements. The effect of different calibration slopes in 
the 6 = 0 and 180° positions can be taken into account in the mean flow response 
equations by replacing equation (C.7.a) by: 

E 12 =(B ! /B 1 ) 2/m (Si/S 2 ) (FI) 

and by using the appropriate value of B in equation (C.16) or its equivalent for 
wire Position 2. The turbulence response equations Me unchanged since they 
are cast into a non-dimensional form in which the calibration slope is eliminated 
(see equation (C.25)). By accounting for the different calibration characteris- 
tics of the offset slant-wire, it should be possible to obtain accurate mean flow 
and turbulence data for unskewed flow conditions (V=W=0). The effect of the 
probe prong interference under skewed conditions is uncertain and requires fur- 
ther study. 

F.5 Unskewed flow verification 

The mean velocity and Reynolds stress components in the fully-developed 
pipe flow were measured by means of the Method A and Method B hot-wire 
techniques. The measurements by the Method B technique were performed to 
demonstrate proper operation of the data acquisition equipment and correct pro- 
gramming of the response equations. For the Method A measurements, the offset 
slant-wire data were reduced using the calibration in Position A only (uncor- 
rected) and also using the calibration in Positions A and B (corrected) as de- 
scribed in the previous section. The measured axial mean velocity profile and 
transverse flow angles are shown in Fig. F.7.a and Fig. F.7.b, respectively. In 
Fig. F.7.a, a reference velocity distribution measured with a flattened Pitot probe 
is shown for comparison. The axial velocity results for Method A were calcu- 
lated from the output of the slant- wire only, rather then averaging the slant and 
normal-wire results, in order to illustrate more clearly the effect of the calibra- 
tion anomaly. The corrected Method A results and the Method B results show 
very good agreement with the reference distributions. The uncorrected Method 
A results show the need to include the calibration of the offset prong slant-wire 
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in both Positions A and B. All subsequent Method A results will be based on 
calibrations in both positions. 

The six components of the Reynolds stress tensor for unskewed conditions 
were measured by means of both the Method A and Method B techniques. The 
normal stress and shear stress profiles are shown in Fig. F.8.a and F.8.b, re- 
spectively. The normal stress distributions are compared to the fully-developed 
pipe flow results due to Lawn [64], which were taken at a comparable Reynolds 
number. The shear stresses are compared with the theoretical values. Overall, 
the results for both techniques show good agreement with the reference distribu- 
tions. It should be noted here that agreement of the measured uw and vrU shear 
stresses with the zero theoretical value is really only a partial verification. This 
is due to the fact that when the yaw angle is zero, the response equations will 
yield zero values for these components, regardless of whether or not there are 
probe disturbances present, so long as the probe output exhibits rotational sym- 
metry about the 8 = 0° position. Since the normal-wire probe is used to measure 
the u> 2 stress component as well as the uw component, there is justification to 
assume that interference effects associated with the normal-wire probe are mini- 
mal, inasmuch as the tu 2 stress component is predicted well. The straight prong 
slant-wire, however, is used solely to measure the vw stress component. Full 
verification of measurements with this probe must be made by simulating yaw 
and pitch flow conditions. 

The unskewed flow case is used to determine the reference tangential cooling 
coefficient fc re / for the hot-wires. For both Method A and Method B , k re f for 
the slant-wires is varied until the best agreement with the theoretical uv shear 
stress distribution is obtained. This is accomplished by fitting a straight line 
(least-squares fit) through the uv/U % results and adjusting k re f until a slope of 
one is obtained. To avoid probe prong disturbances, the straight prong slant-wire 
of Method A had to be used in the Method B mode to determine its value of k re f. 
For the normal-wire used in Method A , k rt j is adjusted until the best agreement 
with the Method B or Lawn’s w 2 distribution is obtained. Specification of k re j 
for the normal-wire used in Method B is not critical, inasmuch as the tangential 
velocity component is always much smaller than the normal component for zero- 
to-moderate skewness levels. Typical values for k re f were within the range of 
0.05-0.2. 

F.6 Skewed flow verification 

With good results achieved for the unskewed case, the next step was to 
repeat the measurements under simulated pitch skewness conditions. Verification 
measurements were obtained at pitch angles of ( = -20,-10,0,-1-10, +20°. The vw 
stress component for the skewed conditions was not measured because of the 
aformentioned problem of verification under simulated pitch skewness only. For 
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calculation of the w 2 and uw components (step 3 in Section F.2), the vw stress 
component was set to its theoretical value of zero. 

Measurements of the mean flow field under simulated pitch conditions are 
shown in Fig. F.9. The axial velocity profiles shown in Fig. F.9.a have been 
transformed back to the unskewed (pipe flow) coordinate system to facilitate 
comparison with the Pitot tube results. Overall the agreement is fairly good. 
The spread in the data increases slightly as the pipe wall is approached (r/R = 
1.0) and is probably related to the higher turbulence intensity in that region. The 
measured transverse flow angles are shown in Figs. F.9.b and F.9.c. Generally 
the technique predicts the set flow angles (solid-line distributions) well, but slight 
deviations as the wall is approached are again evident. 

Five components of the Reynolds stress tensor were measured (the vw com- 
ponent was not measured) using the Method A technique under simulated pitch 
skewness conditions. The results are shown transformed to the unskewed co- 
ordinate system. The normal stresses are compared with Lawn’s data and the 
shear stresses are compared to the theoretical distribution. The, results of the 
turbulence measurements are shown in Fig. F.10. These results show that good 
agreement is obtained near the centerline of the pipe, but that systematic de- 
viations occur as the wall is approached. Some conclusions can be drawn from 
the manner in which the stresses are calculated. Since the uw and vw stress 
components are zero for all of these measurements and the coordinate rotation 
is about the 2 -axis, the only stresses which are affected by the transformation 
are the u 2 ,v 2 and uv components, which are computed independently of the 
remaining stresses (see step 1 in Section F.2). The output from the normal-wire 
probe in the 6 = 270° position (from which the streamwise stress component 
is primarily derived) was observed to be the same for all pitch angles. This is 
expected since probe prong disturbances are unlikely for this configuration and 
the effective cooling velocity is invariant with pitch angle. This implies that the 
deviations are attributable to the offset prong slant-wire probe. This, of course, 
shouldn’t be surprising since a prong interference effect was clearly present in 
the calibration results. The w 2 stress component deviates from Lawn’s data 
because it’s calculation is coupled with the it 2 ,v 2 and uv components. It was 
considered that perhaps the deviations from the reference distributions for the 
skewed flow case were due to an error either in the derivation of the response 
equations or in the implementation of the equations into a computer program. 
However, since no errors have been detected and the fact that the results agree 
well near the center of the pipe, it is believed that the progressive deviation as 
the wall is approached is a result of the increasing turbulence intensity, coupled 
with the probe prong interference effects. 
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F.7 Conclusions and recommendations 

A hot-wire technique for measuring the mean flow and Reynolds stresses 
in a low turbulence intensity flow with moderate mean flow skewness has been 
proposed. Preliminary verification measurements have been obtained under sim- 
ulated pitch angle conditions. The results are encouraging but further work is 
needed. In particular, the technique is adequate for determination of the mean 
velocity components for pitch angles up to ±20 degrees. Indeed, the results in- 
dicate that this range can probably be extended. The technique is also adequate 
for determination of the Reynolds stress tensor for small skew angles (less than 
±10 degrees), but systematic deviations were observed for moderate skew levels. 

It has been concluded that the source of errors is attributable to the design 
of the offset slant-wire probe. Specifically, the hot-wire is apparently sensing the 
local acceleration of the flow around the longer probe prong when it is positioned 
upstream of the center of the hot-wire. The effect of the prong disturbance has 
been at least partially compensated for by calibrating the wire in two positions. 
In order to reduce the effect of prong disturbances, a new offset slant- wire probe 
design is proposed here. This new probe configuration, which is illustrated in 
Fig. F.ll, has the following advantages: 

1) The offset distance is greater while still remaining within the clearance en- 
velope. 

2) In the rotational positions 9 = 0 and 180°, the probe prongs are not upstream 
of the sensing element. Even if some prong influence does occur, it should be 
the same in both positions, thereby eliminating the need for two calibrations. 

3) By closely spacing the prongs and filling between them with epoxy, a more 
aerodynamic shape is achieved which will reduce the chance of prong vibra- 
tion due to vortex shedding. 

4) The relocation of the longer prong changes the rotational sweep where prong 
interference will occur. An estimate of the interference region is shown in Fig. 
F.12. The need for the straight prong slant-wire probe may be eliminated. 

The primary disadvantage of this probe is that it is considerably more difficult 
to construct. 

It is felt that further refinement of this technique is justified inasmuch as 
there are many flow configurations where it is not practical to align the rotational 
axis of a hot-wire probe with the direction of the primary flow. The simple cali- 
bration and the increased sensitivity make this technique an attractive alternate 
to the technique proposed by De Grande and Kool [61] . Should the new probe 
design yield improved results, the verification procedure should be expanded to 
investigate the best way to obtain the secondary shear stress vw . Once accom- 
plished, a study of the performance of the technique for unknown yaw skewness 
and moderate turbulence intensity flows should be undertaken. 
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Fig. F.l. Schematic of Method A hot-wire verification setup. 




Fig. F.2. Hot-wire orientation for Method A. 
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All dimensions in millimeters 
4.57 dia clearance envelope 



a) Normal-wire. 

Fig. F.4. Detail of probe prong geometry for Method A. (Continued . . . ) 
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All dimensions in millimeters 


4.57 dia clearance envelope 



b) Offset prong slant-wire. 



All dimensions in millimeters 
4.57 dia clearance envelope 




c) Straight prong slant-wire. 


Fig. F.4. Detail of probe prong geometry for Method A. 
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Fig. F.5. Offset slant-wire calibration in Position A and Position B. 



Fig. F.6. Conjectured probe prong influence for offset slant-wire. 
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a) Axial velocity profile (from slant-wire only). 




b) Transverse flow angles. 

Fig. F.7. Fully-developed mean flow for unskewed conditions. 
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a) Normal stresses. 

Fig. F.8. Fully-developed Reynolds stress profiles 
for unskewed conditions. (Continued . . . ) 
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r/R 


b) Shear stresses. 


Fig. F.8. Fully-developed Reynolds stress profiles for unskewed conditions. 
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a) Transformed axial velocity profile (from slant-wire only). 
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r/R 

b) Pitch angle. 



c) Yaw angle. 

Fig. F.9. Fully-developed mean flow for skewed conditions. 
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a) u' stress component. 

Fig. F.10. TVansformed Reynolds stress profiles. (Continued . . . ) 
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b) v' stress component. 

Fig. F.10. Transformed Reynolds stress profiles. (Continued . . . ) 
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c) w* stress component. 

Fig. F.10. Transformed Reynolds stress profiles. (Continued . . . ) 
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e) uw stress component. 

Fig. F.10. Transformed Reynolds stress profiles. 
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All dimensions in millimeters 


Fig. F.ll. Proposed offset prong slant-wire probe. 



Fig. F.12. Interference region for proposed offset prong slant- wire probe. 



APPENDIX G 


UNCERTAINTY ESTIMATES 

Data presented in this dissertation axe deduced from measurements with 
surface static pressure taps, static pressure probes, total pressure probes and 
hot-wire anemometry. The uncertainly estimates quoted for a particular variable 
represent error bounds associated with the peak value observed for that variable. 

G.l Pressure Measurements 

Uncertainties associated with pressure probe design are based on data pre- 
sented by Chue [41] Static pressure coefficient ( C p ) values based on wall tap 
readings and the static pressure probe are estimated to be accurate to within 
±0.002. Skin friction coefficient (C/) values deduced from Preston tube results 
at stations 1,5 and 6 are estimated to be accurate to within ±0.0001, ±0.0005 
and ±0.0002, respectively. Normalized velocity components U + = U/U T based 
on total pressure probe results at stations 1 ,5 and 6 are estimated to be accurate 
to within ±0.2, ±0.7 and ±0.3, respectively. 

G.2 Hot-Wire Measurements 

The uncertainties in hot-wire based results were estimated from an error 
analysis presented by Al-Beirutty [45]. In that analysis, the effects of uncertain- 
ties in the slant- wire angle a, binomial and tangential cooling coefficients ft bn 
and ft, normal and slant- wire calibration intercept values E 0i n and Eq iS , slant- 
wire calibration slope B c , and the measured normal and slant-wire voltages E n , 
Eg, e 2 and e 2 are taken into account. 

If the above considerations are applied to the present results, then the uncer- 
tainties in hot-wire based values of U jUy, and V r /Ub are estimated as ±0.01 and 
±0.002, respectively. Uncertainties in the Reynolds stress components u 2 /U ^ , 
w 2 /U 2 , uv/Ul , uw/Ul and vw/U% are estimated as ±0.0001, ±0.0002, 
±0.0002, ±0.00015, ±0.00015 and ±0.0001, respectively. 
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